

v-^TUDIES OF HEAT TRANSFER IN THE FIRST 


AND THE SECOND ORDER BOUNDARY 

LAYER FLOWS 


BY 

NOOR AFZAL 


3 > 



9/p 

DEPARTMENT OF AERONAUTICAL ENGINEERING 

INDIAN INSTITUTE OF TECHNOLOGY KANPUR 




1969 



^rnniES OF HEAT TRANSFER IN THE FIRST 
and THE second ORDER BOUNDARY 

^ A ’.rr-'T^ 'Cl n \\ 7 ^ 


In 


A Thesis Submitted 
Partial Fulfilment of the Requirements 
for the Degree of 


doctor of philosophy 


BY 

noor afzal 


- C i? 

fK ft 


f 

i 







to the 


department of aeronautical 

INDIAN INSTITUTE OF TECHNOLOGY KANPUR 


1969 



I.!.T. KANPUR 

C&MTHAL LIBRARY 


. If -Jftil 


W A f CIS 




CERTIFICAT E 


This is to certify that the present work STUDIES 
OP HEAT TMSFER ID THE FIRST iJJD THE SECOND ORDER BOTTUD/.RY 
LAH'R Flows has been carried out under my supervision and 


ha.s not been submitted elsewhere for the aAvard of a degree 


or diploma. 


Madan Mohan Oberai 

Assistant Professor 

Dept, of Mechanical Engineerinp; 

Indian Institute of Technology 

Kanpur, India, 



ACKNOVyLEDGEMBWT 


It is a great pleasure to thank Dr, Madan Mohan Oherai 
for his invaluable guidance, comments, criticisms and encourage- 
ment throughout the course of present work, I also appreciate 
his liberal attitude towards the freedom of thought which I 
enjoyed during his guidance. 

I also wish to thank Dr, Roddam Farasimha for introduc- 
ing me to this subject and for many stimulating discussions 
during his brief visit to Kanpur,. I also take this opportunity 
to thank all the members of ’ Aerod'naroic Group' particularly 
Dr- A. C. Jain for many interesting discussions,. 

further, I acknowledge with gratitude the continuous 
interest end encouragement of Dr. P.N, Murthy and Dr. S.M. 
Ramachandra during the course of this work. 

I am Gulso grateful to Mrs. Rachel Oberai for her 

generous help. finally, thanks a.re due to Mr. J.U. Siddiqui, 

Mr. M.Z. ■■ Siddiqui and Mr. B.K. Gupta for their interest in the 

work and due to Mr. J. K. Misra of Graphic Arts for his coopera- 

tion in bringing out the work in the present form. 


No or Afzal 



TiLBLl. OF CONTENTS 


Pa.a:e No. 


NOIvJENGLATUNE vi 

LI^T OE EIGUPES xiii 

SYNOPSIS xvii 

CHAPI'ER 1 INTRODUCTION 1 

CHAPTER 2 EORiroLATION OE THE SECOND ORDER BOUNDARY 

LAYER THEORY 

2.1 Introduction 17 

2.2 Coordinate System 18 

2.3 Governing Equations 19 

2.4 Eormulation 21 

2.5 Governing Equations for the First and 

the Second Order Boundary Layer Problems 32 

2.5.1 First Order Boundary Layer Problem 32 

2.5.2 Second Order Boundary Layer Problem 53 


CHAPTER 3 FIRST ORDER COMPRESSIBLE: BOUNDARY LAYER 
ON A FLAT PLATE WTTH ARBITRARY THERJ/IAL 
PROPERTIES 


3. 1 Introduction 35 

3.2 Governing Equations 37 

3.3 Solution at Low Prandtl Numbers 39 

3.4 Solution at High Prandtl Numbers 50 

5.5 Solution for Prandtl Number of Order Unity 55 

5.6 Discussion 60 



CHAPTER 4 HEAT TRANSEER IN PALKIffi'R-SKAN BOUNDARY 
LAYER ILOWS 

4-1 Introduction 63 

4.2 Gcverning Equations 64 

4.3 Solution of Palkner-Skan Equation 66 

4.4 Solution at High Prandtl Numher 69 

4-5 Solution at Low Prandtl Numher 80 

4.6 Solution for Prandtl Numher of Order Unity 84 

4.7 Discussion 87 

CHAPTER 5 SECOND ORDER EFFECTS ON LAMINAR BOUNDARY 
LAYER IN NON-SIMILAR FLOWS OF A CONSTANT 
PROPERTY FLUID 

5.1 Introduction 92 

5.2 Governing Equations 93 

5.3 Boundary Layer Characteristics 95 

5.4 Choice of Intiependent Variables 97 

5.5 First Order Boundary Layer 101 

5.6 Second Order Boundary Layer Theory 110 

5 . 6.1 Transverse Curvature Problem 111 

5 . 6.2 Longitudinal Curvature Problem 119 

5.6.3 External Vorticity 124 

5 . 6.4 Stagnation Enthalpy Gradients 128 

5 . 6.5 Displacement Speed Problem 129 

5.7 Discussion 155 

REFERENCES 143 

APPFiNDIX .A EVALUTION OF INTEGRALS 153 

APPENDIX B EULERISATI0N:(E, q) METHOD 156 

APPENDIX C NUMERICAL SOLUTION OF HIGHER ORDER BOUNDARY 
. , BOUNDARY LAYER PROBLEMS I 60 



NOaffiNGLATURE; 


a.,TD,w, 
B ( yj ), 

C, 


d , 

D( p , 0“ ) -, 



indices of power laws defined by equation (3.7) 

Blasius solution for flat, pla^te 

Bernoulli function defined by equation (2.11a) 

compressibility factor 

skin friction coefficient 

specific heat at constant pressure 

a constant defined as = b(l + m) 

heat transfer rate at the wall defined by the 

equation (4.26a) 

a constant defined as = _(8+w)(1+ni) 

Eckert number defined as = 

dimensionlesr velocity in s direction due to 
first Order boundary layer 

correction to f due to second order boundary layer 


2 ^-^ 5 hypergeometric function defined as 

2]?pa,b,c,x) =E (a)n(t)„^"/[(c)„ «'■] 
n — 0 _ . j — -tt 

where (a)^ = j ( a+nj / j(a) 

gj dimensionless first order bounda^ry layer entha.lpy 

G, change in g due to second order boundary layer 

h , enthalpy of the fluid 

h^^jhg) .. first and second order enthalpy in the inner 
(boundary layer) flow defined in (2.17) 



vii 


), 
3 , 

k, 


first and second order entha.lpy in the outer flow 
defined in (2.7) 

HO f O 

a function defined as = f + pf' (1-f ). 

a number equal to zero for two dimensional flow 
and unity for axisymmetric flow 
thermal conductivity of the fluid 


k 


1 ’ 



K, 

L, 


m 


> 


M ,M , 
r’ rr’ 

n. 


n ( p ) V 
N , 

N ( a 5 X ) 


N^(a,x) , 

P, 


Pl,P2»- 


longitudinal curvature parameter defined as = 

K f(2^ )/(U^r^) 

transverse curvature parameter defined as 
= (j cos 0/r) \f(2^ )/(U^r^) 
longitudinal surface curvature of the body 
characteristic length of the body 

operator defined as = u^ 6/5s + v^ 6/3H', r = 1,2,.. 
c'- constant defined as = (b-a-w)/( 1 -b+a+w) 
functions defined by equations (4.23) and (4.34) 
coordinate normal to body 
a function defined 3.s=Il^(p, °°) 

1 /2 

boundary layer (inner) variable defined as =R ^ n 
related gamma function defined as 

= a at 

° ^X 

function defined as = / t™14(a, t) dt 

0 

pressure of fluid 

first and second order pressures in the outer flow 
defined in (2.7) 

first and second order pressures in the inner flow 
defined in (2. 17) 



< 1 , 


the heat transfer rate at the wall 


viii 

fi^rst and second order heat transfer at the wall 
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Q( )? differential operator defined as 
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SYNOPSIS 


The present work deals with the study of the second order 
effects to the classical boundary layer theory when characteristic 
Reynolds nuiriber of the flow is only moderately large. However, 
in the beginning, we have studied some aspects of classical 
boundary la.yer theory itself, which have either not received 
sufficient attention in the literature or the treatment is not 
exhaustive or coherent. 

The second order effects studied in the present work are 
due to transverse curvature, longitudinal curvature, external 
vorticity, stagnation enthcolpy gradient and displacement, which 
are described by linear equations. For each of these second 
order effects a general theory of Gortler type is developed. 

The flow ma.y be two dimensional or aixisymraetric. Attention is 
focused on low speed and small temperature changes, so that 
the velocity field is es'"-‘ent ially that of an incomnres^ible 
fluid. We calculate temperature profile too, dissipation is 
neglected. 

In addition to principal velocity and thermal functions 
a principa-1 shape f-unction is defined for axisymmetric flows 
in the study of first order boundary layer and a new principal 
function is introduced for each of the second order effects ; 

(i.e., principal transwerse curvature function, principal ; 

longitudinal curvature function, principal vorticity interaction; 
function and principal displacement speed function). These new i 
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principa,! functions are of fundamental importance in determin- 
ing the structure of the second order effects. In the present 
work we have studied the first two term.s in the second order 
Gdrtler power series for the velocity and temperature profiles, 
hy choosing the proper universal functions. The resulting sys- 
tem of total differential equations and the accompanying boundary 
conditions constitute the two point boundary value problems and 
are solved hy Runge-Kutta. method with Gill improvement on IBM 
704-4 Computer at Indian Institute of Technology Kanpur, Solutions 
arG obtained for both accelerating and decelerating flows with 
Balkner-Skan pressure gradient parameter |3^ between -0.198838 
and 2,0. and the Prandtl number 0.7 to 5. The classical concepts 
of displacement and momentum thicknesses are extended to higher 
order boundary layer flows. Second order perturbation to dis- 
placement and momentum thicknesses are evaluated separately for 
each of the above mentioned second order effects. The results 
for the variation of various second order quantities such as 
changes in velocity profile, temperature profile, skin friction, 
heat transfer, displacem-ent thickness and momentum, thickness due 
to various second order effects are represented graphically. A 
critical comparison with various available resu'Lts is presented. 

Our solutions for the second order effects sho?(;' that the 
transverse curvature increases the skin friction and heat trans- 
fer. The convex longitudinal surface curvature decreases the skin 
friction and heat transfer. The effect of external vorticity 
is to decrease the skin friction. The heat trans-fer is a.lso 



decreased except for > 1 where it increases. The effect of 
stagnation enthalpy gradient is to decrease heat transfer, Finally 
the displacement speed increases the shin friction and heat trans- 
fer. Lastly, the results indicate that the convergence of Gortler 
power series for the first and the second order "boundary layer 
functions "becomes poor as favourable pressure gradient diminishes. 
Furthermore, in adverse pressure gradient (near = -0.198838) 
some of the second order boundary layer quantities beoomes very 
large compared with the corresponding first order boundary layer 
quantities. This indicates that the boundary layer (inner) ex- 
pansion is not uniformly valid near the separation.. 

The other two problems studied in the beginning of present 
work arc of first oider boundary layer theory. First one is the 
study of compressible boundary layer with no pressure gradients^ 
with arbitrary thermal properties (usina the power la,w relations 
for the dependence of viscosity, density and Prandtl number upon 
temperature) separately for small, In-ge and order unity values 
of Prandtl number. For limiting cases of low and high Prandtl 
number, the method of matched asjrniptotic expansions is used. 

In low Prandtl number problem it turns out that the outer flow 
is I governed by a simple nonlinear equation of second order, 
which is correct to all orders in Prandtl number and can be solved 
once and for all. For the case of insulated wall this cuter equa- 
tion can be solved a,na.lyticailly. However, for heat transfer case 
it is solved numerically (to yield an outer solution which is also 



correct to o21 orders in Prandtl number) by Runge-Kutta method. 

This numerical solution is matched asymptotically with the corres- 
ponding inner solution. For high Prandtl number case, the inner 
flow is governed by coupled nonlinear differential equations, with 
the outer boimdary conditions to be determined by matching it 
with outer solution (which is singular) and, in general, is very 
difficult to solve exactly. In this case an approximate solution 
of the ?/eyl type is carried out. For Prandtl number of order 
unity, the Meksyn type method of inversion of variables is used. 
Further it is shown that if Prandtl number is constant, equal to 
xinity the classical Croccos integral holds even with arbitrary 
properties. In all three cases of low, high and order unity values 
of Prandtl number, the leading terms are obtained for recovery 
factor and heat transfer rate at the wall. The results a,re dis- 
played graphically. 

The second problem in the study of first order boxmdary 
layer is the study of heat transfer in Falkner-5ka..n flows for 
small, large and order unity values of Prandtl number (dissipation 
is included in the energy equation). For low and high Prandtl 
numbers, the method oF , matched asymptotic expansions is used 
while for Prendtl number of order unity a Meksyn type method is 
used. The probleni has been studied extensively in the past and 
most of the authors have either ignored the presaure gra.dient or 
dissipation or both. 
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In luw and high Prandtl mxmher flows, many higher order 
perturbations are obtained and results are good even up to 
Prandtl number of order unity. Detailed examination of the 
results indicate that in favourable pressure gradient the 
results of both the analyses put together fill the gap in the 
intermediate region where Prandtl number is of order unity. In 
adverse pressure gradient, there is a gap in the intermediate 
regi&n 0.1 < CT < 1 due to poor convergence of both the series 
in the above region. However, if the series are properly 
Eulerised results fill the gap. 

Thus we see that in some sense it suffices to analyse 
the two limiting cases: of low and high Prandtl numbers, .pro- 
vided large number of terms are computed in each expansion 
and series are properly Eulerised. 



CHAPTER 1 


INT RODUCTION 

The classical houndary layer theory of Prandtl represents 
the leading term in an asymptotic expansion of the Navier-Stokes 
eoLuations for high Reynolds number R, which becomes asymptotically 
exact only for very large Reynolds number. Neveig-theless, these 
asymptotic solutions are exceedingly useful from a practical point 
of view owing to the fact that in many cases these have been found 
to hold with surprising accuracy even under distinctly nonasymptotic 
conditions. It appears desirable therefore to study a general tech- 
nique for constructing the corrections to these asymptotic solutions 
in order to extend their range of validity. The practical utility 
of such a refinement is suggested by an estimate of lagerstrom and 
Cole (1955) that the second approximation may in some cases predict 
the resiilts accurately down to Reynolds number ten or even five. 
Moreover, it is shown here that in some cases (see Chapter 4) it 
suffices in some sense to study the two limiting cases of very 
small and very large values of the parameter, as the results obtained 
from above two cases overlap in the intermediate range, where 
parameter is of order unity. 

further, the Prandtl boiindary layer theory does not account 
for what Rott and lenard (1959) call the secondary effects. These 
effects may be classified as those due to longitudinal curvature, 
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transverse curvature, external vorticity gradients, stagnation 
enthalpy' gradients, displacement, slin and temperature jump, At 
lower Rejmolds numhers, however, these various (above mentioned) 
effects are of the order and necessiate corrections to the 

classical boundary layer theory. However, for flow situations in 
which these effects are of order unity (for example if is 

of order unity which implies E the longitudinal curvature is of 
the order as R becomes very large), some authors have modi- 

fied the tangential momentum equation of thp Prandtl so as to 
include such terms but have ignored the fact that a correspond- 
ing change becomes necessary in the normal momentum equation too — 
3p/8n will now he of the order of and hence the surface pres- 

sure cannot be equated to its value at the edge of boundary layer. 
This statement may appear too strong, but we feel that under these 
circumstances the treatment should he called empirical rather than 
one based on Prandtl’s boundary layer theory. Furthermore, these 
authors have studied only the individual effects by formulating 
a nonlinear total differential equation. Thus in a problem deal- 
ing with various effects simultaneously we cannot use the results 
of the above authors as the superposition is not valid. In the 
present work only those flows will be studied in which these 
effects are really secondary and appear only in the form of 
corrections (described by linear equations) to the classical 
fCfirst order) boundary layer. 

I f, ■ ■ ' 
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A systematic formulation of the scheme to calculate the 
higher order corrections to classical houndary layer theory is 
giTen by Van Dyke (1961, 196 2) for compressible and constant 
properly fluids respectively'-. The method used is the techniq.ue 
of matched asymptotic expansions developed "bj Kaplun (19545 1967), 
Lagerstrom and Cole (1955), etc, 1 feature of the perturbation 
technique is that the resulting second (and higher) order boundar 3 ’^ 
layer (inner) and outer equations are linear. Following the sugges- 
tion of Rott and lenard (1959), we can subdivide these higher order 
problems into a number of simpler problems, each of which has a 
clear physical interpretation and then superpose. 

The second order boundary layer theory' developed for com- 
pressible flow of a fluid by Van Dyke (1961), is used by him to 
discuss the flow past an axisymmetric stagnation point. He uses 
the first term of the Blasius series expansion in the coordinate 
along the body and finds the solution of the first and second order 
boundary layer equations in the stagnation region with viscosity 
proportional to temperature, Prandtl number is 0.7 and the ratio 
of wall to stagnation tdnperature is 0,2. Davis and FlUgge-Lotz 
(1964a) have further calculated the vorticity effect (due to en- 
tropy gradients only) for various values of the ratio of wall to 
stagnation temperature when viscosity varies as the square root 
of temperature and the Prandtl number is 0,71. JiSe Davis and 
PlUgge-Dotz (1964b) have solved the partial differential equations 
of the second order boundary layer theory with the aid of finite 
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difference method for studying the second order effects (except 
stagnation enthalpy gradients) for hypersonic flow past a para- 
boloidal and hyperboloidal at free stream Mach number of infinity, 
and a sphere at free stream Mach number of ten, with Prandtl number 
0,71. The corresponding problero.s of boundary layer around a cir- 
cular cylinder and slab with circular leading edge have been 
analysed by Pannelbp and PlUg.ffe-Lotz (1966). 

for the case of constant propeity flow Van Dyke (1962b) has 
studied the flows past the stagnation point of axts 3 nmrietric or 
plane body and also of cusped leading edges at ideal incidences. 

He also studied the heat transfer with Prandtl number 0.7 and 
1.0, later Van Dyke (1964) applied the above theory to a para- 
bolic cylinder in uniform stream. Devan and Oberai (1964) for- 
mulated an integral relation of the Karman type for the second 
order momentum boundary layer, the flow being two dimensional, 
later Devan (1965) formulated the integral relation in Crocco 
form (for the above problem) and applied the Dorodnitzyn method 
of approximate solutions. 

In addition to Van Dyke, the second order boundary layer 
theories have been given by lenard (1961) and Malsen (1963). These 
authors have provided the solutions which are valid in the stagna- 
tion region only. In general, they are in agreement with Van Dyke 
in their results with the exception of the solution for vorticity 
interaetion. Honda and Kiyokowa (1969) have also developed a sys- 
tematic second order boundary layer theory for the flow of an 
incompressible fluid in stream function coordinates. 
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,.A critical examination of the various second order and 
related theories and comparison with experiments up to that time 
have been discussed rather thoroughly by Van Dyke (1962c). Cheng 
(1965) has also given a good account of various second order 
theories. Eurthermore, each- of the secondary effect has been 
dealt with by several authors. The results of most of them- are 
in disagreement with each other. One source of confusion was the 
lack of realization that all the second order effects are not 
independent in the strict sense; the difficulties resulting from 
the arbitrary division of these effects can be avoided only by 
treating them concurrently. However, whatever the division, the 
effects relegated to the category of displacement are global in 
nature (and as such are most . difficult to calculate) while the 
others are local. Second source of confusion was that the analy- 
sis of most of the, authors was not systematic upto the second 
oidler as they have omitted some of terms of the same oi^er as those 
retained. Thirdly, in some cases the matching of the inner and 
outer flow has not been carried out correctly. 

The subject of longidutinal curvature has been studied by 
Tani (1949, 1954), Murphy (1953, 1962, 1965), Yen and Toba (1961, 
196 2), Hayasi (1963), Massey and Glaj^on ( 1965 ), Schultz-Grunow 
and Breuer (1965) and Narasimha and Ojha (1967). The results 
previous to Van Dyke are controvertial. Thus Tani and Murphy 
concluded that.- the skin friction will decrease on the convex side 
because of the longitudinal curvature, Yen and Toba arrive at just 
the opposite conclusion. However, the formulations of Murphy, 
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Yen ana Toba and Hayasl are not consistent up to the order consi- 
dered as these authors have omitted some of the teimiB of the second 
order. Sohultz-Srunow and Breuer studied the effect of longitu- 
dinal curvature with no pressure gradient by formulating a non- 
linear differential eauatlon. Their equation contains all terms 
of order unity and yan Dyke's equation and Include 

in addition one viscous term of order R'^. The above analysis is 
not systematic in the sense that using their formulation one can 
not obtain hj^her order approximations. Similar comments apply 
to the work of Massey and Clayton who also constructed a nonlinear 
equation. Harasimha and Ojha using Van Dyk^s formulation studied 
only the similarity term for the effect of longitudinal curvature. 

The transverse curvature problem for. a longitudinal flow 
of an incompressible fluid over a circular cylinder has been studied 
by Seban and Bond (1951). Clauert and Lighthill (1955), Stewartson 
(1955) and Jaffe and Okamura (1968). Bor compressible flow of a 
fluid past a slender body of revolution, the transverse curvature 
effects are studied by Brobstein and Elliot (1956) and Tasuhara 
(1956)., Burther, Stewartson (1964), Dr ®rid Gross (1964) and 
Bllinwood and Mlrels (1969) have also studied the transverse cur- 
vature effects in hypersonic flow. 411 the above authors have 
formulated the problem by Including the transverse curvature 
terms in the first order equations* 

The effect of external vortlolty on the boundary layer was 
first pointed out by Berrl and Libby (1954). The first Investigators- 
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of the vorticity effect Li (1955» 1956) and Glauert (1957, 19^2) . 
engaged in a prolonged arguments, whether or not the pressure 
gradient on the' bounda,ry layer is altered. Hayes (1956) argued . 
that at hypersonic speeds the external flov/ has to pass through 
c. strong shock wa.ve, and, therefore may contain vorticity compa- 
rable to that in' the boundary layer. He studied the vorticity 
interaction in such a flow by modifying the outer boundary condi- 
tion in the classical boundary layer theorj'. With the assump- 
tion that the boundary laj^-er displacement does not appreciably 
change the pressure and velocity distribution along a streamline 
in the outer flow, he arrives at a modified outer boundary condi- 
tion for the velocity field f^ in the form 

f| = 1 + 2 H f , h -» 

v/here 44 is the vorticity interaction parameter, defined as, the 
ratio of the external vorticity to some average value of vorticity 
in the boundary layer. Kemp (1959) used the above condition to 
study the vorticity interaction. Here also, we wish to comment 
that such an analysis is empirical as mentioned in second para. 
Rott and Lenard (1959) have invariably disregarded the pressure 
change associated with vorticity, which is immertial for high ■■ 
cooled wall (with which these investigators were mainly concerned) 
he cause the ef fect is then slight . However, in general, this 
induced pressure gradient should be accounted for while studying 
the vorticity interaction. Herri, Zakkay and ling (196 1) studied 
the problem on the basis of local flat plate similarity. In par- 
ticular, they follow Lees (1956) in altogether neglecting the 
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tangential pressure gradient and external vorticity. Alfhougl 
the neglect pressure gradient has been shovm to yield reason- 
able accuracy for the first order boundary layer on a cooled 
wall, the significant role of press'ure in the second oi’der 
theory would suggest careful scrutiny the restilts, Further, 
their procedure of joining the inner and outer solutions for 
velocity and temperature profiles involves a. complica.ted patch- 
ing process along two different lines, rather than matching 
them asymptotically. Cheng (1961) has also' studied the vorti- 
city effect by using the thin shock layer approximation. 

Prandtl (1935) himself suggested, hov/ the boundary layer 
on a flat plate might be corrected by the effects of displace- 
ment. Massey and Cla.yton (1965) in their longitudinal curva- 
ture ana,lysis also include a consideration of displacement 
effect by formula-ting a nonlinear differential equation by a 
procedure which amounts to matching the outer limit of inner 

solution at the edge of boundary layer. This , clearly an over- 

; 

simplification. The displa.cement effects are more difficult 
to ca.lculate than others. The reason is that they are global 
in nature, being dependent upon the classical boundary layer 
solution in the whole range and its influence upon the outer 
inviscid solution, where as the other effects are local* To 
be exact the effect of displacement thickness should be solved' 
by perturbing the first order inviscid solution. However, in 
practice the flov; due to displacement thickness is studied by 
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solving the first order inviscid prohlem for a fictitious "body 
v;ith appropriate boundary conditions. Ihe profile of this fic- 
titious body is obtained by a suitable fit to the original body 
and its displacement thiclcnes= Thus Lenard (1961) has suggested 
approximating a sphere (locally near the stagnation point) and 
its displcioement thickness by a slightly larger concentric sphere. 
Malson ( 1961 ) proposes using a dilated and shifted, sphere. Davis 
and Dlugge-Lotz (1964b), DannelSp and Pliigge-Lotz (1966) and Honda 
and Kiyokowa. (1969) have also approximated their bodies and dis- 
placement thicknesses by such fictitious bodies. 

In global problems involving shock wa.ves, the flow in 
the shock layer, i.e., the region between the shock wave and 
the body is generally divided into an inviscid (outer) region 
and a viscous (inner) boundary layer near the body, and is 
solved according to the following procedure. ¥e first calcula- 
te the first order inviscid flow using the Rankine-Hugoniot 
relations as the outer boundary conditions, and then solve the 
classical boundary layer equations using the boundary conditions 
at the wall, matching it with first order inviscid solution, 

Hext, for the second order inviscid flow (due to displacement) 
we still use the Rankine-Hugoniot relations as the outer boxand- 
ary conditions. This displa.cement flow, in turn determines the 
second order boundary layer. Thus the second order boundary . • 
la 3 ''er ‘depends indirectly upon the Rankine-Hugoniot relations. 

The next higher order, i.e., the third order approximation to 
the outer flow requires that for the outer boundary condition. 
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the Rankine-Hugoniot relations he corrected due to curvature 
of shock wave (these corrections involve the study of the 
shock wave structure). The modified shock conditions and the 
associated analysis of the shock t ransition zone have been given 
hy Ting and Chow (1961) and Germain and Guiraud (1962), Oherai 
(1963) has calculated the corrections in the Rankine~Hugoniot 
relations in the stagnation region for the case of sphere in 
hypersonic flow. Further with the constant density assumption 
he has studied the third order effects in outer flow. Finally 
this outer flow will determine the third order boundary layer. 
Kao ( 1964 ) has calculated the effect of these coirections on the 
third order boimdary layer. 

Due to recent interest in high speed flights at high 
altitudes (for example, a space vehicle entering in the earth 
atmosphere encounters a wide range of flow conditions varying 
fromftee molecular flow at high altitudes to continuum at low 
altitudes) some authors attempted to fill the gap between the 
continuTim and free molecular flow regimes. The solution to des- 
cribe the entire (transitional) regime between continuum and 
free molecular flow, though very desirable, was foreseen to be 
a formidable task. So far the efforts have been made to obtain 
solutions for the following sub-regimes. (I) Fear free mole- 
cular flow: here, the procedure adopted has been to include 
first order collisions, i.e., collision with a free stream 
molecule by one that has already hit the surface. (2) Jlear con- 
tinuum; here, the (empirical) approach is that the governing 
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equations are those of the continuum regime, hut the houndary 
conditions near the wall are velocity slip and temperature 
jump as obtained from kinetic theory of gases. It is well- 
known that the velocity slip and temperature jump are propor- 
tional to the mean free path or in some sense to If the 

results of the two above analyses, put together could fill the 
gap the problem is completely solved. However, the range of 
validity of the solutions in these subregimes has not been de- 
finitely established. 

It is to be recognised that the final justification for 
the continuum analysis has to come from the study of the Boltz- 
mann equation. In clas=‘ical Chapman-Enskog expansion of the 
solution of Boltzmann equation (ref, Chapman and Cowling 1939) 
the leading term gives the inviscid Euler equations, the first 
two together provide the Navier-Stokes equations and the 
first three the Burnett equations. Though it would be natural 
to expect that the Burnett equations will provide better results 
than those by Navier-Stokes, It is found that in many situa- 
tions the results obtained by the use of Navier-Stokes equations 
compare better with the experiments than those by the Burnett 
equations; in some problems the use of Burnett equations gives 
anomalous solutions. Stewartson (1964) has observed that, 
provided the temperatures are not extreme and the Kundsen num- 
ber is small, there is no known case of disagreement between 
experiment and theory which can be traced to a fundamental 
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deficiency in the Na.vier-S.tokes equations. Even the flow pro- 
perties inside the moderately strong shocks are given fairly 
accurately by the llavier-Stokes equations, even though their 
thickness is only a few mean free path' and the corresponding 
Kundscn number -not very large. Of significance in this res- 
pect is the recent study of shock transition zone by liepmann, 
Earasimha and Chahine (1962) in which Boltzmann equation with 
Krook's collision model is used. The solutions obtained show 
small differences from Navier-Stokes solution on the high pres- 
sure side of the shock even for a Mach ninnber of ten. The con- 
vergence of the Chapman -In skog expansion has not been established. 
Moreover, it appears that the full depth of the properties of 
Eavier-Stokes equations have not been plumbed and that their 
usefulness extends beyond the theoretica,! range. Thus it is 
preferable to regard the Havier-Stokes equations as axiomatic 
with the constitutive properties of the fluid to be determined 
from experiments. Then develop the theory from the Havier-Stokes 
equations in as rigorous way as possible. The solution of Boltz- 
mann eqxiation (even for simple situation ) is a formidable task; 
so the further comparison with experiments and physical observa- 
tions can help us to decide the validity of above axiom. 

The main aim of the present work is to develop a general 
theory of the Gortler (1957) type for higher order boiuodary layer 
flows past ^ plane or axisymmetric bodies. In the spirit of what 
has been said earlier we shall be studying only the second order 
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effects .due to longitudinal curvature, transverse curvature 
externa'',! vorticity gra^dients, stagnation enthalpy gradient and 
displacement. Attention will be focu.sed on low speed and small 
temperature gradients, so that the velocity field is es'^entially 
that of an incompressible fluid. ¥e calculate the temperature 
profile too, (dissipation is ignored). We solve for first two 
terms in the second order G-ortler power series for velocity and 
temperotiire profile, by choosing the proper universal functions. 
The resulting system of total differential equations and the 
accompanying boundary conditions constitute the two point bound- 
ary value problems, and are solved by Runge-Kutta-G-ill method 
on IBM -7044 computer at Indian Institute of Technology, Kanpur, 
Solutions are obtained for both accelerating and decelerating 
flows with Balkner-Skan pressure gradient parameter between 
-0.198838 and 2 and the Prandtl number 0.7 to 3. Various other 
boundary layer characteristics such as displacement and momen- 
tum thicknesses etc. are also studied. 

However, at the beginning we have studied two problems 
of the first order boimdary layer, which have either not received 
sufficient attention in the literature or the treatment is not 
exhaustive or coherent. Eirst is the study of first order 
compressible boundary layer with no pressure gradient (using 
the power law relations for the dependence of viscosity, density 
and Prandtl number upon temperature) separately for the small, 
large and of order unity values of the Prandtl number. The 
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probl-em of low Prandtl number is of interest in liquid metals 
and plasma. For air the Prandtl number is 0.73, hut the flights 
at high speeds and altitudes lead to high temperature and low 
pressure at which the Prandtl number can became very small (see , 
inset" in Pig. 3.6 after Mitchell 196 2). The high' Prandtl 
number problem is of interest in liquids and vapours. It is 
shown in Chapters 3, 4 that the study of these limiting cases 
of small and l3.rge Prandtl numbers is of great use even under 
distinctly non-asymptotic conditions. ' ' 

The general problem of compressibl.e laminar boundary layer 
has been studied extensively in the past and a good review is 
given by Stev>/artson (1964). In general, it is foimd difficult 
to account for arbitrary fluid parameters and Mach number, and 
most of the authors assumed the Prandtl number of order unity. 
Except, in specific numericai calculations, accounts of which . 
are given by severs,! authors, notably Kuerti (1950), Young (1953), 
Van Driest ,( 1959) , Moore ( 1964 ) and Stewartson (1964); in most 
of the anaiyses it is assumed that the viscosity of the fluid 
varies linearly with the temperature. This is true from kinetic 
theory of gases, only, if the gas is made up of the so called. 
Maxwell molecules obeying the inverse fifth power law of 2 ?e pul- 
sion. Most of the gases, however, are not Maxwellian^. The 
viscosity - temperature relationship is better represented by 
simple power l&wP a T^ where w is a constant usua,lly between 
1/2 and 1. Though the Sutherland law can also be employed when ■ 
the interest is limited to a limited range of temparat-ure. 
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Here, we have used the method of matched as 3 nnptotic 
expansions for the low and high Prandtl number problems. In 
low Prandtl nimaber problem it turns out that the- outer flow 
is governed by a simple nonlinear equation of the second 
order, which is correct to all orders in Prandtl .number. In 
solving this outer equation one can proceed in the routine 
way by expanding the outer variable in powers of Prandtl num- 
ber. Because of certain simplifying features (see, 'Section 3.3) 
v/hich enables us to adopt a most efficient technique, i.e., we 
expfind the arbitrary constant of intigration instead of solu- 
tion. Por the case of insulated v/all this outer equation can 
he solved analytically, Howevei-, for the heat transfer case 
it is solved numerically hy Runge-Eutta method. Por high Prandtl 
number problem the inner flow is governed hy coupled nonlinear 
differential equations, with the outer boundary conditions to 
he determined hy matching it with the inner limit of outer so- 
lution (which is singular) and, in general, is very difficult 
to solve exactly. In this case an approximate solution of the 
Weyl t 3 rpe is carried out. Por Prandtl number of order unity 
we use a certain asymptotic method developed by Meksyu (I96l), 

In all above three cases of low, high and order unity Prandtl 
numbers, the leading terms are obtained for the recovery factor 
and the heat transfer rate at the wall. 

In the second problem, we have studied the heat transfer 
in the Palkner-Skan flows when dissipation is included in 



16 


the energy equation. The prohlem has been studied extensi-^-ely 
in the literature hut most of the authors he.ve either ignored 
the pressure gradient, the dissipation or both. The high 
Prandtl number flow is studied by Narasimha and Vasantha (1966) 
by ignoring the pressure gradient, and by Lighthill (1950), 

Adams (1962) and by Goddard and Acrivos (1966) by neglecting 
dissipation. In general, high Prandtl number case leads to 
a singular perturbation type problem; however, if dissipation 
is neglected, the problem becomes one of regular nerturbation 
type. The problem of low Prandtl number is studied is studied 
by Adams (1962), Goddard and Acrivos (1966) by ignoring the 
dissipation and Morgan, Pipkin and Warner- (1958) by neglect- 
ing that portion of the dissipation which is due to comipressibi- 
lity and by Afzal (196?) by neglecting the pressure gradient. 

Por Prandtl number of order unity, the problem is studied by 
Merk (1959) and Evans (1967). In the present work, we have 
used the method of matched asymptotic expansions for the low and 
high Prandtl number problems. For Prandtl number of unity the 
leksyn type method is employed. 

In the low and high Prandtl number flows, many higher 
order perturbations are calculated and the results obtained 
are good even up to Prandtl number of order unity, so the re- 
sults of both the low Snd high Prandtl number analyses overlap. 
Thus in a sense it suffices to analyse the two limiting cases: 
of low and high Prandtl numbers. 
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CHAHTER 2 

FO BMILATIOH OF SECOIE ORDER BOUNDARY LAYER TIEOIg 


2*1 Introduction : 

Y>Fe now describe a "brief resume of fhe derivation of the 
scheme for calcul8.ting the higher order corrections to the clas- 
sical "boundary layer theory^ designed for unreparated flows. "We 
shall follovv the treatment of Van Dyke (196 2). It consists of 
in developing the genera,! second order boundary la.yer theory, the 
governing equations for which are obtained from 15Ia.vier~Stokes 
equations by an expansion procedure in inverse power of the 
square root of Reynolds number. The procedure used, is the 
technique of matched asymptotic expansions developed by Kaplun , 
(1954, 1967), Lagerstrom and Cole (1955), etc. In this method 
we generally study the inner and outer limits of the equations 
o.nd try to match their solutions in. an overlap domain (see, for 
example Yan Dyke 1964, Kaplun 1967, Cole 1968)., ‘The inner and 
outer expansions (obtained by repeated a,pplication of these 
liro.it processes) result in replacing the Navier-Stokes equations 
by tv/o sGiparate sets of equations, one set which is valid in an 
outer region a.nd other valid in the inner viscous region (near 
the wa.ll in the problems discussed here). It is well known that 
none of the above two sets of equations are determinate in the 
sense that the number of boundary conditions are not enough. 
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The additional To oimdar;'.’- conditions required to maJce the system ' 
closed are obtained by what Kaplun (1957) calls the extension 
theorem. , ■ 

In the spirit of what has been said in the first chapter 
we shall be considering only the low speed flows with small tem- 
perature differences (so that the heat added by conduotion and 
dissipation are small thereby the density is sensibly constant, 
as are the transport properties). The governing equations for 
such a flow are IJavier-Stokes and the corresponding energy equa- 
tion. Furthermore, we shall be studying only plane or axisymme- 
tric flow past a semi-infinite body described by an analytic 
curve . 

2*2 Coordinate System ; 

We shall em^ploy an orthogonal coordinate system (s, n) 
shown in Fig. 2.1, where s is the distance along the surface 
(measured either from the stagnation point or the leading edge) 
and n is the distance normal to surface. The cori«sponding 
velocity components of the velocity vector U are (u, v), . Let 
E(s) be the curvature of the body section reckoned positive 
for a convex body. In axisymmetric flow, 0(s) is the angle ' 
between the axis and the tangent to the meridian curve at any ' 
point and r(s) is the distance of the point from axis. It is 
to be noted that all these quantities are. not -independent, but 
are connected according to 

sin e dr/ds, cos 0 = - K”^ d^r/ds^. 


( 2 . 1 ) 
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Tine leiigfh element dl in space is now given by 

dl^ = (1 + Kn)^ ds^ + dn^ + (r + n cos 0)^^ d0'; (2,2) 

Here j = 0 for plane flow and 1 = 1 for axisymmetric flow, and 
jZf is tbe cylindrical coordinate (linear distance normal to tbe 
plane of flow in two dimensional flow and the azimuthal angle 
in tlie axi symmetric flow). 


2.3 Q-ovemina equations: 

I'be basic equations of motion for a, steady plane or axisy- 
mmetric flow in the present coordinates are given in Van Dyke 
(1961). These are quite lengthy and will not be, given here, 
Alternatir^ely, y/e shall be using vector notation. The basic 
eduatlons of motion, for such a fiow of constant propeiiy fluid 
can be yo^itten, in nondimensional vector form (here all distances 
X are nondimensionalised with respect to characteristic length 
L of the body, velocity U and enthalpy h by the characteristic 
reference values and h^, and pressures P by p ) as 

div Tj = 0, ^2.3) 


U* grad IT + grad P - -R curl curl U, 


(2.4) 


U-grad (h + ETf/a) -E-'' ’ 9 O^h + BU^/S) 

= IR”** [ (grad U)*(grad (2.5) 

where we have written R for the Reynolds number (=pU^I,/g ), cr 
for the Prandtl number ( - CpA) and E for Eckert number 
(= R ) . Here is the coefficient of viscosity, p the den- 

sit y,:k the thermal conductivity and 0^ the specific heat at 


constant pre ssure . 
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Some times it is preferable to reduce the momentum eq^ua- 
tion (2.4) to two scaldr equations which are the generalisations 
of Bernoulli's and Kelvin's theorems, First forming the dot 
product of (2.4) with U, we get, 

U'grad (P + uVo) = R"'' U* (2.4a) 

where left hand side is the rate of change of total pressure 
P + U /2 alonv a. stream line. Taking the curl of (2,4), we get, 

(r + n cos 0)^ U' gradf *2 (r + n cos 0)“^ i (2,4b) 

‘ U- 

= R“^ div (r + n cos 0)“'"^ grad)o.(r+n cos 6)^ jj 

here again the left hand side is the rate of change of vorticity 
n along a stream line (it is recalled that for plane or axisymme- 
tric flow the vorticity;"^ has only one component which is always 
perpendicular to the section). 

The boundary condition at the solid surface for the 
velocity profile (v^ is the speed of suction or injection of the 
fluid and n is the unit vector in n direction) is 

U = n v^(s) (2.6a) 

which is the famous no slin condition for pejuie able surface, 
and for enthalpy profile is 

h = 9h/6n = 0, (2.6b) 

Only one of the two conditions at the surface in (2,6b) is 
required, the first prescribes the surface enthalpy to h^(s) 
and the second implies an insulated surface. For upstream the 
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flow is to approocb prescribed (possibly nommiform) velocity 
and enthalpy fields and respectively. 

2 * 4- for mulati on : 

The problem defined by equations (2.3), (2.4), (2.5) and 

(2.6) is highly complex and nonlinear. It is fore seen to be 

a formidable task to find the solution for entire range of 

Reynolds number. So far, the efforts have been made to obtain 

solutions for the limiting cases. We shall study the boundary 

-1 /2 

layer flov/ because this provides a natural small parameter R ' 
for studying the asymptotic solutions at large Reynolds number 
by the method of matched asymptotic expansions. 

2.4.t Outer Expansion : 

following Kaplun (1954), we define an outer limit (Euler 
limit) as the process R ^ with X fixed. A repeated application 
of the outer limit process in conjunction with appropriate se- 
quence of gauge functions gives the outer expansion. In general 
the outer expansion is straight forward. Eor unseparated flows 
past analytic bodies, it has been foiond that the sequence of 
gauge functions consist of inverse power of squaxa root of 
Reynolds number. The outer expansion therefore, has the form, 
TJ(X; R) = U^CX) + r-^/2 ^ ^ 

P(g; R) = P^(X) + R"^/^ + . . . 

h(X; R) = H^(3C) + R^v^^ H2(X) + . . . 

R) =r^.,(x)^^^+ * * * 

X (X5 = 'k ^(X) + R~^^^ ' * V. 


(2.7) 
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Here'i.|/ is the stream fxinction, defined in usual v/ay from equa- 
tion of continuity (2.3). We no?/ substitute the outer expansion 
(2.7) in equations (2.3), (2.4) and (2.5), and collecting the 
coefficients of the various powers of R give the equations for 
succespive outer approximations. Thus the first order equations, 


div U. =0 

(2.8a) 

U^’grad + grad P^ = 0 

(2.8b) 

U^'grad (H^ + = 0 

(2.8c) 


are the inviscid Euler equations of m.otion, and the second order 
equations are 

div U = 0 (2,9a) 

U^'grad JJ 2 + grad Pg = 0 (2.9b) 

U^'grad (H 2 +EU^-U 2 .).+ Wg'grad (H^+EU^ V2)=0. (2.9c) 

These are the displacement flow equations and are also inviscid 
in the sense that there do not contain the viscous term (of 
order R"'"^ in outer approximation). 

It is well known that Euler (outer) equations fail at 
the surface, through loss of the highest order derivative in 
the lavier-Stokes equations, in that, these can not satisfy the 
no slip condition. The outer expansion is invalid with in a 
distance of the order from the surface. To overcome 

this difficulty we need a different approximation, what Eaplun 
calls the inner limiting process^ l.e., the Prandtl boundary 
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layer type vrtiich is valid in the region of order near the 

body. Also it is observed by lagerstrom (1964) and Kaplun (1967) 
that for a given configuration the Etiler solution is not unique. 
The question is than which Euler solution for a given geometry 
of flow is the Euler (outer) limit of the Kavier-Stokes equations, 
^^nless ?i/e can guess it directly, finding the relevant solution 
would presuppose a Icnov/ledge of the solutions of Kavier-Stokes 
equations for infinite R; but this is exactly the solution to 
be found by perturbation procedure. The potential flow solution 
of Euler equation (also satisfies Kavicr-Stokes equations) is 
only one such possible solution. If the on coming stream con-' 
tains vorticity, a solution of Euler equation does not in general 
satisfy the Kavier-Stokes equations. However, it does so in spe- 
cial cases, for example in -plane flow with constant vorticity and 
axisymmetric flow with vorticity? proportional to radius. In 
other wot-ds out of a family of solutions of Euler equation we 
choose one which satisfies the far upstream boundary conditions 
and at the same time matches with the corresponding inner solu- 
tion (term matching will be discussed later in detail). Further?-: , 
more, the serious difficulties arise in finding the Euler solu- 
tion for non-analytical shapes (see, Goldstein I960). For se- 
parated flows the appropriate Euler solution is not kno'vm and 
for unseparated flows past finite bodies it has a complicated 
mathematical description. This motivates the restriction of 
our study to un separated flow past a semi- in finite body described 
by an analytic curve. 
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In the present work we shall not study the solution of 
the outer equations, assume that the first and the second order 
outer equ8.tions (2.8) and (2.9) with inner boundary conditions 
determined by matching ha.ve been solved (numeric ally or analy- 
tically), which provides at least, say, the slip velocities at 
the wall U^(s, 0) and U 2 (s, 0) respectively. How the pressure 
distribution etc. at the Yr&ll are foimd from the study of kine- 
matic equations (2.4a) and (2.4b). Substituting (2.7) in equa;^'" 
tions (2.4a) and (2.4b) and collecting the like powers of R, 
the first order terms are 

U^-grad (P., + U^^/2) = 0 (2.10a) 

‘gradl^Xk. ^ (r + n cos 0)“^ _ = 0. (2,10b) 

Now each of the equations (2.10a), (2.10b) and (2.8c) represents 
the conservation of total pressure, vorticity and total enthalpy 
respectively, along the first order stream lines described by 
7-:^ (obtained from first order continuity equation 2.8c) Further, 
the equations (2.10a) and (2.10b) are connected by Bernonlli 
equation (as a consequence of Euler equations 2.8), 

j!rL^(r + n cos 0)~^* = - d(P^ + ^ 

!»Keae equations (2.10a), (2.10b) and (2.8c) are integrated 
t o give , 


P.J + l]^/2 = B.j( .T,|) 


( 2 . 1 la) 

yL.|(r + n cos 9)"^ = - 


(2.11b) 

H., + ETI^^/2 = 


(2.11c) 
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where B^( and ( Tf^) are the Bernoulli and total enthalpy 

functions to he evaluated from upstream conditions. The second 
order kinematic equations are, 

IJ^«grad (P 2 + ^ ~ ^ ■ (2.12a) 

U^*grad[xi2(^ + ^ ] 

+ U 2 " gradjfl, ^ ( r + n cos J (2.121) 


Similar comments apply to the equations (2.12a), (2-12h) and 
(2.9c), i.e. , the flow quantities (total pressure etc.) are 
actually conserved alongr perturbed stream lines rather than 
the first order streamlines. When the on coming stream is 


prescribed, i.e., independent of Reynolds number the first inte 


gral of (2.12a), (2.12b) and (2.9c) is given by, 

n 2 (r + n cos 0)"^ = - 2 (^1^ 

Hj+EUg-Sl = ^"2 


( 2 . 13 a) 

(2.13b)‘ 

(2.13c) 


Row the pressure, vorticity and enthalpy at the wall from 
the first order equations (2.11) are 


T.j(s> 0) ~ B.|(0) ~ (®’ 

ni(s, 0) r-3 = - B'(0) 

H^(s, 0) = S^(0) - EU^ks, 0)/2 


( 2 , 14 a) 

(2.14b) 

( 2 . 14 c) 


and from the second order equations (2.13) we get, 

-B^isy 0) = B., '(0)'f'2(s, 0) - U.|(s, 0) 
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1 - 1 . 2 ^^’ »(0) ->^ 2 ( 8 , 0) (2.15b) 

HgCs, 0) = S^‘( 0 )'lJ^ 2 ^®’ 0) -EU^(s, 0) U 2 (s, 0).(2,15c) 

Here U^(s, 0) and UgCs, 0) are known from the solution of outer 
equa.tions and '^’" 2 ( 8 , 0 ) will he determined later. 

further, we shall he needina: the normal derivatives of 
the outer solution at the surf&ce, Evalusting the equations 
( 2 . 8 a) and ( 2 . 8 h) at the surface gives t?ro of the required nor- 
mal de rivat ii- e s , 

0 ) = -r“^ d l^r^ U^(s, 0 ) | / ds ( 2 . 16 a) 

F^^(s, 0) = K U^2(s, 0). (2. I 6 h) 

Ihe other normal derivatives are found by first differentiating 
( 2 . 11 a) and ( 2 . 11 c) with respect to n and then e vaHuating 
at the surface (using the definition of ^ •|n(s> * r^U^(s, 0 ) 

from 2 . 8 a) we obtain, 

Uin (s, 0) = - KU^(s, 0) + r^ B^'(O) (2.16c). 

Hin (s, 0) = U^(s, 0) [r^ S^>(0) - r% B^'(O) 

+ EKU^ (s, 0 )]. _ ( 2 .l 6 d) 

2 . 4 . 2 . Inner expansion: 

It is now clear that a different inner limiting process 
is to he constructed, which is valid in the region of order 
near the surface, following Prandtl we introduce the 
inner variable, 

• N = R^/^ 
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and study an inner (Prandtl) limit as the process R with R 
fixed. Again the repeated application of Prandtl limit process 
along with appropriate sequence of gauge functions will produce 
an inner expansion. For present case the appropriate inner 
expansions are, 

u(s, N; R) = u,j(s, E) + u^^s, E) + ... 

„5(s, E; R) = R"''/^w,^(s, E) + R"'' v^Cs, E)+. .. 

P(s, E; R) = p,,(s, E) + R"^/^ Pg (s, E) + . . . (2.17) 

h(s, E; R) = h^(s, E) + hg (s, E) + ... 

•\f(s, E; R) = R"^/^^,,(s, E) + R~'* -^2(3, E) + . , 

Substituting these expansions in equations (2,3), (2,4) 
and (2.5) 3-3ad collecting the eoefficients of like powers of R; 


we get for the first order approximation 

(r^Ui)g + (r^u,j)^ = 0 (2.18a) 

^^(u^) + Pis - 0 (2.18b) 

P-jP = 0 (2.18c) 

L^ (h^ + Bu^ 2/2) - ( (r-'’’h-, + Eu^V2)pj = 0 (2.18d) 

where (operator) = d/5s + 9/6N, 


and the subscripts s and E indicates the partial differentiation. 
These are the well-known Prandtl boundary layer equations. Simi- 
larly for the second approximation we get. 
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^ cos 9/r)lTu^ 

-fr^(K + 3 cos e/r) 1^ (2.19a) 

L^CUg) + ''■ ^28 “ "^ 21 ^ “ ^ 1 N 

- - ^i"^-}) (j COS e/r)u^jj (2.19'b) 

P 2 k = (2.19c) 

L^(h 2 + Eu^jUg) + L 2 (h^ + Eu^V 2 ) - ( 

= KEu^(h^+Eu^V2)3 + K( a^''Ti^-Eu^^/2)jj (2.19d) 

+ (3 cos e/r) ( cr^h^ + Eu^^/ 2 )jj ^ 


The above equations are the second order boundary l 8 .yer equations. 


The corresponding inner boundary conditions obtained from 
( 2 , 8 a) and ( 2 . 8 b) for Prandtl boundary layer equations are 

u^(s, 0 ) = 0 , v^(s, 0 ) = v^(s) ( 2 . 20 a) 

h^(s, 0 ) = h^(s) or h.,jj = 0 ( 2 . 20 b) 

and for the second order boundory layer theory are 


U 2 (s, 0 ) = 0 ^ 2 ( 3 , 0) = 0 ( 2 . 2 la) 

h2(5, 0) — 0 or h 2 jj-(s, 0 ) = 0, (2.21b) 


2 . 4.3 Matching: 

i /p 

The inner expansion is valid in a region of 0(R ) 

near the body and in general violates the far upstream condi- 
tions. The outer expansion is valid out side this region of 


violates some of the boundary conditions at the 
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surface (through the loss of highest order derivative of Uavier- 
Stokes equation in Euler If'dt). Hence, neither expansion has 
sufficient houndary conditions. This means that these missing 
boundary conditions must he replaced hy something which makes 
the problem determinate. This is done possibly by what Kaplun 
(1957) calls the extension theorem, which states that, if we 
know^ 

lim u = 0 

for some C ■> “tben these are other limiting process ((■ , in the 
neighbouihood of but not equivalent to it, such that 

lim u = 0 

also. Starting with this extension theorem we can easily estab- 
lish various forms of the matching principles used in practice, 
■first one is generally called the asymptotic matching principle, 

which states, 

The m-term inner expansion of (the p-th term outer 

expansion) = the p-texm outer expansion of (2.22) 
(the m-th term inner expansion). 


The second one states. 


Outer limit of the 


(2.23) 


Inner limit of the 

outer expansion inner expansion. 

Van Dyke (1962) applied the first matching principle (2.22) 
with m = p and m = p-1 and obtained the appropriate matching 
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conditions. .However, we sba.ll prefer to use the second one, i.e, , 
(2.23). We illustrate the procedure by the following example. 

let c|(s, n) be a t:^.TDical variable, whose outer expansion 
is given by, 

^(s, n| R) = §^(s, n) + n) 

+ R-^ ^^(s, n) + (2.24) 

Its inner asymptotic expansion is given by Taylor series (assum- 
ing it exists) around n = 0 and using the inner variable this 
becomes, 

^ (s, n - 0) = f/s, 0) + R-l/® [fg(s, 0) + 

2 

+ R"^ #3(3, 0) + 

On the other hand the outer limit of inner expansion (2. 17) 
is 9 

rj)(s, H; R) = N) + ER“^^^ 02(s, H) + 

+ R"'' ^ 2 ^ 3 ( 3 , H) + ... , N - «>. (2-25) 

Using the matching principle (2.23) we get 
0.^(s, N) = ^^{a, 0) 

02(s, H) = H 0) 4- <| 2 (®» as N - «> 

03(s, H) =^flnn O) +N$ 2 (s. 0 )+<f 3 (s, 0 )^ ( 2 . 26 ) 

Apply these results (2.26), the first one gives the matching 
condition for the Prandtl boundary layer equations, 
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(s, N) = (s, 0) 

N) = ^'-|(s, 0) 

as N - (2.27) 

P-](s, = P-jCs, 0) 

h^(s, U) = H^(s, 0) 


and the middle one gives for the second order boundary layer 
U2(s, N) = Nl^(s, 0) + UgCs, 0) 

V2(s, N) = ^2*^®’ 

%is, N) = 0) +-^2^s, 0) 

asn-*- (2.28) 

P 2 (s, H) = (s, 0) + P 2 (s, 0) 

hgCs, N) = M^^Cs, 0) + (s, 0). 


Various normal derivatives required here are already evaluated 
in section (2.4.2). Moreover, we also need 

0 ). It can easily^hown from above that 


•\J^p(s, 0) = lim ( 


(2.29a) 


V 2 (s, 0) = lim (v^ - ). (2.291 

The quantities on right hand side can be evaluated from the 
classical boundary layer theory. Physically, the above equa- 
tions represents the fact that the effect of the classical 
boundary layer upon the outer flow is that of distribution of 
sources over the body whose strength is proportional to the 
slope of the displacement thickness. 
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2^5 G-Q-y e rning Equati ons fo r Firs'b and Second __Orde_r ^ou^ 


Laver P ro blems : 

\Ye nov/ write the final fonr*. of the governing eq.uations 
for the first and second order honndary layer problems. 

2.5.1 First Order Boundary layer Problem: 

Integrating the normal momentum equation (2,18c) with 
respect to H- and eva-luating the function of integration from 
the matching condition (2.27c) we get, 


p^(s, N) = P^(s, 0), (2,30) 

This shows that the pressure p.j is constant across the first order 
boundary layer. Differentiating (2.30) with respect to s and 
using (2.14a) we get, 

Plg(s, B) = 0) = - U^(s, 0)U.jg(s, 0). (2.31) 

Thus the final form of the first order boundary layer 
equations is 


(r^u^)g + = 0 

Bl(u.,) - =TJ^(s, 0) U^3(s, 0) 

L^(h.j + Eu.j^/2) - ( = 0. 


(2,32a) 

(2.32b) 

(2.32c) 


and the corresponding boundary and matching conditions are 


U .J ( s , 0) 

= 0 = V^(s, 0) 

(2.35a) 

h.,(s, 0) 

= h^(s) or h.,jj(s, 0) = 0 

(2.33b) 

U^(s, N) 

= UJs, 0) 

1 [ as IT -♦ oo 

(2.33c) 


h^(si N) + Eu.l^(s, lif)/2 S^(0) j 
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2.5^2. Second Order Boiindary Layer ProLlem: 

Integrating the normal momentum equation (2.19c) with 
respect to U a-nd using (2.15)> (2.l6h) and (2.28d) gives, 

p^Cs, IT) = KmJ^2(s, 0) + K / 0) - u^^(s,M) ] dl 

+ B’(0) ~ U2(s,0)^ ( 2 . 34 a) 

Differentiating (2.34a) with respect to s and using the defini- 
tion of op-g we get, 

PpgCs, N) = [KNU^^(s, 0) + ^ U^^(s, 0) -u^^(s, N) dUlg 

- r^ B'(0) V2(s, 0) - [U^(s, 0)U2(s,0) J g 

Substituting the above into (2.19) gives the second ordejs 
boundary layer equations, 

(r^U2)g + cos 0/r)Nu^ ] ^ 

- [r^(K + D cos 0/r)Bv^ (2.55a) 

L^(ti2) + - ^2M = 0)Uis(B,0)7 

+ U^JJ - U^v^ ] + [kTJ^^(s, 0)+K / CU^ -» U| ) dlT ] g 
+ (3 cos 0/r) u^^ + r3B’^(0)V2(s,0) + [u.,(s,0)TJ2(s,0)]g 

.. (2.35b) 

L^thg + Eu^Ug) + E2^^1 ^ E^-|V2)-(cr^^i2‘^®'^1^2W 

= KIu^(h^+Eu^?/2)g + K(0 ^^E^-Eu^V2)h 

4- (3 cos 0/r) ( + Eu^^/2)jj , 


(2.350) 
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The second order boundary conditions are, 

U 2 (s, 0) = 0 = (2.56a) 

h 2 (s, 0) = 0 or h 2 |j(s, 0) = 0 (2.36b) 

and the matching conditions as obtained from (2.28) using (2.29) 
and (2.16) are, 

U 2 (s, N) = 17 [ r^B' ^(0) - K U^(s, 0) ] + U 2 (s, 0) 

as 37 

H) + Eu^(s, N) U 2 (s, 17) =^^ 1 ( 3 , 37) S*^(0) 


(2.36c) 



CHAPTI'R 5 


Fli^T 0RJ3ER CO&JPRESSIBLE BOmTDA RY LAYER OH A F IAT PLATE 
WITH ARB ITRARY- FLUID P ROPE^HOIES 

5*1 Introduction: 

Boundary layer flows of a viscous compreesiBle fluids 
with no pressure gradients depends mainly on tbe Prandtl number 
and the viscosity-temperature relationship. Baring, certain 
specific numerical calculations (see Chapter 1) all the avail- 
able results assume either the Prandtl number of order unity or 
the viscosity to be directly proportional to temperature. The 
aim of the present study is to fill this gap, by studying the 
flow of a variable thermal properties fluid for three different 
ranges (i.e., low, high and of order unity) of Prandtl number. 

To account for the variation of viscosity v/ith enthalpy 
we use a power law relation of the type d< h^. Further, as 
shown by Peng and Pindroh (1961), the behaviour of a gas at 
high temjjerature departs from that of a perfect gas. The gas 
law, then may be written as 

P = 

v/here Z is the compressibility factor and R_ is the gas cons- 

C “■ o 

tant. To simplify the problem, we will assume that this gas 
law can be wi’itten as p — ^ ^g~ ^ is approximated by a 

simple power law relation, for a given pressure ck T , ’a’ 
being a constant to be determined by curve fitting. Another 
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advantage of the law 9 at T®* is . that the seme study can he 
applied to liquids (with careful interpretation). 

Thus in the present study we use the following power law 
relations for viscosity p- , density ^ and Prandtl number <3r 
with enthalpy h 

p a h^, ^ a h®", 0 " a h^ 

Here, the indices w, a and b are to be determined from the curve 
fitting to experimental data or by some theoretical considera- 
tions. liirther, we study the flows for three different ranges 
of the Prandtl number, i.e.,, low, high and of order unity. The 
problem of low Prandtl number is of interest in liquid metal 
and plasma. Por air the Prandtl number at ordinary temperatures 
is 0.73. Due to high speed flights at high altitudes the tem- 
perature increases and pressure falls thus the Prandtl number 
can become very small (see inset in Fig. 4.6 after Mitchell 1962) 
The problem of high Prandtl number is of interest in liquid and 
some vapours. Moreover as said in the first chapter, the study 
of limiting cases is quite useful even under non asymptotic con- 
ditions. The method used for these limiting cases is that of 
matched asymptotic expansions. For Prandtl number of order 
unity we use certain asymptotic method developed by Meksyn (1961) 

The present work is not aimed at getting the exact asymp- 
totic solutions (though, in one - low Prandtl number - case the 
exact results are given) hut at finding the analytical depen- 
dence of the boundary layer characteristics on various 



parameters, in order to obtain a quick overall viev\^ of the 
relative importance of these parameters as well as a guide 
when looking for correlations. 


3.2 G-ovemina 


First order boundary layer equations for a steady flow 
of a viscous compressible fluid past a .flat plate in usual nota- 
tions are (see, Stewartson 1964? Van Dyke 1961) 


( 9 u) g + ( 

Q(uUg + vu^) = 


n 


= 0 




n'n 


Q(uh^ + vh^) 


R' 


■1 r 


{ M- fT” h ") 
vr (T -n-'n 


2 1 

] 


(3.1a) 

(3.1b) 

(3.3c) 


and the corresponding boundary conditions are 


n = 0 


n -> oo 


u 


u 1 


0 = V, 


h = h^ or h^ = 0 (3.2a) 


h 


1 . 


(3.2b) 


Pie re tlie coordinates s along the plate and n normal to it a.re 
non-dimensionalised b^/ some reference length 1 (as there is no 
characteristic length in the problem, L will not appear expli- 
citly). The velocity components u, v (in s and n directions), 
enthalpy h, viscosity jm. and density are nondimensionalised 
by their free stream values h^, and respectively. 

The number R is the Reynolds number U^l 

Introducing the stream function defined by 
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uhe equation (3«1a) is automatically satisfied while the equa- 
tions (3. Ih) and (3.1c) reduce to total differential equation 
if we assume ''f and h to be of the form 





n) = 

(2s/R)^/^ f(yj ) 





h(s, 

n) = 

h ( rj ) 


(3. 

3) 

with 


1 

= (R/2s 

/ 

n 

? . 







0 



The 

reduced 

(tot a 

1) differential 

equations are 






!! , f 

f ) + 

f f 

ff 

0 

(3. 

4) 



( cr” 

'pp.h ) 

+ fh' 

+ C (^juLf = 0 . 

(3. 

5) 

Also 

the boundary 

conditions (3. 

2) become 




f(o) 

= 0 

= f'(o), 

f (°°) = 1, 

(3. 

6 a) 


h(0) 

“ ^w 

f 

or h 

(0) = 

0, h(^) = 1 . 

(3. 

6b) 


Here dashes denote the differentiation with respect to . 

To account for the variation of ^ and CTwith h, we 
use the power law relations 


II 

ll 

= h\ 

Eliminating |ju 

5 ^ and cr in favour o f h 

from (3.4) and (3 

get 



III tf 

f ■•+ f 

[^-(a+w) + (a + w) h'/b 

1 = 0 

T« » 

h- + h- I 

, i 

[fl-. - (b - a 

- w) h'/h ] 


b "2 

+ G 0^ f = 0 . 


(3.7) 
.5) we 


(3.8) 


(3.9) 
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Further introducing the transformo-tion. 


h = , m = 

the equations (3.8) and (3.9) reduce to 

f + f" (fe® - e e V©) = 0 


b - . 

i + a + w — b 


( 3 . 10 ) 


^ ^ TTl ^ d — TTl 

9 + .r- f0 0 + fT- 0 0 f 

U oo ^ u DO U 


(3.11) 

0 ( 5 . 12 ) 

where e = - (a + w) ( 1 + m) , d = b( 1 + m) and C = C/(l+m), 
The boundary conditions for the velocity profile are same as 
(3.Sa) and those for temperature profile are 

2(0) = e^, 

3 . 3 Solu t i on for Low Prandtl Number 


= h or 9* (0) = 0, Q(<^) = 1.(3. 13) 

W 


Edward and 'Tellep (1961) ha,ve analysed the heat transfer 
in low Prandtl number fluids with va.riable therma.l properties 
by ignoring the viscous dissiuation and pressure gradient. To , 
a,ccount for the varia-tion of thermal properties with tempera- 
ture these authors have used for the term energy 

equation a power low dependence of the type T^ (in the nota- 
tions of Edward and Tellep), They have solved the problem, for 
the range of m between zero and unity (0 _< m _< 1). However, for 
a perfect gas and most of the fluids m is usually between zero 
and minus one (-1 2 ^^2 O)- this range that we '.-tie re 

study the head transfer problem with full account of dissipa- 
tion. In terminology of the method of matched asymptotic ex- 
pansions (as employed here), Edward and Tellep have solved the 
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so called, the zeroth order outer ec^uation; their solution is 
however, uniformly valid up to the order for heat 

transfer case. To get the higher order terms it is essential 
to formulate a different inner limit of the e questions. In the 
present work we have formulated and solved the equations for 
the two regions with full account of dissipation. Also the 
study of an insulated wall .in the flow of a fluid with arbi- 


trary pi'-operties is included. 


We now solve the equations (3.11) and (3.12) under the 
boundary conditions (3.6a) and (3.13) a-t low CTlo* Following 
Lagerstrom (1964) we proceed as follows. 


3.3.1 Inner limit: The inner limit is defined as 0 with 

fixed, and we write 



n=0 ' 

(3.14a) 

9( 1 . (r«) 

= s e„(r]) erf 
n=0 ^ ! 

(3.14b) 


Substituting the above series (3.14) ir equation (^i'. 11 ) and 
(3.12) and collectinf?, the coefficient of the like power of 
the momentum equation (3-11) gives, 


in ti. . ^ 

f + f '■ ( f 0“^ 
o o 0 0 




0 


(3. 15a) 


e-1 


- ee’/Bg + ee'ege^) = o, (3.i5b) 


1 
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and the energy equation (3-12) givej 


0 


0^=0 


■f e“ - -5 e f"2 

0 0 O 0 


The solutions to a.'bo"''e three equations (3.16) are 


fl 4- b 

^ I 0 


a . ^ + b . 


^ii ' -1 

( ■"! - x)(a^x + h^)™ fQ(x) dx 


- C /'(■>! - x)(a^x + b^)^ “ f^^(x)dx + ®-2»1 ’’^2 


(3- 16a) 
(5. 16b) 
(3. 16c) 


(3. 17a) 

(3.17b) 

(3.17c) 


Here a's and b's are constant of integration either of them to 
be determined from inner boundary conditions.. 


All these solutions (3.17) are singular for large y| a.nd 
do not satisfy the boundary conditions at infinity. This singu- 
larity is rather similar to one encountered in improving Stokes 
solution for low Peynolds number flow. 


3-3.2 Outer limit; We now need a different outer limit. From 
an order of magnitude analysis we introduce the outer variables 

and study the limit -* 0 with ^ » S' and ^ fixed- With this 
outer limit the equations (3-12) and (3.13) became 
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S’ f + cry^ = 0 (3.19a) 

+ C =0. (3.19t) 

It can be shown easily that the solution o.f the momentum 
equation (3. 19a) which satisfies the boundary condition at 
infinity is 

F = a + 0(cr-’'"). (3.20) 

.*-.00 

Miere cr' denotes the exponentially small terms in the limit 
CJ" Nov/ with the help of (3.20), the outer energy equation 
(3.19b) may be v/ritten as 

+ (C - = 0(e-“). (3.21) 

The outer boundary condition is 0(°°) = 1. The nonlinear outer 
equation (3.21) is correct to all orders in i.e., the 

error is exponentially sma.ll. This is due to the fact that 
there is a small momentum, layer inside the thick thermal 
layer and all that the momentum layer does far av/ay is to 
displaice the stream lines from their inviscid position hy an 
amount a. The outer equations of all orders can be obtained 
from equation (3.21). When m = 0, the solution to this equa- 
tion (-5 .-21) which satisfies the boundary condition at infinity 
is 

0 = A + (1 - A) eTf[(C^- 

where erf (x) is the well-known error function defined by 

X 2 

erf (x) = f e~^ dt 

o • 

end A = E 
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It appears that for a non zero m, in general, equation 

(3.21) has to be solved numerically. However, it is interesting 
to see that for an insulated wall the solution of zeroth order 
outer equation = 1 matches with the corresponding inner 

solution; so that it is not necessary to solve the equation 

(3.21) numerically. Therefore, we now study the outer equation 

(3.21) separately for the cases of the insulated wall and heat 
transfer. 


3.3.4: Insulated wall case: Before solving the nonlinear outer 
equation (3.21), let vlb recall that for an insulated wall 
(O) = 0, and the inner equations (3.17) gives a^=a^=:.,=0. 

New the outer limit of the inner solution for large oj may be 
written an 

y OO 

0(^ co) = +cry^ (^> 1 + C / f^Cx) dx)+..^ (3.22) 

Fe now proceed to solve the outer equation (3.21). let 


us assume 


^ n=0 ^ 

Substituting this expansion (3.23) in equation (3.21) and 
collecting the coefficients of the like powers of cn„, we 
get 




(3.24h) 
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+ om0-|0Q^ /0Q (3.24c) 

-(m/2)(m-l)^02 0^^ /02). 

The "bo-undary conditions at infinity are 

0^ (-) = 1, 0^(°-) = 0, 0^(-) = 0, (3.25) 

The solution of (3.24a) v/hich satisfies the boundary 
condition 0 ^(°°) - 1 and matches with (3.22) upto the zeroth 
order is 

0y^) = 1 = . (5.26) 

Q'he first order outer equation (3.24‘b) no?\r hecoiDes 

= ° 

and its solution v/ith boundary condition 0.^(°°) = 0 is 

0^ = A erfc (^/Y'2) (3.27) 

where erfc (x) is the complementary error function defined by 

X ,2 

erfc (x) = 1-(2/)r-ji;) / e“ dt. 

0 

Tor small ^ the inner expansion of (3.27) is 

0^(^- 0) = A (l-(2/ 7t)'*/\ +...). (3.28) 

Matching (5.28) with the first order terms in (3.22); we get 

b^ = A = "C ('n;/2)00^ / f ” ^ (x) dx, (3.29) 

Now with the help of the relation (3.26) and the solution (3.9) 
of the outer momentum equation (5.15a) reduces to the BlasiTis 
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equation; we denote the solution "by B()|). 

The adiabatic wall temperature is 

b„ = 1 + C f b”2(x) dx, (3.30a) 

o 

and the recovery factor is 
r^ = .2(h^ - 1)/C 

= t(2%c^) f b "^^ (x) dx = 0.9255 (3.30b) 

o 

This last result is independent of thermal properties of the 
fluid. 

5.3.5, Heat Transfer Case: In this case the nonlinear outer 

equation (3.21) can be solved exactly. It is convenient to 

/ 

change the following i/independent variables 

2 = a. 

The equation (3.21) reduces to 

0^2 + 20’®0^ = 0 (3.31) 

and the corresponding boundary condition is 0(°°) =1. 

To solve the equation (3-31) (which is still correct 
to all order in CT^) numerically, it is first to be transformed 
into a marching problem, by setting 

0{z) = 1 T(t), 2 = q^t (3.32a) 

where 1 and q,^ ar« constants. Substituting (5.32a) in (3-31), 
v\re get. 
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t = 0. (5.32b) 

If 1 - 1 

the equation (3. 32b) reduces to 

+ t t“ = 0. (3.32c) 


This equation (3.32c) is similar to equation (3.31) but the 
boundary condition at infinit' is changed, which is 1 ( 0 °) = I/I. 
This makes possible to generate the family of solutions. Now 
for a given 'm' assume a value of T(0) = 1 and set the values 
for T^(0). Tor each value of T^(0) the equation (3.32c) is 
solved numerically by Runge-Kutta method (on IBM '7044 computer 
at Indian Institute of Technology Kanpur) to get T(°°). Knowing 
T(°°) the constants 1 and q^ are determined by 


1 = (T(oo))-'’ 

and q^ = (T(°o))®'/^. 


(3.32d) 


The equation (3.32c) and its solution may now be transformed to 
original variables 0 and z with the aid of (3.32b5 d). Solutions 
are obtained for a. range of values of m between zero and unity 
(0 < m < 1). The tabulation of all the results of the numerical 
computation is rather difficult. Instead we shall be displaying 
in Tig. 5.3 the initial values of each of the solutions. Only 
some of the solutions for outer temperature profile are shown in 
Tigs. 3.1 and 3.2. 
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The results of numerical solution, neces'^arj?' for matching 
are plotted in Fig. 3.3, which displays 7^0(0) for 

various values of m. To Xffect the matching, we write for- 
mally the inner expansion of the , (numerical) outer solution 
as ■ ’ 

0) = 0(0) +(ry^ (■»]- «) 0^(0) + (3.33) 

Here 0(0) and 0 (0) are constants independent of z hut are 
functions of cTlo? 

0(0)= £ ^ c er”/®. (3.54) 

n=0 ^ ^ n=0 

Substituting (3.34) ,.:in (3.33), we get the inner expansion of 
the outer solution as 0 

0(^- 0) = +<ry2 (A^ - aO^ +l|0^) 

+ Ol. (Ag -<>«., + 'pep + 0(0-^^) . (3.35) 

First matching the zeroth order inner solution (3.17a.) with 
the zeroth order outer solution (3.35) we get a^ = 0- Next 
the outer expansion of the inner solution (3.17) may he written 
as (hy the inner boundary condition h^ = 6^, h;^ ~ ”*^2 ~ 

e( -rj - =o) = 6^ + + o-yt^a^+Cey” / f^Cx) dx) 

- 0 / X ^ (x) dx ) + 0 ( (3.36) 

O 

Matching (3.35) and (3,36) we get, . 

% " -^O “ 

”^1 = °o ■ = °o“ 
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a„ = C ~ C / f ^(x) dx 

Ag = C^a-’C 0^““ / x,f”2(x)dx. 

It can be sbov/n ea,sily that the matching result (3.37) together 
with the solution of the outer equation (3.20) and the transfer 


mat ion 

foCvj) = B(^) ey/2 , ^ 9-n/2 (5_38) 

reduce the zeroth order inner momentum equation (3.15a) to the 
well-fcnovff:i Blasins equa,tion with the solution B('>l). 


The heat transfer rate at the wall is 
t, i/o r _ d-m+(n/2) 

0 (O) ^ cT" ^ C I + G0 

^ ^ oo 0 i 1-0 Q 


+ • * • 


oo -n 

/ B*^(x)dx 1 

Pi -I 

(3.39) 


Introducing the reco-^^ery temperature (3.30) we get 


b (0) 


1/2 

<T' 

OO 


m 

'pr 


+ f(2/n) (h -1) h^ 


m , b n 

2(a+i) 2 


+ 0(crj (3.40) 


let us rewrite the above relation as 


h'(0) = +-/I (h^-1) 


m , b 

^JTTSI 2 cr-J/2 (3.41) 


m 

= (m+1) 0 JO ) o - y ^. 


2 


(3.42) 


where 
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The quantity may be recognized as the heat transfer rate if 
viscous dissipation is ignored. This factor W^/ (T^ ' is 
plotted against h^ in fig. 3.6 with m as a parameter.* for a 
given value of h^, we first determine 0(0) = Next 

from fig. 3.3 we obtain for a given 0(0) and m. This 

procedure determines W^. Knowing the recovery temperature, 
heat transfer can be determined from (3.43). 

for the case of perfect gas (a = -1) and viscosity pro- 
portional to temperature (w =1) and constant Prandtl num- 
her (b = 0) this becomes, 

\ = 0^(0) = (1 - h^) r(2cr/7t) 

and 

l^'(O) = (h^ - h^) f(2cT-/7i) . 

This last result is same as given in Stewartson (1964). 

The coefficient of skin friction is defined as 

1 /2 

With the help of (3.38), (3.7) and (3.10) the above result 
becomes 

Cj. = fz h^(''+®V2 b"(o) + 0 ( 

(5.4+) 

The uniformly valid solution to enthalpy profile is 

) = \ + C o;; ^ 

.+(1+m) .(Y|l a)^ -0 


45 ) 
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The above soltition shows that the temperature increases like 

as 0. Thus it not only fails to satisfy the boiuidary 
condition at the wall, but is singular there. 

3.4.2 Inner limit: 


We need now a different inner limit process. Irom an 
order of magnitude analysis we introduce the following inner 
variables 



= 1 /e 


X 

_ Q ^ -l/(5m-4e+3-3d) 

- 3. ^ . 

(3.49) 

p 

^ (4m-3e-2d+1)/(5m-4e-3d+3) 


e 

_ -(2m-e-d+l)/(5m-4e-3d+3) 



and study the lir;it CT^ -'OO with , P and 7^ fixed. Assuming 


P = ■ Pq + G P.| + . , . 

X = + e X| + . . . 


the zeroth order inner equations are 


P - e P X / y 


0 


m ^ f , M d'^m 


r + P Y V + C -V 

■O^^ 0^0 Ao^ A 


P 


0 




= 0 , 


(3.50) 


(3..5la) 

( 3 . 511 ) 


These equations (3.51a) and (3.511) are coupled and nonlinear 
with the outer boundary conditions to be determined by matching. 

It turns out that the outer boundary condition for TCq 
singular at large . Thus the equa.tions (3.51) snd the accomp- 
anying boundary conditions lead^ to, very difficult botmda ry . 

value problem. As our purpose is to find the anaijfti.^^ 1 .• 
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dependence, it is instructive to carry ont just the first 
step in Yfeyl’s iteration scheme. To accomplish this, we 
first (formally) integrate the equations (5.52a) and 


(5.521)) to yield 


^0 = (y- 


X 


2m-d+1 


o 


2m-d/ 


2m-d+1 


.-P(x) 


= A. + A.J f \ (x) dx 


o 


0 


(3.52a) 


(3.52-b) 


- C (t) 

0 0 it 




0 

dt dx 


where F(^) = d^ ^ 

Upto this stcige the analysis is exact. To stud 3 )' the first 
step in Weyl’s iteration scheme, we use only one term in 
the expansion of Xq wall, i.e. It is observed 

in view of the fact, that the Prandtl number is high, the 
thermal boundary layer is very thin and is contadned in the 
viscous boundary layer which itself ha,s small thickness, the 
above approximstion should be 'very good in the inner region 
of thermal boundary layer. 


for this approximation 
a V® 

"sXw ,,5 


F 

o 

_ . 2m-d+1 

Xo 


;» 


f = fs ^m+e 

3 1 Xw “ 


y 2m-d+1 

Aw 


2 A 


2m-d+1 


(36a^) 


1 X 

T7f w 


w 

5m-e-5d 


i( 


2G 

1+S 


F) 


(a^/6)'’/^ N.,/3( l/5y f) (3, 53) 



where 


, 53 


x) = / e"^ dt 

0 

is the WGll-known incomplete gamma function and 

X 

IT (a x) = / t^ lT(a, t) dt 

ill 5 0 ' 

and uCa, x) = af e“^'*'('1-t)°^ dt 

0 ■ 

a.re certain related gamma factions studied in Appendix. A, 


Further, the outer expansion of this inner solution is 


^2ro— d+1 
Xo 


,2m-d+1 

"w 


2ca 


s 


A 


5m-e-3d 


2m-d+1 


(36ag) 


1 X 


_ 12c , % U/3 X 5 
m+1 ^ ^ ' 






( 1+m) 

4m4-e 

t 

Now the inner expansion of the outer solution (3.48) in inner 
variables (3.50) may be written- as 


X 


2oa V 


■W 


2m-d+1 


(l+m)' 


1 

_ 2m-d+l 


• cJ(o) 

, m+1 


^ 5m-4e-3ii+3 

oo 

(3,55) 


where , ■ 

°° II o 

is regular near the wall and a„. = . 

S 1 e w • 

3.4.3 Results: 

Various vinknown constants will now be determined by 
matching the inner and outer solutions in the overlap domain. 
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I'or the case of an ins'ulated wall' X (O) = 0> requlires 
= 0 and the matching condition (2.24) gives the wall 
recovery factor 

■llatl)- 


""f = C 


^2. 

V 


(aV6)1/5 


ir 


.V5 


5m-4e-3d+3 


1+m °°) cr;® 

. -i/t; i /■7:T3/(2w+2a-h+3) 

2C(l-a-w) (a^/6)^/^ N_^y^(l/3-, “) 


= ••• (3.56) 

X o 

When wall temperature is prescribed to a constant value, 
the matching principle (2.24) determines 

1/5 

l ] c-lii— u. ' I 


2m-d+1 2m-d+1 

A. = ( X - X ) 

1 V /^r ^ ' 


(2m-d+1) fri/3) 
2m-e-d+2 


(e-5m+3d)/3 


- - f|©Vfr^ 5m-4e-3d+3 

" =0 " °o 


The heat transfer rate at the wall is 

m-d 

m+l 1 (36a ) 

h^(hyh^) -h, I = 


h (0) = 


W 


1/3 


e+m 


V ■^ ••• 




-a-w 


■h. 


w 


a |ri/3) 

(36a„)''/^ 4(l3-2a-2w) 


h: 


Ti-sXwXIT^TXT "w 

To study the above result (s-S'Ta) in detail let us write 

% 1 


CToo •• 


h (0) = Wj^ 


V5 j JjeapX^MCb-Sa-Sw) 

w 


w,(...) = (.— - V A 


(1-a-w) |T^1/3) 


1/3 ■^(b-2a-2w) ^y 3 


.-•(3.57) 

crX^ (3.58) 


where 
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is the contributiojn to heat transfer when dissipation is neg- - 
lected. Knowing W^(m, h^) and the recovery temperature actual 
heat transfer may he obtained from (3-57). 

low for a perfect gas (a = -1) vrith viscosity propor- 
tional to temperature (w = 1) and constant frandtl number 
(b = O), the result (5.56) for the recovery factor reduces to 

r^ = 1.922cr/^ (3.59) 

and (3.56) for the heat transfer becomes 

h’(0) = 1.6227 (h^ - h^ 0-^/5 (3-59b) 

There results (3- 59a) and (■3.59b) are same as given in 
Stewart son (1964). 

3.5 Solution for Prandtl l\xmber of Order Unity ; 

In the previous sections (5.3) and (5.4) it is seen that 
the low and high Prandtl niunber problems are of singular pertur- 
bation type. However, when Prandtl niimber is of order unity, 
the thickness of momentum boimdary layer is of the same order as 
that of thermal bounda.ry layer, so, in this case the problem 
is expected to be a regular perturbation type, low to solve 

f’"+ f" (fh® - e h /h) =0 (5.60a) 

and (o^ ffi'* - qh'/h) + 0 = 0 (3.60b) 

where e = -(a+w) and q = b-a-w, with the boundary conditions 
(3.6a, b) , when Prandtl number is of order unity we utilize 
a certain asymptotic method developed and used extensively, 
on many problems , by Meksyn (1961) . 
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The equations (3.60a) and (3. 6 Oh) may formally be 
integrated to give 

f ' = a /^(h/h^)® exp (-F^) d y| , (3.6la) 

b = Aq + A^ exp(^5Fj^) dr| -ca^a ^(yj) (5*62b) 

where 

*'1 ^ ^ exp(4ffP ) d 

r 0 wv 

exp 

and ¥(y\) = fh^ dn 

, r-‘ ' 0 t ' 

In writing down the equation (3.6 la) we have used the boundary 
condition at the wa,ll and the condition at infinity is to be 
satisfied by a proper choice of a. The constant of integration 
A^ and A^ in (3.62b) are yet to be determined, 

I 

For an insulated wall, the condition h (0) =0 leads to 

A^ = 0 and condition at infinity gives 

_2 

h = 1 + C (j;; a (oo) _ (yj) j 

The wall recovery enthalpy is 

-2 . 

\ = 1 + C(3r^a ^(-)\ . (3.63) 

When wall temperature is prescribed h(0)=h^=lQ and the 
condition at infinity determines A^. Hence the solution to 
enthalpy profile is 

-0(T;,a4(Tj) 





2e-q+d 


0 


h. 


w 
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Heat transfer rate at the wall is 

oo 

h’( 0 ) = (h^ - h^^) / / exp (-crS'^) . ( 3 . 64 -) 

As our aim is only to obtain an expression for the 
analytical dependence of boundary layer characteristics on 
various perameter it is enough to use only one term in the 
expansion of f and h near the wall, viz . 

f = , h=h^ + .,.^ (3.65) 

Here h = h for an insulated wall and h„ = h , when the wall 
0 r 0 w 

enthalpy is prescribed. 


How to obtain the value of ’a' we make use ( 3 . 6 la), 
(3.65) and boundary condition at infinity gives 

OO 

1 = a / 

0 

r 9 h 


u X ^ 

[2 

The coefficient of skin friction is given by 
Of If R3 = T2 a 

= . 0.684 

where is local Reynolds number. 

The integral^ is evaluated by using the following 

relation proved in Appendix A 

2 

/ 2® exp (-PZ) / t^ exp (Pt - at) dt dz (3.67) 

■ . 0 : 0 

. m+n+2, n+2, 


f 3.-0 

exp ( - g- 


e 

o 


1/2 


) drj 


1 / 3 ) - 


= a h®'^^ 
so* 


( 3 . 66 ) 



58 


where (a, h, c; x) is the hyper geometric function. Now 

the expression (3.63) for the recovery enthalpy becomes 


^ ^ 1 I 3 [r2/3) 2^1 ( '1/3, 2/3, 4/3 

^ P^(l/3) 




The recovery factor ±£ 


r^= 1.29 'gJ’-t- (V3, 2/3., 4/3, 


?/hen b = 0 and 0^ 


1, this gives r = 1.11 


For heat transfer case the equation (9.64) becomes 

1 /^3 

l^l>-2a-2w)/3 g-y? _ ( 3 _TO) 

Here again to separate the effect of Mach number from that of 
wall temperature we rewrite (3.70) a 

■ (36a 


h (0) = W 


(b,a.,w; h^) aV5 - (l-hj 


2n’l/3) 


(b-2a-2w)/3 ^1/3 

^ oo 


(5.71) 


where 


= 


2fri/3) 




(lD-2a-2w)/3 <ry^ 


(3.72) 


W' W 


is the contribution to heat transfer from wall temperature. 


Lastly, apart from the approximate analysis described 
in this section (3.5), the energy equation (3.60b) can be 
integrated exactly when <^= 1 and b =0. To accomplish this, 
^ we eliminate the factor fh® - eh/h from equations (3.60a) and 
(3.60b) to' obtain 
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!! ! tn ff 9 

h - h f + C f ^ = 0 . (3. 73) 

The selution to this equation is 

h = + Ig f - I (3. 74) 

v^here and Ig are constants of integration The "boundary 
condition h (°°) = 1 gives 

+ I 2 ~ 1' + C/2 „ 

For an insulated 'wall Ig = and the entbalpv profile is 

h = 1 +1 (1 - f'^). (3.75^ 

The wall enthalpy is h^ = 1 + C/2 and the recovery factor is 
r = 1 . ( 3 , 76 ) 

For heat transfer case 1. = h„ and 1^ = h_ - h„ and 
the enthalpy profile is 

h = tw + (hj, - hp f' - Of'^/2. (3.77’) 

The heat transfer at the wall is 

fe’(O) = (h^ - h^) f‘’(o) (3,78) 

and the Reynolds analogy factor is 

m = — — = 1 . ( 3 . 79 ) 

(h^-h^f ( 0 ) 

These results (5.76) and (3.79) are same as clas!=^ical Crocco’s 
integral. Thus at CT" = -j ^ the recovery factor and the Reynolds 
analogy factors are independent of thenrsl properties. The 
results ( 3 . 68 ) and ( 3 . 72 ) for the perfect gas (a - - 1 ) and cons- 
tant Prandtl numher (h = 0) reduces to those of Yasantha and 

: ■ 

R aras imha ( 1 96 9 ) » 



The main resnlts for low, high and order unity ranges 
of Prandtl numher are contained respectively in (3.30)5 (3.56) 
and (3.69) for the recovery factor and in (3.41), (3.57) and 
(3.70) for the heat transfer. The main interest centers on 
the leading terms, as these indicate the extent to which the 
results are useful. Por the case, of a perfect gas (a = -1) 
with viscosity proportional to temperature (w = 1 ) and constant 
Prandtl numher (h = 0 ) our results reduce to those given by 
various authors,, It is interesting to note that the result 
(3.30) for the recovery factor at low (j^ upto the order 

is independent of the thermal properties of the fluid. 

It seems that the precise dependence of ^ and (f on the 
enthalpy h is not crucial when the Prandtl number is low. This 
is plausible, because the plate is thermally insulated and the 
fluid is highly conducting (low CT^) ; therefore the variation 
of temperature across the boundary is small. Further, the 
result (3.69) for the recovery factor when the Prandtl number 
is of order unity, is independent of w and a and depends on b 
only. Thus, here it also appears that precise viscosity and 
density laws a,re not crucial. The numerical solution of 
Crocco (1941) for Prandtl number of order unity indicate that 
the corresponding numerical values using different viscosity 
lav/s are close together. Further, Kuerti’s ( 1950 ) review 
makes it clear that when wall is insulated and the Prandtl 
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number held fixed near unity, the temperature variation in 
the boundary layer for various viscosity laws is approximately 
the same function of velocity. However, the result (3.56) for 
the recovery factor when Prandtl ntunber is high depends strongly 
on w, a, b (the form of property relations of the fluid) as 
well as on Mach number. The effect of the thermal properties 
on r^ (proportional to recovery factor) is shown in Pig. 5.5. 

It is seen that "r^ increases as m increases. 


Prom Pig. 3.4 which displays the recovery factor for 
above three cases it is seen tha.t low Prandtl number analysis 
gives 0.9255 (the exact result of Pohlhansen (1921) is unity) 
while at cr= 15 gives r = 3.59 (the exact iresult of Pohlhansen 
is 3 . 54 ). It is fortuitous that quite useful estimates of 
recovery factor can be obtained even for CQ= 25 with in 8 per 
cent of error. The expression (3.69) for the Prandtl number 
of order unity is close to the result of Pohlhansen at 
cr= 0.7 and 15. The result of Spence (i960) for a = -i, w = 1 
and b = 0; viz, 

r = 1^1 + 0.0095 (cr- 1)^ + '•*** 1 

is also shown in the Pig. 3.4. .It is seen that the Spence's 
result show a marked divergence from those of Pohlhansen for 
appreciably different from unity; thus, the region of validity 
of Spencels result is rather limited. 

^ ^ ^ ire suits ( 5 . 41 ), ( 3 . 57 ) 2 nd (3.70) for the heat 

transfer for the three cases of low, high and order unity 
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Prandtl nxunbers - depends upon w, a and "b, and also on aiach 
number through the recovery enthalpy: the form of these results 
are quite complicated. However, we can get some insight hy dis- 
cussing the factors and Wq defined respectively by (3-42), 

(3,58) and (3.72). These quantities are proportional to heat 
transfer rate if dissipation is neglected- The contribution 
of dissipation may be added to respective W’ s to get the total 
heat transfer. 

For low flows the quantity W^/ vs. h^^ is dis- 
played in Fig. 3.6 for various values of m. It is seen that as 
m. increases the quantity decreases for h^ < 1 , increases 

for h >1 and remains zero for h„ = 1 . This increase and de- 

W 

1/2 

crease in '' is more pronoxmced as we go away from ^^=1* 

The similar trend is observed in vs. h^ for the Prandtl 

number of order unity as shown in Fig- 3.7, The last figure 
also displays W^/ for high and it is seen that mag- 
nitude of decreases as m increases, lastly, our re- 

sults (3-43) and (3.66) for the skin friction for low and order 
unity Prandtl numbers respectively differs by v. numerical factor 
■of 1.02. and tbe_ former is displayed in Fig. 3,8. 

Thus we see that simple analysis of the type described 
here can be vary useful in showing with sufficient accuracy, 
the analytical dependence of the boundary layer characteristics 
on various parameters. 



CHAPTER 4 


HEAT TRAITS EE R IH FALKITER-SKAH ■ 
BOOTDARY LAYER ELO?/S 


4. 1 Introduction : 

In previous Chapter 5s we obtained the leading terms in 
the asymptotic expansions, corresponding to the two limiting 
•values of the Prandtl number of various quantities in the flow 
of a compressible fluid when pressure gradient is zero. In 
general, it is found that the zeroth order equations are coupled 
and nonlinear but those of higher order are coupled and linear. 
They all lead to two point boundary value problems with missing 
boundary?' conditions at either point (determined by matching). 
Therefore, the process of obtaining a few higher order terms 
in the perturbation expansion is quite labourious and involved. 
These higher order terms, as a matter of fact represent small 
corrections to the previously obtained terms. Therefore, it 
is desirable' to develop at least some of them, in order to 
estimate theoretically the range of their validity. To accom- 
plish the above, we consider here in this chapter, the similar 
two dimensional or axisymmetricfflow with pressure gradient at 
limiting values of the Prandtl number for the case of constant 
property fluids. In addition to limiting Prandtl number, the 
case when Prandtl number is of order unity is also studied. Heat 
transfer in Elakner-Skan flows has been studied by various 
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authors (as described in Chapter 1). Invariably all of them 
have either neglected the pressure gradient or dissipation. 

4.2 Governing Equations: 

The governing equations for the first order, steady 
plane or axisyinmetric boundary layer flows of a constant property 
fluid as described in Sec. 2.5^1 are 

(r^Ui)g + = 0 (4.1a) 

b.|(u.|) - (4.1b) 

L.|(ii-.|) - (j~ ~ ® "^Is J (4.1c) 

The last equation (4.1c) may be written in a more convenient 
fashion in terms of total enthalpy g = (h^ + E u^/2)/gg-. 

Adding Eu^ times (4.1b) to (4.1c) we get, 

i'.i(g) - cr^ ^ -cr'^) (^-t^in)]! (4.id) 

and the corresponding boundary conditions are 

H = 0, u = 0 = V. , g. = or g'(0) = 0 

^ ^ (4.2) 

N = ~, U.J - 1, g ^ 1. 

Introducing the stream function defined by 

= fs = ; 

the equation of continuity (4.1a) is autonatically satisfied 
and if we assume the following form cfq^^ and g 

.*•' ' ■ ' . 



with 


6 5 

'Lp^Cs; M) = f(K|), g(B, S) = g(r|) . (4.3) 

S / 

1^= f ds, ^ = r3H/(2^)^/2 

0 

the equation (4.113) and (4. Id) reduce to the following total 
differential equations 

ft! H ^ to 

f + ff + p (1 - fl = 0 (4.4) 

g" +crfg' + C(cr- 1) (f’f" )‘ = 0, (4.5) 

Also the Doundary conditions (4.2) become 

f(0) = 0 = f’(0), f’(°°) = 1 (4.6a.) 

g(0) = or g’(0) = 0, g(“) = 1. (4.6b) 

Here, 

p = (dU.^/d^) is-- t-be-preasure gradient para- 

is the p.ressure gra.dient parameter and the coefficient 

0 - 1 ) 

Q r= 

2 + (^- 1 ) 1 ^ 

should be constants for the similarity to exist (M^ is the Mach 
number at the edge of boundary layer). The quantity G is indeed 
constant in the following three flows: (1) no pressure gra- 
dient, i.e., M = constant, (ii) incompressible M = 0 ■ 

(iii) high Mach number 1/M^ -► 0. However, in general G is a 
function of s. In our analysis we have retained G and treated 
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it as a constant. In the absence of more general results - i.e., 
without the use of similarity -- one may calculate C using local 
value of Mg and hope that this will not lead to significeut 
discripancA' - particularly for the flows in which the proper- 
ties do not vary rapidly. 

The momentum equation (4.4) is the falkner-Skan Bound- 
ary layer equation, and is independent of g and hence can he 
solved once and for all* 

4 . 3 Solution for Falkner-Skan Equation : 

The momentum equation (4.4) along with boundary condi- 
tions (4.6a) is the well-known Falkner-Skan (1930) equation. 

The solutions to this equation have been obtained by Hartree 
(1937) and recently improved by Smith (1954) for various values 
of |3 in the range (= -0. 198838) < p < 2, However, the solu- 
tions can be obtained for all p > 0. The mathematical exist- 
ence of solution of (4.4) subject to (4t6a) has also been proved 
by Tfeyl (1942) for 0 _< p _< 1/2. The uniqueness of solutions 
assuming they exist, has been established by Coppel (i960) 
for p > 0. 

Further, for p < 0 the solutions are not unique and for 
every p in a certain range of values, it is possible that, 
there is a family of solutions (i.e., a continuous range of 
value of f (o)) of (4.4) with f(o) = 0 = f (0) and f -* 1 as 
To exhibit this property we first introduce a new 
dependent variable 
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y(yj) = 1-f'(Kj) or 
in (4.4), v^e get, 


f = yj ~ a + f y(x) dx 
' 1 


y + (vj _ a + / y(x) dx ) y - py (2-y) = 0. 


(4.7a) 


(4.7b) 


Now ciS we approach the edge of boundary layer, ^ becomes 
la,rge and y becomes smaller, thereby, 


/ y(x) dx << '*'] 

.yv ' 


and y(if|) << 2 


the equation (4.7b) degenerates to 

y + (/]- a) y - 2p y = 0 (4.8) 

or y'--a^(Y|_ a)“^^"^ exp( ( l/2) ( '>^ - a)^) + ( vj _ a)^^(4.9) 

where and Ug are undetermined constants. 


For p > 0, it is clearly necessary to take Ug = 0, if 
the boimdary condition y ( ■’^•* °°) = 0 is to be satisfied, and 
then the solution is unique. However, for p < 0 the boundary 
condition y(r| — o®) = 0 does not require to be zero. This 
nonvanishing of introduces a lack of uniqueness in the solu- 
tion and in fact there exists a family of solutions of equa- 
tion (4.4) satisfying (4.6a), To select an appropriate solu- 
tion^ of equation (4.1) satisfying (4.6a) various arguments 
have been put forward. One is to suggest that the algebraic 
terms in (4.9), which introduces nonuniqueness arises essen- 
tially from the inertia terms and is therefore associated with 
the inviscid flow. So, it cannot arise with in the boundary 
layer and must be due to the boundary condition at the outer 



68 ' 


edge. Another argument is based on the consideration of dis- 

I 

placement thickness, i.e. , if f -» 1 algebraically and p 2 -'I? 
the displacement thicknes'- is infinite. Hence for p < 0, the 
condition at infinity to be used (as by Hartree 1937 and Smith 
1954-) is 

as f -» 1 ( 4 ^. -lO) 

exponentially from below. In numerical work it is not easy to 
decide whether the asymptotic behaviour ' is purely exponential 
or contains a grain of a.lgebra. However, it is reasonable to 
expect that with respect to the neighbouring solutions the 
exponential approach to be the most rapid. 

It is shown by Hartiee (1937) that there are no solutions 
of (4-. 4) subject to (4-6a) and ( 4 , 10 ) if p < Stewartson 

( 1954 ) found the duality of solutions for 0 > p > Pg, one of 

ff f« 

mhich has f (,0) < 0 and other f (O) > 0; both exhibit? the 
exponential behaviour ( 4 . 10) at infinity, Further, for the 
values of p slightly greater than p^ 

f” ( 0 ) a (p - Pg)^'^^, 

If 

At p = Pg, the skin friction f (0) vanishes. Experiments at 
high, but finite, values of Reynolds number indicate that the 
separation takes place at or very near the point of vanishing 
skin friction. 

We now consider the solution to the Falkner-Skan equation 
as being known. We will be using them later for small and large '''1, 
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For small Vj , it is 

f(v|) = a YjV2 - 13/|V3 + (2p-l) >|V5! 

+ (4-6p) ^av|^/6! - 2(2-3§)p^ vjV?! 

+ (l-2p) (11-IOp) a^ yj^'/8! +0 ( (4.11) 

and for large is 

f(^ ) = ''''] -a-' + f* ("^ ) , f*('^)s:0 ( ■'j (4.12) 

y V —CO 

Here a = f (0), a the displacement thickness and Tj means 
the exponentially small terms in the limit The valnes 

of a and a for Pg £ P ^ 2 are taken from Smith (1954). 

4 • 4 Solution for high Prandtl Humher : 

Following section (3.3) we now study the energy equation 
(4.5) subject to (4.6b) at high by the method of matched 
asymptotic expansions, further, as shown later, it is necessary 

If 

to stuc3.y the ■’problem seperately for the ordinary point (f (0)?^0, 

ft 

p p ) and at separa.tion point (f (O) = 0, p = p,^). This is 

due to the fact that near P = P„, the dominant term in the 

s 

5 2 

inner solution for f(''^) is of order of rj compared to '/I at 

ordinary point. Thus the stretching required to form the 

inner variable for the energy equation (4.5) is of different 

order in the above two cases: p * p„ and p p_. Lastly, it 

s s 

is shown that a logarithmic singularity is observed when 
p = Pg in the solution. 
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4 . 4.1 Solution at Ordinary point (p f Pg) : 

a) Outer Limit: The outer limit i.s defined to be £( = 0-“ ^ 0 , 

with fixed; and assume the follov/ing outer expansion 

g = 2 g (4.13) 

n=0 

Substituting the expansion (4. 13 ) in ( 4 . 5 ) and collecting the 
terms of various powers of £<, we get, 



^ f n ? 

- C(f f ) 


II 

0, fgg = 0., 


fg3 = 

tf ^ r fi t 

- go + C!(f f ) 

( 4 . 14) 

Tg^ = 

ff 

5 ^ 4 > 5 0 * 



The solutions to first four equations in (4,14), which satisfy 
the bounda,ry conditions at infinity are, 

So = 1 + c r ( f f ) /f , dvj 

! O 

= 1 + C (f '^-l )/2 + C / J/f\dy] 

V ’ ( 4 . 15 ) 

g-j ” 0, g2 ~ 0,1. 

g3 = - C / (J/f)'/f dv^ 
where J(v| ) = f''^ + p f‘ (1 _ f'^) ^ 

Now this solution for g shows that the total enthalpy g 
increases as -* 0. Thus this solution not only fails to 
satisfy the bo\nidary condition at the wall, but is singular 
there. To study the detailed behaviour of this singularity 
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near tlie wall we first write the inner expansion of the outer 
solution. The behavior of the integrand,.. J/f as 0 is 
(using 4- 11 ) , 

j/f is 2a/Yj^ - l6p/(3>|)+ , (4.16a) 

'Thus (4.15) indicates that the total enthalpy is singular 

•= oo 

-4 

like and In as 0. Now to make the integral ^ j/f d>| 
bounded we ha-^^e to add and substract singular part (4.16a), 
Thus the integral 

J/f - 2aVy| + 16P/3 y( drj 

is rcgultu; near ^ =' 0, but this integral diverges at upper 
limits. To avoid this difficultj^ we can combine the second 
and third terms of (4.16) binomlally a,s ^ 0, which is in 

some sense equivalent to writing 

J/f ^2a/-rj^ - l6p/(3 (i|+1)) ♦ (4.16b) 

The inner expansion of the outer solution (4.15) as 0 may 
be written as 


g 


1-C/2 + C Ip (0) + 2Ga/r| +{l6pC/3) lnv| 
. . . + 0(€^ lni| ) ^ 


+ • 


(4.17) 


where 


Ip(y|) = / d/f > 2 a/ + I6p/(3rj^ + 3>| ) drj 
is regular nea,r = 0. 

b. Inner Solution: It is now clear from (4.17) that a different 
inner limit is needed for small V| . from an order of magnitude 
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oiia.lysis we introduce the inner variables 

y = ^ /e , , e = eg (4. 18) 

and study the limit as e -» 0 with end 0 fixed. I'rom (4. 18) 
and the inner limit of the outer solution (4- 17), 3 it is 
obvious that the logarithm^' should" appear in the inner 
expansion 

e = So S2n)- ( + 

Substitutiny (4.18), (4-19) a-nd (4. 11) in (4.5) and collect- 
ing the coefficients of like powers of 6; we get 

Q(0^) + Ca^ = 0 

Q(e.^) = 0 

QCsp - jf- 5 ® y = ° 

Q(e3) - 5?- y = 0 

9.(64) - ^f- y^ Sgy + I cp2 = 0 

QCij) - 3^ y^jy = 0 (“^-so) 

Q(e6) - ^ y^ 6^y + (2P-1) I? y5 e^y 
+ ^ (2p-1) a^y^ - Ca^ = 0 

Q(9^) _ Q 

where Q is a differential operator defined as 
Q()= d^/dy^ + ~ y^ d/dy . 
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Each of the above equations is linear differential equation 
of second order. In present case each of the above can be solved 
in terms of gamma and related gamma functions. However, as the 
solutions of the higher 03?der equations require the solutions 
of the subsequent lower order equations, it is very laborious 
to solve many higher order equations. The solution to the first 
seven equations is given below 


0^ = - 12 Ca2 (1/3, x)/(36a)2/5 

02 = Mg + (5/3,x)/(48a5)''/5 „ pO(6N^( 2/3 ,x) 

- (1/3, x) + M^(5/3,x)/5)/2 

65 = 1 ^ + pb., (5/3,x)/(48a5)^/^ 

0^ = - pbg (5/3,x)/(48a5)V3_p2^^ ^ 

+ p2c (6/a^)''/5 N,^^^(i/3„x)/16 

+ (3,x)/l44 -H,^/3(5/3 ,x)/40 

(4.21) 

- 314 ^/ 3 ( 2 / 3 , x)/4 - 114 / 3 ( 2 , x)/4-H^/3(1,x)) 

83 = M3 + pb3 (5/3,»)/(48a5)V5 (5,x)/(l6a2) 

^ = \ + Pb 3 (5/3 x)/(48a5)''/3 + :• (3^x)/(l6a2) 

+ (l3/3,x)/(32(36a''^)'’/^) 

+ (1-2p)b^ (7/3,x)/(lO(36a)''/5) 

+ 2C (a^/6)’’/^ N _^/3 ( 1 / 3 , x) 

5 
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_N^^^(1/3,X)/10 + N5/5(7/5,x)/70)-Cp^(l2F5/3(l,x) 
-36ll3/3(7/3,x)/7 - 9 Nq/ 3(2/3 ,x)/2 + 313/3(2, x) 

-918/3^'’^/^’^^/''° N^^/3(1/3,x)/4- 3N^ ^/3(5/3 ,x)/20 

+N^^/3(3,x)/ 12 - N^^/3(l3/3,x)/52)/(3072a®)'’/^. 

In writing down the above solutions we have used the relation 

y 2 

J z” exp(-a2V3^' ) / exp (et^/r! )dt dz 
0 o 

= N^^^^^2.j.)/r[^^'^^'^2)/r,ayVr!] /(rfm) 

X 1 ^ 

vbere iT (a,x) = / t^N(a,t)dt, ^(ajx) = a/ e ^(l_t)*^“ dt (4.22) 

0 o 

are certain related gamiBa functions (see Appendix A). Further, 

IKI 

^(a,x)=/ e”"*" t‘^“'^ dt is "^he well-known incomnlete gamma function. 

and X = (a/3')y^, = a^ + 2h^ (l/3,x)/(36a)'''^^,r=0, 1 . „ , (4.23) 

and a. and h^ are constants of integration to be determined from 
r r 

the inner boundary conditions and by matching. 

c ResrCLts; Yarious unknown constants in (4-21) and (4-18) will 
now be determined by the use of matching principle. 

A. Insulated y^'all Case: Here the inner boundary condition g (O) = 0 
requir'es b^ = b^= ...= 0 and the matching principle gives the 
recovery factor 

r-f = (2g^-2+C)/2 

= 4(aV6)'''^^ n(-l/3,1/3) -32p ln0"/9 

+ 21^(0) - p m^ - 32p In (a/6)/9 
+ 2p^(6/a'^)^/^Cr‘^'^^ [(n(7/3,l/3)/l6+n(7/3,3)/144 
-■nil/^r5/3)/A0-3n{A/3,2/-5)/A+n{A./3,2)/A 
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+ n(l/3,1) [_4(aV6)''/^n(-1/3,l/3) 

+ (1-2p)(aV6)''/^ (-5n(2/3,4/3) + n(5/3,l/3)/5 

-n(5/3,7/3)/35) - P^(1 2n(5/3 , 1 )-36n(5/3, 7/3)/7 

-9n(8/3,2/3)/2 + 3n(8/3,2)-9n(8/3, 10/3)/l0 

+n(ll/3,1/3)/4 - 3N( 11/3, 5/3 )/20 + n(ll/3,3)/l2 

-n(11/3,7/3)/52/(384a®)''/^ ]+ 0(cr'' Incr). 

Where m. = / 6N(2/3,x)-xN(1/3,x)+x1T(5/3,x)/5-32/(x‘^+x) dx 
0 

and n(m,a) = ^^(^5 °°) • This can be evaluated in terms of gamma 

function from (4.22) by changing the order of integration. 

For the case of no pressure gradient p = 0, the above 
resiflt (4-24) becomes, 

r.= 1.922 - 1.341 + 0.468 0-“^/^ + .. (4.25) 

X •* 

Here in the above expression, the first term is due to Meksyn 
(1961) and first two due to Narasimha and Vasantha (1966). 

B. Heat Transfer Case: In this case the inner boundary condition 
g(0) = g^ requires a^ = a^ = 0, ^2 = g^, a^ = a^ = . . . = 0. The 
matching principle gives the heat transfer rate at the wall as 

g'(0) = (g^ - g^) I) (p,o^) 

where 

I)(p,cr) = 0.61627 (acr)''/^ - 0.09434 p/a 

-0. 03303(pVa'^) + [o. 013695 (2p -1) 

-0.020966 p^/a^] + 0( . (4.26) 
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4.4*2 Solution at Point of Vanishing Skin Friction: 

At the point of vanishing skin friction, the a.syinptotic 
relations just derived in ( 4 . 4 , 1 ) do not apply, as explained 
in the beginning of this section, I'he correct expansions can 
again he derived hy the present method. 


a) Outer limit: Following, section ( 4 . 4 . 1 ) we define the outer 
limit as CT with -q fixed, writing 

g = 2 g^a~“^ (4.27) 

in ( 4 . 5 ), we get 

fgo = - C(f f ) 

+ C(fVf") ’ (4.28) 

fgp = - , r = 2,3, .... 

Ihe solutions to the first two equations which satisfy the 
boundary condition at infinity are 


00 

go = ^ ^ ^ + 0 f j/f dy| 

. ^ (4. 29) 

§1 = ■ ^ -(j ^ ' 

Ihese solutions also fail^ at. satisfy the boundary condition ■ 
at the wall. To write the inner expansion of this outer solution 
( 4 . 35 ), we first evaluate the integrands from ( 4 . 11 ) as 0 


J/ f - 9 P A| 

( J/f) Vf - - 


(2+3p) p^ ^ V 4 ! + 0( ) 

+ (51 + 171 P) ^ + 0( r|5). 
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Employing, the same arguments as in (4.4. la), the inner 
expansion of outer solution ( 4 . 35 ) is seen to he 

g(Y^ - 0) = + a9|3C lny| + VC + CI^(O) + 

t cr"”^ [+(51-17ip)pc Inv^ /7 - 0 X 2 ( 0 ) _ 

+ 0(cr"^lnY|) (4.30) 

where 

00 

I^(^ ) = / J/f + 9pAj dy| 

12(1 ^ Vf + 54/^5 _ 

are regular near ’^ = 0 . 

h. Inner Solution: To get the inner solution, we first need the 

proper inner variables, Erom. an order of magnitude analysis 

1 /4 

the inner independent variable is found to be 2 ^ . 

Ho?;ever, if z is introduced in ( 4 . 36 ), then the leading term 
starts with lms| and this gives the proper inner va,riables as 

i e = g/incr . (4.31) 

we now define the inner limit as CT — “ with 6 and z fixed, and 
writing 

" = Jo <® 2 n^ (^*52) 

Ql(%) = 0 

Q.,(e,,) + -^ = 0 


we get 
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Ci^ce^) - ^ (2-3P) e^- = 0 

Q^(e3) - ^ (2-3P) pV 9^g + 560(2+3p) p^ z®/6! 

- 3Cp^z^/2 = 0 , 

Here is a differential operator, defined as 

-3 

= d^/dz^ - d/dz . 

The solutions to above equations are 

0 = M 

0 oo 

+ 3Cp (3/4, x) 

®2 " M 22 + to(2-3p) (24/-P)'^'^'^ i (|, X)/140 

'e^ > -^Cp|(2-3P) [SN^ (7/4, x)- 3/7 ^ 2 ( 3 / 4 , x) 

+ 9/77 ^2 (11/4, x)]| / 5 - 3 C p Nic(V 4 ,x) 

+ b,C2-3B) C2-4/-6)'V(9/4,'^)/H0 

where ' ■ 

X = -pgV4! 

Mrr = +[^r' ^24/P) i (l/4 , x) / 4 

a and h are constants of integration, 
r r 

c. Results: ' 

a) Insulated wall case: In this case h^ = "b^. - , . . = 0 and 
the matching the inner and outer solutions we get the 
enthalpy recovery factor is 
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= - 9p ln(r/2 - p m2 - 9ph(-p/24) + 21^(0) (4.35) 

- (51-171P) per"'' lncr /14 

+ [6 + Cp (2-3p)m3/5 + p(51-171p) ln(-p/24y48 

- 2.12 (0)] cr"'' + o( cr 2 in<r ) 

00 

where = f N(3/4, x) - 3/(4x + 4) dx 

0 

^00 

m = f 8xI(7/4, x) - 3x^ N(3/4, x)/7 + 9x%( 1l/4,x)/77 

- 187/14 + — — — - dx . 

28(x + 1) 

b) Heat transfer case: Here the inner boundary condition 
g(0) = requires = a, a^ = g^, = a^ = 0, and 

the ^latching principle gives the heat transfer rate at the wall 

g' (0) = (g^ - g^) D (pg, cr) 

where 



At 6=8 = -0.198838 the above result reduces to 

^ '^s 

D(pg ,cr) = 0.33285 0.007716 + 0( <7-"'^/^) . 

... (4.36) 
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4 . 5 Solution at. Ijow Pra.ndtl Number : 

In this section; we study the solution of the energy 
equation (4.5), along with (4.6h) for small values of <r . The 
procedure follow'ed is similar to that of (4.4. 1)* 



The solutions to first six equations (4.38) are 
% = Of' 2/2 + 0^^ + d^, 

S1 = °iT^ 

gj = - - rj2 a/2 + /' (i| - x) f*(x) dx | 

- C /^('vl - x) J(x) dx + + dgs 

° Y] i 

g^ = - a/2 + /'( x) f*(x) dxj + 

+ C^vj + d^ 
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g/ 


, 5 1 1 ^2 

C( ^1^/6 - yra/2) f J(x) 6.x - G f - xa)J(x)dx 
-Of (xv3 - x^a/2) J(x) dx + C / (v^ - x) f (x) f 

n r\ •' 


J(t) dt dx + C 


Y|^/40 - y^^a/8 + rj^aVs 


+ ( T^/5 - ^^a/2) S^f (x) dx 
+ /'(Tj-x) f*(x) •/ f*(t) dt dx 


/^(x^/6 - x^a/2) f* (x) dx 


1 * 1 
+ / ("y] -x) f (x) dx + + d 


- C, 


y|V6- yj^a/2 


4 


gc 


= C. 


2 v,3 


/6 


-^"/ 4.0 - a yy'/8 + 

+ ( - a’lVs) r f*(x) 


dx 


r1 


o 


+ / - x) f (x) / f (t) dt dx - 


0 
•X- 

where f is defined as 
* 


- f^(x^/6 ~ x^a/2) f*(x) dx 
o - 

+ -x) f*(x) dx 


- C|[^V6-or'f/2 


+ G 


5l 


+ d. 


f (t^ ) = a - + f( y| ) , f 5- 0( 7j" ) as v| 

and J(v|) = f"2 + pf’ (i^ f’2) ^ (4.40) 

It is seen from the above that the solutions become singular 
at largeV| and do not satisfy the boundarj'' condition at 
infinity. This singularity is similar to one encountered in 
Stokes solution in low Reynolds number flows. 
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4'. 5.2 Outer limit: It is clear from (4-55) that a different 
outer limit is needed for large Yj . From ar order of ma.gnitude 
analysis v^e introduce the outer variable ■ 

^ (4.41) 

and study the limit CP — 0 with ^ fixed. 

From, (4.12) for large , we have 
f('y|) =7y| _ a + 0( ■v|‘"°°) 

and the equation (4.5) ie present outer limit reduces to the 
follov?ing outer flow equation 

^ - a ) == o(cr“) ^ (4.42) 

This equation is correct to all order in or, i.e. , the error 
here is exponentially small, as explained in section (3.1.11). 
The outer equations of all orders may he obtained in a routine 
way from this simple equation. But this equation can be solved 
once and for all-. The solution to equation (4.52) which 
satisfies g('>°) = 1 is 

g = A + (1 - A) erf [( ^ - a(r'^^y^2 ] (4.43) 

where erf (x) is the error function defined as 

X 2 

erf (x) = / e""^ dt 

and A is a constant of integration, to- be expanded in powers 
1/2 

of cr '^ instead of solution, 

; . ; P : n=0 


Either the constants C’s or d's can be determined from the 
inner boundary?' conditions and the others alon^ with D’s are to 
be determined hj the application of the matching principle. 

4.5.3 Results; 


a. Insulated wall case: For a thermaJly insulated surface 
g (0) = 0, so = . . . = 0, and the matching 

princip3>e gives the enthalpy recovery factor 


■f(2mcr) / J(x)dx + 2<rS (a-x) J(x)dx 

^ oo ® oo 

+f(%<r^/2)f(a^-2 f f*(x) dx) / J(t) dt 

- ° . °o 3 “ 

+ 2 f J(x) / f(t) dt dx J + 2a'% / J(x) dx 


°° X 

+ / J(x) / (a-t) f(t) dt dx 

*^oo oo __ , 

+ f (a-x) f*(x) dx / J(t) dt I + (4.44) 

o o - 

b. Heat transfer Case: In this case we prescribe a constant 
wall enthalpy g(0) = to study the heat transfer at the surface. 
In order to satisfy this inner boundary condition we must have 
d = g,„, d., = do = . . . =0. Further, the matchinv gives the 
heat transfer rate at the wall 

g (0) 

whe re 

Here if dissipation is neglected i.e, g^ = 1, the first three 
terms in (4.45) are due to Coddard and /icrivos (1966). 


= (g^ E^J 33 (p,cr ) 


(2crA)'^/2 _^a-+ {2cr%)^/^ ( Af*(t) 


dt 


+ a^) + cT^ (- — f (a+x)f*(x) dx + 

£_ 7X TC ^ ^2 ci 


+ 0 ( ) * 


3m: 
(4.45) 
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4.6 Solution for Prandtl Nmiiber of Order Unity : 

To solve the enei^y equation (4.5) along -with the 
boundary conditions (4.6b), when Prandtl number is of ord-e-r- 
unity, vre use the method of inversion of variables developed 
by Meksyn (1961). The equations (4.4) and (4.5) are integrated 
respectively once and twice to give 

f ’ = y( v| , p) exp ( - P) ('4. 46) 

(4.47) 


g = 1 + B / ‘ exp ( - cr P) dy 

o 

A ^ 

^ C(<ri. 1) / exp (-(tF) dx J y.(t) exp ((cni2)F)dt 
o 0 


X 


where 


P 

y 


/■ f(x) dx 
0 

p p (1-f’2 (t) ) dt 


cl 


(4.47a) 


o 

f tf ? 


y^ = (1 + 


Here the constants of integration A and B are to be evaluated 
from the boundary conditions. 

t 

Por an insulated wall, the condition g (0) =0, requires 
S = 0 and using the boundary condition at infinity g(°°) = 1, 

the equation (4.47) becomes 

oo i: 

g = 1+ C(<r-i) / exp (-crp) dt / y.j exp ( (ci.2)P)dx^ (4. 48) 

o 

The re cove rj)" factor is 

oo „ _ 

r„ = 1+2(cr-l) / exp (-<rP) dt f y. exp ((<r -2) P)dx . 

^ o o 

(4.49) 
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Next, when the wall enthalpy is prescrihed, i.e, 
g(0) = = 1. In view of the boundary condition at infinity, 

the solution (4.47) becomes 


gC^n) = + (g^^-g^) / exp (-crP)dv^/ /exp(- (TF) dn 

0.0 
1 _ ^ — 
-C(.0~-1) / exp (-crF) dt / y. exp ((cr-2)F)dx, 

o 0 

(4.50) 


Heat transfer rate a.t the wall is 

g (0) = (g^-g^p/HT, HI = / exp (-crp)d^ ^ (4.51) 

o 

Nov; to evaluate the integrals (4.48) and (4.51) we 
need in terms of "F. -To obtain this relationship, we follow 
the usual procedure by assuming a power series of the type 
( 4 . 11 ) for f in terms of Vj . Now with this series (4.11), the 
expressions ( 4 , 47 a) become 

F = (a^V5!)- p v|'V4! - s.^ /S\ + 2ap(2-3p) ^Vv! 

-p2(4-6!3) + 0( ■)^^) (4.52) 

y = a-p^i + 2pa2 Tj^/3! _ pa( 1+6 p) y|V4! (4.52) 

+ p^(l+6p vj5/5! + 10p(a^(l+2p''^V6! + 0 (t<|^) 
y^exp(-2F) = a^ _ 3ap-rj+ 3p^7^^/2 -5a^(l-2§) ')^V6+ ” 
Inverting the first of the (4.52) we obtain the following power 



( 4 . 53) r = 3 or 4, accordingly as a ^ 0 or a =0. So, 

these two cases are considered rseparately. 
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4.6,1 Ordinary point Case a ^ oi In this case r = 3 and the 
relation (4.53) "becomes 

^] = (SP/a)''/^ + p( 6 F/a) 2 /V(l 2 a) + (8a^) 

46/3) V3 r 1 _ _ _ 

-- ^ 864a‘^ J 


V5 I ^ 

86 4a-^ 


Expressing in terms of E we get 

y. = - 313 (61&,2)V3 ^ 5p2(g^^g,^2/3/4 + ... (4.54) 

• * 

Substituting (4.54) in (4.49) and using the relation (3.47) we 

get the following expression for the recovery factor. 

r^ = 1 + 2 (cr“-i) [3(36a^)''/5 ^2/3) | ,^) 

-9P ?42l'i(l/3, 1, 4/3. ^)- 


p/8)(6/a-^) 


V3 P/l 


P(^) ^3925’^ (1/3, 4/3, 4/3,^^) 


+ 8 2F,,(2/3,4/3,5/3,^)/3+2l'.,(l,|,2,2^)J 0 ^^/^+ * *] 


For the case of no pressure gradient |3 = 0, the above result is 

r = 1+6(cr_.i)(36a"^)^/5 [7|) 2F., (5, |, | , ~)<F^/^+ •*: 


(4.55) 


Similarly the heat xtransfer result (4.51) reduces to 

HI = (g^ - g^)/g'(0) 

= 1.62265 (e/cr)^^^ - 0.248400 (p/a^/^) 

+ 0.12500 - {0.08764 pVa^ + 

+0.03606 (1-2P) ] (a"''.cJ-~^)^/^ + 0 (cr’’^) . 
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4.6.2 Case of separation point e. = 0: Por the present case 


when e = 0 , p = 13^, r = 4 end the equation ( 4 . 53 ) becomes 


"I = (24r/-fi) 


1/4 r 


1+(2-2p) P/1+0 + (3p+2)(202-f 


+ 145p) P^/77160 + .... 


Purther the expression for becomes 


(4.57) 


1 


= (-24p^P) 


1/2 r 


3/2 + (465 p - 64) P/82 + .. . 


(4.58) 


Sovi using ( 4 . 57 ) and ( 4 . 58 ) the result ( 4 , 49 ) for the recovery 


factor becomes 

r = 1 + 57.2653 2 ^ 1 '^''/'^’ » 5/4, -^)(1-^"^)+ .... 

(4.59) 

end the head: transfer result (4.51) become 

HT = 3.0043 + 0.06965 

-0.004899 + 0( 5/4)^ 


4,7 Discu ss ion 

The main results of interest in this Chapter are ( 4 - 26 , 
4 . 36 ), ( 4 , 45 ) and ( 4 . 56 , 4 . 60 ) for the heat transfer and ( 4 . 24 , 
4 . 35 ), ( 4 . 44 ) and ( 4 . 55 , 4 . 59 ) for the recovery factor respectively 
for the causes of high, low and order unity values o.f Prandtl 
number. 

The heat transfer results (4-26, 4.56) for high CT" are 
shown in Pig. 4.1. The convergence of higher" series (4. 26 ) is 
good even up to cr'= 0.5 with -1 < p ^ 2- In the region p 4 -1 
the convergence of the series (4.26) becomes poor owing to the 
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v/as not always successful, a.n extension described ly Hardy (1949) 
and cailed by Him ttie (E, q) process gave very satisfactory results. 
In this method whose application is described in Appendix A-2, we 
must choose a suitable q, possibly by trial and error procedure 
(in conventional Eulerisation process as q is ta.ken as unity) v/hich 
gives the desired convergence. 

We can improve the convergence of either the low or the 
high T or both the series in the intermediate range. In the pre- 
sent work we have Eulerised the low series (as its convergence 
is poor com.pa,red to high results) and the results are shown 
in Eig. 3.4. Ihe results are extremely good in the sense that 
the gap is filled even for p = -0.16 and -0.198838. 

In the above para the low CT series is Eulerised in the 
region 0.1 < CJ” < 1. However, the series can be Eulerised even 
for much higher values of To illustrate this we again study 

the no pressure gradient case and compare the present results 
with those given by various authors in Eig. 4.5. If the first 
term in lov/ 'J ' series is used as given by Stewartaon (1964) there 
is an overestimate of D(0, 0.25) by 40% . If two terms are 
used as given by Morgan, Pipkin and Tfemer '(1958), there is an 
overestimate of E(0, 0.25) by 28%. No comparison is made with 
the three term result of Goddard and Acrivos (1966) as there is 
a numerical error, in their third term for p = 0. Further, the 
result of Sparrow (1955) » 

+ 0.82 0-'’/^), Cn- 0, 



90 


based upon the Karmsji-Pohlhausen method underestimate P(P,a~), 
throughout the range, by 6 % at small values of .0* and by 12% 

at 0 = 1. It is interesting to note that Eulerised low cJ- results 

pass through the exact results of Pohlhausen (1921) even for 
: t " ~ -^ 5 ^ an( 3 _ for.ir of this order overlap with only the one term 
in high-0~ result (4.26) due to Narasimha and Vasantha (1966). How- 
ever, to eulerise low series beyond this point, we need more 
number of terms to get useful results even with large values of q. 

Regarding the heat transfer results (4.56, 4.60) for 
Prandtl number of order unity we wish to say as follows. If we 
just expand D(|3,'T ) series (4.56, 4.60) binomially for large 
we arrive at results similar to (4.26, 4.36) obtained for large 
0", It does not mean that the heat transfer = P(p,(.T') (g^-g^) 
is same from the high and order unity values of the Prandtl 
number. In fa,ct the recovery enthalpy g^ are governed by diffe- 
rent expansions. However, it is fortuitous that if dissipation 
is neglected (g^ = 1), the results (4.56, 4.60) for rr':sx0(1 ) 
expanded binomially for large cJ~ gives the results same as (4.26, 
4.36) for large X - ff dissipation is neglected the results 
(4.56, 4.60) reduces to those by Merk (1959) and Evans (1967). 

Let us now study our result for the recovery factor (4. 24) 
and (4.24) for high and low Prandtl number. Again convergence 
of both the series becomes poor as p decreases. We have again 
improved the convergence of low CT" re suit. The comparison of 
present results with various availables results (only available 
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for (3 = 0) is made in fig. 4.6. Our result (4.25) for high -CT" 
gives r^ = 1.048 at O'" = 1. If the first term of the result 
due to Meksyn (1961) is used, it overestimates r^ hy 92% and 
if first two terms due to Narasimha and Vasantha (1966) is used 
there is an underestimate of r^ by 62% at CJ"= 1. The Eulerised 
results from low •-0“ series show that at CF =1, - 1.0019 (the 

exact result of Pohlha^sen 1921 is unity). If the first term in 
luwcT'series is used as given by Stewartson (1964), There is an 
underestimate r^ of 5% near 0~ = 0.25. If two terras of Morgan 
et.al (1958) are plotted there is an overestimate of r^ about 6^ 
near 0 = 0.25. It is interesting again to see from Eig. 4-6 ''n 
that the results from low and high -T' series fills the gap in the 
interm.ediate range. Similar comments apply to the results with 
pressure gradients. Lastly, our results (4.55) for Prendtl 
number of order unity he.s a graphical representation similar 
to that shown in Fig. 3. .7. 

From the above study, it is now clear that how the simple 
study presented here for two limiting cases: low and high values 
of Prandtl number, fills the gap in the intermediate region where 
Prandtl number is of order unity, provided a sufficiently large 
number of terms in each of the expansions (low and high 'IT ) is 
computed and the series are properly Eulerised. Thus, in some 
series, it is enough only to study the two limiting cases. 



CHAPTER 5 


SEOOITO ORDER EEEECTS ON UiMIMR BOUHDARY LAYER IH 

HOH SIMILAR PLOW OR A CONSTAN T P ROPERTY ELUID 

5*1 Introductio n : 

The Prandtl houndary layer theory provides the leading 
term in an asymptotic expansion for large Reynolds number R, 
a,iid does not account for what Rott a,nd Lenard (1959) ca.ll the 
secondary effects. These may be classified as those due to 
transverse curvature, longitudinal curvature, vorticity gra- 
dient, stagnation enthalpy gradient, velocity slip, tempera- 
ture gump a,nd displacement. These effects are of order 
compared to those of classical boundary layer theory and 
become^ significant where R is only moderately large. Van 
Pyke (1962a) gave a systematic method of improving the , cla-ssi 
cal boundary layer theory for these effects of order 
and is described in Chapter 2. A critica.1 comparison of 
this theory with various others is given in Chapter 1 . 

Van Dyke (1962b) has applied the theory so formula.ted 
to study the special cases of plane flow, ax:isymmetric flow 
and of cusped leading edge. He studied the lea.ding terms in 
the Blasius series expansion in the coordinate a.long the 
body. However, at present no general solution of the second 
order boundary layer theory even for similar flows is avail- 
able. The purpose: of this study is to provide a. general 



analysis of the G-ortler (1957) type for non similar flows 
with heaf tre.nsfer. In the spirit of v/hat has been said in 
Chapters 1 and 2, ?/e shall be studying the second order effects 
of transverse curvature, longitudinal curvature, vorticity 
gradient, totad enthalpy gradient and displacement in plane 
and axisArrninetric flows of a constant property fluid. 

5 * 2 Governing Equations : 

The governing equations for second order boundary la-yer 
theory are obtained from Navier-Stokes equations by the method 
of matched asymptotic expa,nsions. This method results in 

replacing the Naa/ier-Stokes equations by two sepa,ra.te sets of 

-l/2 

equations, one set is vadid in the inner region of order E 

- 1/2 

near the wall a,nd other outside this region of order R * . 

The solutions of these tv/o regions arc matched in an> overlap 
dom.ain. Par the present study of constant property fluid 
flows with small temperature changes, V7e can neglect dissi- 
pation terms in the energy equation (this being easily justi- 
fied). Now the governing equations of motion for the first 
and second order boundary layer theory as given in Chapter 2 
are as follows. 

first order boundary layer equations: 
continuity 

= 0 (5.1a) 


momentum 



34 - 


ene rgy 

= 0 

bound cry conditions 

u^(s, 0) = 0 = v^(s, 0), Ti^(s, 0) = h, 
Hatching conditions 


w 


u^(s, K) .--U^(s, 0) = 


h^(s, N) S^(0) 5 . 

The second order boundary layer equations: 
continuity 

(r^Ug)^ + (r^Vg)!^ = - Kr^(Nv^)j^ 


as lif -♦ °° 


(5 


(5 


(5 


- D(r^”' cos 0 Nu^)g - (r^”"^ cos 0 Hv^) 

momentum 

^ ^"l ^ ~ ^2NU ~ cos 0/r) ^-j 2 ;j 

+ K(Nu^u^g ^ ■^1N ■ 

. + (KU^ + K / u5 - u? (s, N) dn)^ 

. \ N ‘ ^ ® 

+ r^ B (0) V (S,0) + (U^ T> ) 




(5 


energy 


-1 


” CT* '•■‘2KK' 


hoTVTTvT {o oos 0/^r} ^ 

+ + Nu^h^g) (5 

boundary conditions 

UgCsjO) = 0 = VgCs, 0) , hgCs, 0)=0 (5 

matching conditions 

UgCs,!) ^ -KNU^ + r^KB’(O) + Ug 

hgCs,!)^ i^(s, N) s’( 0 ) 


as u -♦ oo 


(5 


.1c) 

. 1d) 

. 1e) 

(5.2a) 

.2b) 

.2 c) 

.2d) 

.2e) 
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where 

VgCs, 0) = ^liEi (v^ - 

L^( operator) = 9/6 s + 9/8N 


and 

= U^(s,0) and Ug = UgCs, 0) . 

5 . 3 Bonndary Layer Characteristics : 


To study the detedled "behaviour of the houndary layer 
it is also essential to study the second ns der effects on 
various boundary layer characteristics such as, the shear, 
stress and heat transfer at the wall, the displacement and 
momentum thicknesses; these quantities a-re defined a,s follows 
shear stress: 

'r = u^ (s, 0) 

= t R-'* ‘ 

heat transfer: 


- cr"^ R h^(s, 0) 


,-1/2 


• 1 


R~ '^ *- q^ + R 'qg + " * 
displacement thickness: 

* -5-1/2 


R 


/ ^u(s, - u(s, h) J /U^ 


dN 


= R”'*/^ + R"'' 6* + 


momentum th iefcne ss : 

r ~'*/2 


f |^u(s, ~)-u(s,N) j u(s,H)/U" 


(5.4) 


(5.5) 


dU 


- r -'‘ I /2 + j {-1 + 




(5.6) 
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Now, the suhstitution, of the inner expansions (2.17) in the 
above a,nd the collection the coefficients of various powers of 
R give the boundary layer chara,ct eristics for successive appro- 
ximations. Thus, the leafing terms gives the characteristics 
for the first order boundary layer as 

= -C3r-~ 0) 

00 ^ 

6* = / (U^-u^)U^ dN 

0 

00 

^ = / (U^-u^)u^/U^ dN. 

The next higher terra.s give the characteristics for the second 

order boundary layer as 

t 2 ~ 

62 = / (^2 - KNU^ + r^B^CO)! - U2)/U^ dIT (5.8) 

o 

A 2 = /[(Ug - + r3B^(0)N~U2)uyu^ 

+ (U^ - U^)U2/U^ jdU ^ 

The various boundary layer paravmeters '*l^( ^ f 2 ^ ^ ? ^^^2 

are all of order unity. 
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5 . 4 Choice of Independent Variables; 

In order to study the houndary layer flows one generally 
resorts to similarity solutions (i.e., the partial differential 
equations are reduced to total equations). However, in majority 
of the engineering problems with diverse applications, it is a 
rare occurance indeed when all the conditions for the simi- 
larity are satisfied. In literature various methods hhve been 
used to study the non similar flows, a few of them a, re men- 
tioned below. 

In the first of these methods instead of solving the 
boundary layer equations, we adopt approximate techniques such 
as integral methods etc., see, Devan 1965 and Devan a.nd 
Obera.i 1964 (the references in this para are of those works 
which pertain specifically to the second order boundary la.yer ). 
L second type of non similar calculation employs a strictly 
numerical approach. In this method the full pa^rtial diffe- 
rential equations of the problem are solved numerically (see, 
Davis and ZLugge-Iotz 1964b and Dannelop and Plugge -Lot z 1966). 
There are va.rious difficulties with this type of numerical 
solution. First is that a new calculation is to he performed 
for each set of data in a specific prohlem and one cannot 
draw general conclusions aLout the behaviour of the solution 
with respect to parameters. Secondly, the niimerical methods 
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are expenst^e. The third method is to use some sort of local 
similarity. In this class vre can mention the method of 
Blasius (1908), G-ortler (1957) and that of Meksyn (I96l) 
which was later developed by Merk (1959). Bush (1964) later 
incorporated the important terms of corrections, neglected 
by Merk (1959). This later method studied by Meksyn, Merk 
and Bush will be cited below as LiMB. 

The main shortcoming of Blasius series method is well 
known, namely its poor convergence for slender body shapes. 

The reason for this is as follows. In Blasius series method 
the outer bounds, ry condition (i.e, , the U^(s) approeiches 
a sym.pt otic ally at large distances from, the wall) is satisfied 
by expanding .U..j(s) in a power series of s end adjusting each 
term, of Blasius series so as to assume the value of the 
corresponding term of the outer velocity series asymptotically. 
However, if U.^(s) does not happen to ^a polynomial of low 
degree (as is the case in slender body flow) , the outer 
boundary condition is satisfied with sufficient accuracy only 
if an appropriate large number of terms of Blasius series 
are talcen, and the rapidity of convergence of Blasius series 
will depend to a considerable extent on the convergence of the 
power series of U.^(s). 

The later methods of Gbrtler and MB have the advantage 
that the leading t erms of t heir series satisfy the outer 
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boundary condition exactly at all cross-sections a-long the 
wall. Therefore all further terms of the series give correc- 
tions only to the inner part of the boiands-ry layer. Secondly, 
the well known simila,rity solutions form^" the leading terms 
ill the above method. But in these methods the leading term, 
itself, gives a good approxima.t ion for the boundary layer 
flow from origin to down stream crosssections surprisingly far 
a,way. Binally, the above methods of G-brtler and MliB are 
not only applicable to all cases of cylinders with round noses 
(cases to which Blasius series also applies) but nay also be 
applied to cylinders with forward a.peX' line. However, very 
near to the point of zero skin friction none of the above method 
gives satisfactory results, nevertheless, these methods can 
be used to calculate the flows ahead of separation. The Ccilcu- 
l8,tions can later be continued with the aid of suitable numeri- 
cal technique. 

The underlying principle in these methods of Gbrtler 
and MMB is the following, *if the boundary layer is at all 
tines very nearly described by a, similar solution, then the 
direct effect of the non similar terms can be calcula-ted by 
expanding the full bounda-iy layer equations interns of a small 
parameter which measure the departure of the solution from the 
similar one'. The main difference between Gortler and MB® 
methods is in the choice of the expansion parameters. 
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G-ortler method the expansion parameter is a variable 
defined by 

< = U as 

0 

and in MMB method the expansion parameter is dp/dX where 

p = 2d In U^/d In^ , X = In . 

No compa.rison of the a.bove two methods is ava,ilable in the 
literature. The ninnerical solution of the resulting equations 
in MMB method is relatively much more involved due to the fact 
that the non- similar equations contain terms of the type 
8f^9p and a^f^dp , v = 0,1,2... . for example even to 
obte.in the first non similar term we need sufficiently accurate 
values of quantities like 3f^/ap , a^f^/ep, aX (f^ being the 
similarity solution) and the procedure for obtaining the succes 
sive terms will be much more involved. 

In view of the above we adopt the Gdrtler (1957) type 
method for our study of the second order bo-undr ry ia-yer theory. 
However, the Gortler (1957) variables ( » '^ ) are designed 
for two dimensional flow. In Schlichting (1962), it is men- 
tioned tha,t V. Saljnikov has extended the Gortler scheme to 
axisymmetric case for first order boundary layer. The work 
of V, Saljniko-^ is not available to the present author. However, 
it is fairly obvious that the following scheme (a combination 
of Gortler and Mangier transformations) works for both the 
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two-dimensional (.1 = 0 ) and the axi symmetric (,-] = 1 ) cases. 

We pose 

s 

^ = f ds, '^= ( 5 . 9 ) 

5 - 5 l irst Order Boundary layer : 

The equation of continuity ( 5 - la) for the first order 
boundary layer is satisfied by introducing the first order 
stream function > i.e. 5 

'Viij = '■Ks = (5.10) 

The stream function and the enthalpy h^ are now expressed 
in terms of our variables ^ defined by (5.9) in the form 
'\|/^(s, N) = r( 2 ^)f(^ ,y| ) ( 5 . 11 a) 

h., (s, N) = (h^ _ S^) ) + S^. ( 5 . 11 b) 

On the basis of (5.10) and (5-11) the velocity components are 

U.J = U.,f' ( 5 . 12 ) 

v^ =-U.,r^ jlf +2 ^f^ + (p + _ 1 . ) >^f']/f( 2 ^ ) ( 5 . 12 ) 

where 


Ag = ( 2 ^ /r) (dr/d^ ) 

= ( 2 dr/ds)/(U..r^^) / U..r^^ ds. 

o 

The first order boundary layer equations (5.1) take the forma 
momentum 

’ 2 ^ 


tf I ft 

f + ff + 


p (1 - f’^) = 2 ^(f’f'^ - f^ f") (5.13) 


cr”* s** + fg'- Af'g = 2 ^(f'g^ - f^ g ) ( 5 . 14 ) 


energy 
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boundary conditions 


f(^, 0) + (^, 0) = 0 

f’(e, .0) = 0, f'(^ , °°) = 1 

g(^ , 0) = 1, g( ^ , oo) = 0. 


(5.15) 
( 5 . 16 ) 


Here the dashes indicate the partial differentiation with 
respect to . 


The quantity p, defined by 


p(| ) = 2d In U^/ d In ^ 

= (2IJ^g/U^ r^^) / TJ^ ds (5.17) 

o 

is ca,lled the principal function by Gortler. This principal 
function is of fundamental importance in the structure of the 
solution. Thus., only the form of p(^ ) differs from case to 
case in the boundary layer problems (5.13) and (5.15). It is 
easily shown (see, Gortler 1957) that if two problems are simi- 
lar in the sense of Reynolds similarity they he.ve same principal 
function. Farther, the behaviour of the solution at singular 
point (say s = O) depends on the behaviour of the principal 
function at ^ = 0. For the special case of t wo dimensional 
flow with a s™, the principal function p reduces to the 
Hartree parameter of wedge flows. This principal function p 
is as meaningful for a non similar flow as for a similar flow, 
the main difference in t wo applications being tha.t while p 
remains constant in the stream direction when flow is similar 
it varies when the flow is not similar. 
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The principal function ) depends not only on the local 
values of the outer velocity U. and outer velocity gradient U. 
hut also on the history of the flow upstream from the cross- 

S o * 

section / r ^ ds. The method, of course, lequires the know- 
o 

ledge of the distribution of the main stream velocity outside 
the boundary layer. This miay be interms of gk formula, fitted to 
an experimental data or to obtain theoretically bj?- solving the 
outer flow. 


The quantity A.( ^ ) defined by 


A(^ ) 


-1 


2 ^ (V- ap 
2(VSP'^ “w 


U^r^3 


ds 


S 

f Ur 

0 


,23 


ds 


( 5 . 18 ) 


will be called the principal thermal function. It is interest- 
ing to note that A arises in the temperature problems in a 
manner similar to that in which p arises in the velocity prob- 
lem. further, comments similar to those mentioned in previous 
para for p apply to A-. In our present study we shall be study 
the ca,se when wall t emperature is prescribed to a constant value. 


We choose to study the situation where the velocity 
U^(s, 0) may be represented by the convergent power series in 
the interval 0 _< s s^ (where s = s^ denotes any suitable 
cross section not further down stream than that of a possible 
separation point of the boundary layer >$!. 
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tJ,{s, 0) = 5 ; 

. m=0 

r 21 

v/Taere x = J r ds. 


( 5 . 19) 


Further, the relationship between ^ and x is 


«; = 


2 ^3? ^r(m+1) 


m+'l r=0 


( 5 . 20 ) 


The form of the principal function p corresponding to 


( 5 . 19 ) is 


Z p. 

i=0 ^ 


(5.21) 


where = 2 ro/(ni+ 1 ) and depends on coefficients 7 \|y^of 

equation ( 5 . 19 ). 

How to study the first order boundary layer problem (5.13) 
to ( 5 . 16 ) we assume that f and g are represented by the following 


power series 


f(^,^ ) = + 5, + ... 

g('^ »"^ )= go(^ ) + ^ (■»!) + ^^g^^^('^) + .. . 

which converge in the interval 0 ^ Substituting (5.21), 

( 5 . 22 ) in equations (5.13) to ( 5 . 16 ) and collecting the coeffi- 
cients of the like powers of ^ we obtain the equations for successive 
approximations. Thus the terms independent of ^ are 

f + f f + p^(l_f ^) = 0 

0 O 0 ^0 0^ 


(5.22) 


cr-"! g" + :^go = 0 

f^( 0 ) = 0 = f^( 0 ), 

go(0) = 


fQ(°°) = 1 

go(“) = 


( 5 . 23 ) 


0. 
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The terms of order ^ are 

(l_f^2) 

Z^ig^^h = -3f^''^g’ 

f(0(o) = f'^^^(o) = = 0 

^(O) = 0 = ^(°°). 


(5.24) 


The terms of order ^ are 

= ^ 2 (^ 0 ^ - 1) + P^(2f^f’^^^~3f^’'^f"*'”'^ 

+ (pQ + ) 

= 2f'^''^g^‘^^ - - 5f^^^g' 

f(2)(o) = ^ ^ Q 

g(2)(o) = 0 = 


( 5 . 25 ) 


Here 

Y^(f) = f + f^f - 2(p^+r)f^f + (l+2r)fQf 

„ , , (5.26) 
Zp(g) = cr- g ■*■ ^oS - 2rf^g . 

It is to be noted that equations (5.23) governing the func- 
tions f^ and g^ are the same as those of an equivalent wedge of 
angle p^it. further, the' independent variable''^ on which f^ depends 
reduces to the variable of similar solutions when p(^) = const. 
(U^oCx®^). In general case our variable ''^ approaches the variable 
of similar flows for sufficiently small x or ^ . It can therefore 
be expected that the first order approximation f^ will be good not 
only in the immediate viscinity of the origin but even rather far 
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downstream. Moreover, the functions and satisfy the outer 
"boundary conditions exactly, thus all further terms f^^^ etc repre- 
sent only the correction to inner part of the boundary layer. 


To solve the equations f^, g^ we need to know while 
for g^^^ we need _ _ _ _ Thus for each set of 

data of a specific problem we have to solve the entire set of 
equations. However, we can take the advantage of the linearity 
of the higher order equations by rephrasing our problem in terms 
of some functions which are independent of the explicit data of 
a particular problem. Thus, for any fixed value of p^, it is possible 
to reduce the equations f^ l)(r, , Pq, P^, Pj^) to a linear combi- 
nations of functions not depending on p^ , Pg, Pjj_ and therefore 

for any fixed value of p^ the universal functions can be tabulated 
once and for all. The stated reduction is achieved by setting: for 

a fixed value of 8., 

'^0 



^ 1^11 ^ 2 ^ 2 * 


(5.27) 


Similarly for fixed value of p^ and cr~the corresponding energy 
equation may be splitted: into universal functions by setting 

^ (5.27a) 

( 2 ) 2 

'' = p., g^^ + p2g2- 

How the final form of the expansions (5.22) which leads to desired 
universal functions are 
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^ ^ t + ^ 2 ^ 3 )+ ••• 


g = go P^g^ + <^(P?gii + P2S2)+ ■■■- 
Thus, the final form for the terms of orders is 


(5-28) 


Vf^) = fo^-1 


i(gi) 




(5.29) 


f. ,(0) = f^(0) = f^(~) = 0 

g. |(0) = 0 = g.|(“). 


2 

and the terms of oilier ^ give the following two set of equations 

Z^Cgg) = -Sfggo 

t (5.30) 

f2(0) = f2(0) = f2(~) = 0 

g2(0) = g2(~) = 0, 

and = 2f’f^ - 3f.,f" + (p^ + 2) 

2i(g^^) - 2f^g^ - 3f.|g.| - 5f-j.|gQ 

f..(0) = f’ (0) = f’ (co) =0 


(5.31) 


g-ii(O) 


(“) = 0 - 


Now the above set of equations (5.29), (5.30), (5.31) 
etc, are independent of p^, P 2 . . . and can be solved once and 
for all for fixed values of p^ and CT" . The results may be 
tabulated for the universal functions. Knowing p.^, P 2 • • • 
specific problems,^ The results may be obtained immediately 
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from the tabulated -universal functions. 

•■I T- 

liirthef, our equations (5.29), (5.30), (5.31) etc. obtained 
for the plane and the axisymmetric flows are of the same form as 
those of Gortler (1957) for plane flows. Gortler (1957) studied 
the nonsimilar momentum equations in (5.23), (5.24), (5.25) etc. 
for pQ ~ ^ only. Thus, the validity of the Gortler series 

for negative is not known. Sparrow (1958) studied the energy 
equations of the two-dimensional nonsimilar flows for p = 1, and 
0"= 1. Hence, even the study of the above equations (5.29) to 
(5.31) is not exhaustive in the literature. In the present work, 
in addition to the study of the second order boundary layer problem 
some of the nonsimilar equations for the first order will also be 
studied in detail. 

The first equation in (5.23) for the first term of the first 
order boundary layer theory , is nonlinear and is well known falkner- 
Skan equation, whose solution is studied by various .authors. Some 
general remarks about the existence^ uniqueness of the solutions 

of this equation has already been made in Chapter 4. In the pre- 

II 

sent work the initial value f (0) is taken from the tables of 
Smith (1954) and solution is generated at a step size of 0.05. 

This solution is shown in Fig. 5*1. How the energy equation in 
( 5 .23 ) and the equations (5.24) are linear and const itute the 
two point bo-undary value prob'lems. The solution to above mentioned 
equations are obtained by integrating them (the details of pro- 
cedure are given in Appendix C ) by Ronge-Kutta-Gill method 



on IBM 7044 Computer at IIT Kanpur, The solutions are obtained 
for both the accelerated and decelerated flow in the range 
-.198838 < pQ < 2 and the Prandtl number of the fluid in the 
range 0.7 to 3. The results of the numerical solutions are 
displayed in Pigs. 5.2 to 5-4. 


First order boundary layer characteristics: 


Substituting our variables (5.9) and (5.11) into the first 
order boundary layer characteristics (5.18) we get 


f"(r, 0) 






(S^-h„)o-Uir 


3 = g '(5 . 0 ) 


+6 a + 


(5.31) 


(5.32) 


(5.33) 


6* U.rVr(2f) = / (1-/ )dn 

0 

= ^ ^*,1+ ••• 

oo 

A. U.rVy'(2^) = / (1-f )/ dn 

0 

= ^1.0 tf'l '^1,1'^ 

low, substituting the variables (5 -.25) and collecting the coeffi- 


oients of the like powers of ^ , the leading terms are 
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*1,0 = { 'I'l 

^1,0 = { ^0 'I’l 

while the next higher order terms are 

^1,1 = ^i' (0) 

q.1 1 = gi (0) 

’ i , (5.36) 

*1,1 = / - fi 41 = -fi(“) 

°° . t 

Ai,i = { fi(i-2fo) • 

These characteristics are also evaluated for various values of 6 

'^o 

and CT and are displayed in Pig. 5 . 5 and Pig. 5.6. 


5.6 Second Order Boundary Layer Theory ; 

A feature of the second order boundary layer equations (5-2) 

- — as a matter of fact, of the second and higher order terms in 
all perturbation schemes are linear and hence, allow superposi- 

tion. following, the suggestion of Rott and Lenard (1959) we can 
divide these higher order problem (5*2) into a number of simpler 
problems, each of which has a clear physical interpretation and 
then superpose. However, this division is some what arbitrary and 
it is possible to divide these second order effects in many ways 
(particularly those of vorticity and displacement). Van Dyke (1962c) 
and Oheng (1966) have discussed how various authors (mentioned in 
Chapter 1) have subdivided the second order terms. The way in 
which these effects are divided here, in the present study, is 
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"bhe same as those of Van Dyke (I962a)« The terms on the right 
hsnd side of equations (5.2) fall into the following categories. 

In momentum, equation (5.2a) the terms proportional to (j cos e/r) 
arise from the transverse curvature, K from the longitudinal cur- 
vature, Bj(0) represent the effect of vorticity in the hasi-. on 
coming stream and U 2 (s, 0) arise from the displacement speed, i.e, , 
the change induced in the speed at the outer edge of classical 
houndary layer. Further, "in the matching condition (5.2e) for 
the energy equation (5.2c), the term proportional to s’(0) re- 
presents the effect of stagnation enthalpy gradient in the on- 
coming stream. 


Thus, the second order problem (5.2) can be decomposed the 
into five constituents by writting 


h. 



(5.37) 


(5.38) 

= (0)h(''hs^ (0)h( . 

(5.39) 


Here, the superscripts t, 1, v, s and d identify respectively 
the contributions of the transverse curvature, the longitudinal 
curvature, external vorticity gradient, stagnation enthalpy 
gradient and the displacement. We now proceed to study each 
of these effects separately by applying the G-ortler type method. 

5-6.1 T ransverse Curvature Problem ; 

Keeping only such nonhomogenous terms in (5.2) as involve 
(d cos 0/r) gives, the problem for v^"^^ and (now • 
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onwards these superscripts will be dropped in this section, 
however, the superscripts will be used out side this section 
wherever necessary) is governed by 

] s l]i “ ^ 

^ -j (^2 ^ ^ ~''^2l!n!T ~ 9/r 

(^2 ^ "*"^2 ^^1 ^ ^2Nl!l ~ 9/r (5*40) 

U2(s, 0) = V2(s, 0) = U2(s, ~) = 0 
h2(s, 0) = 0 = h2(s, «>), 

The equation of continuity (5.40a) is satisfied by introduc- 
ing the second order stream function 2> defined by 

r^(u2+jNu. cos 9/r) 

' ^ (5.41) 

14/23 (■^p+Pv^ cos 9/r) * 

Now, the stream function "^2 enthalpy h2 are assumed to be 

of the form 


■^2(3, = ir(2^) F(^,'i ) 

h2(s, N) = • 


(5.42) 


To be consistent with the study of Sec. 5.5 (for the first order 
boundary layer) the relation (5.11) are used for and h^. 
Thus, the velocity components on the basis of (5.'I2), (5.41) 
and (5.42) are 

U2 = U^p' - 3f’N(2^)''/^ cos 9/r^ 


(5 . 43 a) 



(5.43b) 


= 


.-1, 


(5.44) 

f T 


^2 = jF+2(^F^ +(P+DAg-l) ^p'J /(2^ )''/^ 

+ [f+2^f^ +(P+jA 3 -l) if* J(o cos e/r) 

where J\ ^ = (2,^/r) (d^/dr) will be termed as the principal 
body shape function. Making use of (5.10), (5-11), (5.41), 

(5.42) the second order problem (5.40) reduce to 

i-> I t tf U ! t t t !» 

F +fF -2|3f F +f F+2je (f^ F -P f^ -f F^ + f P^ ) 

' (1+ if)f”+ff’-pi (l+f’^)+2 ^f*f^ 

+ 4 i ( f ^ f - f f ^ f ^ 

‘G +fG'~A-f'G+2F (f G‘-f’G )= -Fg'+Ap'g+2<f (f G'-f'G ) 

^ S S ^ ^ \ 

+ (-a-“A3«| f)g'-2^r^ (f g^ g') (5.45) 

F(^, 0) + 2^F^ (^, 0) = 0 

J'*(^ , 0) = 0; f’( ^,T| ) = i;^r| , i-oo. (5.46) 

0(^ , 0) = 0 = G( °°). (5.47) 

Here k^(^ ) is the transverse curvature function defined by 

kt(^ ) = (2^ )''A cos e/CU^r^'^^) (5.48) 

and the function 

/\.(^ ) = (2^/k^) (dk^/d^) (5.49) 

will be called the principal transverse curvature function. 

It is interesting to see that arises in the transverse 
curvature problem in the same manner p arises in classical 
boundary layer equations. 
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In order to study the above problem in detail, we 
assume (in addition to (5.23) for the first order boundary 
layer theory) the following expansions for F, G- and 

F(^ ,vj ) = Y| ) + ) + .. . 

G( ) + ^ ) + ••• (5.50) 

> = ht + •••. 

Substituting (5.50), (5.23) in (5-44) to (5.47) and collecting 


the coefficient of various powers of ^ , the terms independent 
of ^ are 

Zo(4°>) = -4°)g; + k„g-cr--'+ro)8o ( 5 ^ 5 ^ 

p(0)(0) = 0 = f’^^^(O), ^ ) = “ 

(}(o)(o) = 0 = 


and of order f are 


Y^(F 


(1)^ „ 


3f(^)F’'(0)+2(l + p^)f* ^^^f’ ^^^-f"^^^F^^^+2p^f^F' 


( ^ + v| f (0 ) ) f "( 1 )+5 yj f ( "I ) fj^+f^f ' ^ "* ^+3fpf ^ ^ ^ 

-2(2+Pq )1 1 (1+40] [<o+vifp 

- Pol ] 

[-(r'h'<''0r|f^g'O)-2v^4gOh3r]f(''>g; ] 


(-cr\+'nfo)gc 


(5.52) 
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(0) = 0 = ^(oo) , 

Now the above system of equations involve, in addition to • 

^2? •'* quantibies etc. However, we can rephrase 

the problem by introducing the proper universal functions, Por 
the first order problem we still have to use the relations (5,27) 


and (5.28); 

for 

the second order 

' problem; (5,51) and (5-52) we 

assume 





p(°) 

(^) 

= k . 
ot 



(j{0) 

■( vj ) 

o 

!l 

hh) 

(5.53) 

p(1) 

(^ ) 

= 

S', ( -^ ) + 


gd) 

(vj) 

= ht 

a^C'!) + 



^t ^1f 

The final form of the expansions (5*28) for the first order 
boundary layer are 

,>]) = + ••• 

and the expansion (5. 50) for the second order transverse curva- 



ture problem are 

) = I'ot 1 ■; 

= ht }y •• 

^1^9= ht + ^ ht . ^5.54) 
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Fow the -various universal f -unctions are as follows. The equa- 
tions, ( 5 . 51 ) give 




= 




( 5 . 55 a) 



= 

(.o---'+Tfo- 

So 


(5.55h) 


Jo(0) = 

0 

= p;{o), 

(-1) = -1 , 


( 5 . 55 c) 


SoCO) = 

0 

= Gq(°°) . 



(5.55d) 

The 

equations 

(5 

. 51’b) leads to 

two sets of 

the universal 

func- 

tions; the first 

is 





Y,(Fl) 

= 

(l-f-nf^) f” + 


( 5 . 56 a) 



=: 

(-.:r-U-^fo+ 

5>?1>8o 


(5.56h) 


F^(0) 


0 = p' (0), 

piC'-'l ) = 


( 5 . 56 c) 


G^(0) ^ 

= 

(°°) = 0 



(5.56d) 

and 

the second 

set is 







= -5 f^p;’ 4 - 2 







+ (1+ ''jf^)f 

1! ft 

u 5nfifo 






-2-'-](2+p^)f 

1 t ■ to 

0^1-1 (-I 

) 

( 5 . 57 a) 




= - 5^10^0 > 







o«i-2') Vi+3'|fig' 

(5.57h) 


1 - 10 ( 0 ) 

== 


) = 0 


(5.57c) 


Qio(o) 

s:" 

0 = 



(5.57d) 
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The above three sets of equations (5.55), (5. 56) and (5.57) 
which involve the solutions of the corresponding first order 
boundary layer equations are integrated numerically on IBM 7044 
computer at IIT Kanpur by Runge-Kutta-Gill method. The procedure 
of solution is disscussed in the Appendix C . The solutions are 
obtained for both the accelerating and the decelerating flows 
(-.198838 < _< 2) and for the values of the Prandtl number 

between 0.7 and 5.0. The results are displayed in Bigs. 5. 7 to 
5.12. 


The transverse curvature effects on various boundary 
layer characteristics (5.8) are still to be studied. Substitut- 
ing (5.9), (5.12) and (5.42) the contribution of the transverse 
curvature may be written as 


■(2 r(2.^ )/(u2r3) = p"( , 0) 


(5.58) 


= k 


^2 


ot ^2,0 ^^^^1t^2,1 ^l^ot '^-2,10^ 


0(5, 0) 


(S^-h^)n-U^r^ 

^2,0 ^^^1t^2,1 h^ot '^2,10^'^ ^5.59) 

00 ^ ^ 

&2 U.,rVf(2^ ) = / k^v^f -B drj 


o 

* 


( 5 - 60 ) 

^ot ^2,0 C^^lt ^2,1 ^1^ot^2,10^'^ 


* 


* 


I I 


Ag U^rVT(2^ ) = I (l-2f )(B -k^ ^-|f )dK| 


( 5 . 61 ) 


^ot'^2,0 ^ ^^^It'^OjI ^l^ot '^2,10^'*'“ ‘ 
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Now substituting the expansions (5.28a) and (5.54) in the 
above and oollecting the coefficients of like powers of ^ , we 

J 

get, 

T2,0 = 


( 5 - 62 ) 


ii 

0 

C\J 

G„(0) 


A* 

^2,0 “ 

i ’1^0 - ^0 

a,, 

^^■2,0 

°° . 
f (1-2-tf ) 

0 ’ ° 

-VP 

the terms ( 

Df order ^ gives two sets; 

■^ 2,1 = 

^"(o) 


f\D 

il 

&^(o) 


^2,1 

t f 

i '^^0 ■ 

d 

^■2,1 

/ (l-2f )(P. 
0 

-•^f^) dr^ 

the second 

set is 


11 

0 

p;'o(o) 


^^2,10 

Olo(0) 



(5.63) 


^ 2,10 { “ ^10 *^1 




2, 10 


/ (1-2 f^) (N^o - ^'jf^ ) - 2f ;(P^- f^) , dr^ 

o. 


The above mentioned boundary layer characteristics (5.62), (5-63), 
( 5 . 64 ) are also evaluated and are shown in Pigs. 5 13 to 5. 15. 
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5.6.2 Longitudinal Curvature Problem: 

Now keeping only those nonhomogene ous terms in (5.2) 

which involve K, we obtain the following problem for u^^^ etc. 

s 

(Superscript 1 will again be omitted here) 

+ (r3v2 + ENr^v.^)^ = 0 

I^(u2)+l2{up-U2j^ = (5.65) 

^ ^2H1T ” 

U2(s, 0) = 0 = V2(s, 0), U2(s, N) = - KN, N - ~ 

^ 2 ( 6 ? 0 ) = 0 = h 2 (B, °o). 

Here, the second order stream function takes the form 


' t 21 f = ^^’^2 

?2S = 

and the '^2 and h 2 are assumed to be of the form 

Y2(s, N) = f(2^) !'( ) 

la^is, N) = (V^l^ ). 

So that, the velocity components are given by 


( 5 . 66 ) 


(5.67) 


U 2 = U.^e’ (5.68) 

^2 " -U^r^ [^F+2^ +(p+jA^-^ 

+ K [f+2f f^ +(p+0 /Vg-D^ f'j 
The differential equations for the longitudinal curva- 


ture problem ( 5 . 65 ) become 
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in ti t I ft 


It f r It ir 


O 

If f 


p f^ -f p._. 
s 

+f P 

dq + / 1-f ^ 

0 

d9 ) 

1 (l-f’^)+2c.^] 


f 1 

V 

+ f f , dr^ 

)] 


(5.69) 


'G + fG ~A.f G+2F (f^ G 

s 

+ /\ P g+2^ (P g__ _P^ g')+kT ]"(- cfV 'lf)g' 

, t 


where 



-2f)'j (f -f^ g ) J 


(5. 

.70) 


0 ) + 2 P^ (f , 0) = 0 




p'(^ 

0 

11 

0 

11 

H 


(5. 

.71) 


0) = 0 = G( > , 00) 


(5. 

.72) 

ki(^ 

) = K(2P )''A/(u^r3) 


(5. 

,73) 

called the longitudinal curvature 

ftmction and 



Ai( 

e ) = (2^ A^) (dk^/d^ ) 


(5. 

>74) 


the principal longitudinal curvature function. The role of the 
latter function A ^ is similar to that of p,A, explained . 


earlier. 


It is interesting to see that for the two dimensional 
flows (3 = G) with s™ o'Ur k^( ^ ) is constant or in other 

words 

■ K cx. s ^ ’ , 

which represents the similarity surfaces of the studies of 
Murphy (1953, 1965) and Narasimha and Ojha (1967). 
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To study the longitudinal curvature prohlem we have to 
assume that P, G and may he represented hy the series similar 
to the transverse curvature series (5*50) and then to split them 
into proper universal functions. Steps employed here are the 
same as in Sec. 5.6.1. So we shall he writing the final form of 
the expansion for P, G and which leads directly to the uni~ 
versal functions 

^ “ ^ol ^o"^ ^^U^I^^I^ol^lO^ + • — 

^ " ^ol ^o"^ ^^^11^1‘^Pl^ol^lO^ ^ (5.75) 

where k^ = k^^ + ,^ k^^+ * • • ■ , ■ 

Substituting the (5.75) and (5.28 ) in (5.69) to (5.72) and collect 
ing the coefficients of like powers of a we obtain the problems 


for Pq, P.| and P^q. The problem for Pq is 

Yo(Jc’ = - ^ 0*0 * Pot (5.76a) 

zys^) = 

Pq( 0) = 0 = P^(0), p'(v^) = , ^ , (5.76c) 

GqCO) = 0 = G^(o°) . (5.76d) 

Here the equation (5.76a) is due to Narasimha and Ojha (1967). 
The problem for P.^ is 

Y^(Fp = V)p(3-f'2)+2a^ (5.77a) 

z.,(op = (-cr"''+'^V3Ppg^ (5.77b) 
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F.,(0) = 0 = FJ'O) F^(>j ) = _Y|, Y| ^ oo, 

(5.77c) 

G.,(0) = 0 y4^(oo). 

(5.77d) 

Finally the problem for F^q is 


Y^(P^o) * -^flV20+P<,)f>;-fX+2f>; 

(5.78a) 

+( '1 f„-i ) V3 f, <-2 V, ( 1+p^) ^ 

t 

^0 

-*0^1+ ’‘Kl-f7)+2(Vfo^i+5flV5ap(p^+1)/(|3_^+3) 
+ [4(Vf„fX)+2-^ (p2HP„+1)]/(Po+1)(Po+3) 

z^(&io) = -3fX - 3(^10- ’|fi)V2(r;-v|f5g^ 


-( cr-"''- 'iV^o^si 

(5.78b) 

P^q(O) = Jio(o) = " ° 

(5.78c) 

G,0(0) = 0 = 5ioM • 

(5.78d) 


The above mentioned (5.76), (5.77) and (5.78), sets of equations 
are solved numerically in a manner similar to tha.t of the trans- 
verse curvature problem of Sec. 5.6.1. The results are shown 
in Figs. 5.16 to 5,21. 


The various boundary layer characteristics (5.8) for the 
ffect of longitudinal curvature are 





= , 0 ) 

^2,0 (^11 ^2,1*Pl*^ot ^2,10 


(5.79) 
)+ — 


^2 ) 


= &'(f , 0) 

12.0 (>^1tl2.l'"Pl’"otl2,1o)*”’ 


( 5 . 80 ) 
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and 


(5.81) 


6* U^rVfC2^) = - / + p’ d/j 

^ol ^2,0^ ^ <^^11 ^2,1"^ Pl^ol ^2,10^"^ •■ 

oo 

= f (l-2f )P v^- f ' d>^ 

(5.82) 

^01 ^ 2,0"^ C, ^2,1'^Pl ^01^2, 10^'^' 


A. 2 U^rVr( 2 ^ ) 


Substituting (5.75) in (5.79) to (5.82) we get 

= <(0) 

= &'(0) 


'* 2,0 


2,0 

* 

2,0 


- f 1 + T„ an 


A 


2,0 = { 


^-2,1 

^2,1 


A 


2,1 


J’" (0) 

G-,^ (O) 

00 

-/ ■''I + P ' dK-| 

0 oo 

1 {i-2i;) r| f; 


■*'2,10'= bo (0) 


q 

6 

A 


2 , 1 o , 
* 

2,10 
2, to 


®1o 


/ bo '^'1 


J " ( 1 - 2 b ) bo . 


(5-83) 


(5.84) 


(5.85) 
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Here ' a = lira {n- f) 

oo ^ 

= “o + + • • 

with 0 = lim ( wi _ f ) 

^ Yj oo ! ^ 

a. = lim ( -f, ) etc. 

‘ y^-^oo ‘ 

The above characteristics are also evaluated for various 
values of and and are displayed in Figs, 5.22 to 5.24. 

5 . 6.3 External Yorticit yi 

In this case we keep the non -homogenous terms in (5.2) 
which involve B^(0) the problem for Ug^^ etc. reduce to 

(r^Ug)^ + ° 

L^jCug) + l2(u.|) - = r^V2(s,o) 

L^jCh^) + L2(h.,)- h2uj^ = 0 (5.86) 

U2(s,o) = 0 = V2(s,o), U2(s, N) = -r^E, N -* «= 

h 2 (s, o) = h 2 (s, °°) = 0 

where ^..(s, o) = lim (v.-Iv.-«) 

N .-00 ‘ liM • 


The equation of continuity in (5.86) is satisfies by introducing 
the stream function 


^^28 = ^^^2 . 


(5.87) 


Further, we assume the following forms for/^ 2 ^4 h 2 

tygCs, K) = 2^ p( c,. 1 )/CTi 
h^(s, N) = G(^,^ ) (h^- sp(2^)’/7uf 
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the velocity components are 


(2f r^p'/U-i 

^ ^ (5.89) 

- ( 2 ( 1 - p ) P+2 ^ + vj ( p+;j 1 ) p ' ) x-O/xj^ 


and the equations (5.86) hecome 


where 


t ! f t t ! t V ^ , n f ! f t 1 f ' 

F +fF-f P +2(1-p)f P+2 ^ (fj. P - fP _f P + f P_ ) 


a (1 +A^) 


(5.90) 


I! ! 


,.t"Tx +fG ~(1 + A^2§)f G+2^ (f G -f G^ ) (5-91) 

= -Pg'+(l + A-2p)p‘g-2ff- (P g*-P’g^ ) 

2(l-p)P( i , o) + 2^ P^ ( ^ , o) = 0 (5.92) 

p’( o) - 0, p'(^Vl ) - --f 1 “ 

a( ^ , o) = 0 = G( 5 , °°). (5.93) 


(5.95) 


as 

aU^ r^ 0 ' 


will he called the principal functi«n of vorticity interaction 
and arises from the second order changes in pressure induced upon 
the boundary layer hy the interaction of the displacement thick- 
ness with vorticity. 

Assuming the following expansions for P and G- 


P = P^ + ^ p.|P.^ + • • •> 

G = & + .f p . G . + 

0 T " 


(5.96) 
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(Si-h^pr^a~' 


i -I ^ ( c > 0 ) 




"12,0 ’ <12,1 ^ 

oo 

= - / ’>^ + e' dv| 

o ’ i 

^ 2,0 ^ ^ ^1 ^ 2,1 ^ ‘ 

OO 

= f (l-2f’ )p'->^f' dn 

O ‘ ‘ . 

= ^ 2,0 + '^ 2,1 + •••. 


( 5 . 100 ) 


( 5 . 101 ) 


( 5.102 


Substituting the expansions (5.96) into above, the teims inde- 
pendent of i' give 


^^2,0 = 





(V) 

V* 

O 

11 

<(0) 




Af- 

OO 

1 

(5. 

.105) 


Sp O = - 
2,0 

- / V| + Jo h 




^2,0 

[ (i-2f’)p;-^,fo irj. 




of order 

; are 



C 

^2,1 = 

f.!' (0) 




12,1 = 

g' (0) 





OO 

(5. 

104) 


^2,1 




0 

A 

^2 1 ~ 

/ (l-2f^)l’,J - 2f’F^-^^f| dvr 
0 





The above boTuidary layer characteristics are studied in the same 
way as the prexious two problems and the results are shown in 
figs. 5.29 to 5.52. 
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5.6,4 Stagnation Enthalpy Gradients : 

Keeping the nonhomogeneus terms in Eq. (5.2) which are 

t 

proportional to S^(o) we have 

(r^Ui)g + (r^v.^)^ = 0 

l^Cug) + L2 (u^) - XL^-^ = 0 /r- 

/ N / s -1 '(5.105) 

^1 ^^2^ fj ^2NN ^ 

U2 (s,o) = V2 (s, o) = U2 (s, 0°) = 0 

h2(s, 0) = 0 , li2(s, N) = "f.^(s, N) ^ , 

How it is to he noted that the momenttun problem is a homogenons 
one v/ith homogenous boundary conditions and its solution is tri-vial 

U2(s, K) = 0, V2(s, N) = 0 (5.106) 

and the energy equation reduces to 

L.^(h2) - if (5.107) 

Assuming the following form for hg 

hgCs, N) = (2q)''/2 G(^ , v|) (5.108) 

the energy equation (5.107) reduces to 

(j“‘^ g'V fG'+(2A-1)f’G + 2^(f^^ g' -f’o„ ) = 0 ('5-109) 

G(0) = 0, G(w|) = Y| - a, '7 

To study (5.109) we assume the following expansion for G 

G( ^ (5.110) 

Substituting (5.110) in (5. 109) and collecting the coefficients 
of various powers of ^ , the terms independent of give 
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^1/2 (5.111) 

Go &o(0) = 0, Go(oo) = v|_ 

and the terms of order ^ give 

^5/2 ^ (5.112) 

G^ (0) = 0 G^ («>)=- . 

These two sets of equations (5.111) and (5.112) are integrated 
numerically and results are shown in fig. 5.33 to Fig, 5.34, 


To study the various boundary layer characteristics (5.8) 
_ * 

2 “ ^2 ■“2 
1 


we note that ^ = 6^ = ^ r, = 0 and 


0 . 2 / (“tr- r''U.j ) = G( 
The expressions for q's are 


,0 


Ip - + 


0) 

«2,1 + 


(5.113) 


<12,0 = h 

42,1 = h 


(5.114) 


The results for (5.114) are shown in Fig. 5.35. 

5*-6.5 Di splacement Speed Problem : 

The remaining non -homogenous terms give the problem for 
u^'^^ etc; 

(rhg)^ + (Dvg)^ = 0 

iqug) + Ig(up - = ( 11 ^ 2)3 (5.115) 

^ 1 (^ 2 ^ ^2^^1 ^ ^ ^2MI ~ ^ 

Ug^s, 0) = 0 = VgCs, 0), UgCsjUgCs, °°)= 

hgCs, 0) = 0 = 12(3, “’). i; ^ 
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III "this Ccis 6 "ttiG SGcond. opd-Gi? siiirG 0,111 fiinciiioii ip. ^ivsn loy 


'''^2F 
i 2s 


= r^u, 


2 


Vr 


(5.116) 


The functions ^2 assumed to he of the form 

■f^(s, N) = f2^^F( ? , Yj) 

^2^®’ N) = (h^-S^ )&( ^ , Yj ) 

and this gives the velocity components 


(5.117) 


u. 


U^F 


(5.118) 


V. 


^2 = - U^r-MF+2? +v](p+3 Ag+1)P )/Y'(2 f ). 

Now the equations (5.115) with the help of (5.117) and (5.118) 
reduce to 

Ml K fit! Mil t f M ■ 

F + P -2pf f +f P+2F (f^ F -f F^ -f^ F +f F^ ) 


= A (2p + A J 


(5. 119) 


r'r^'G'VfG-'- Af'G+2^(f^ G'-f'&^ ) = -Fg' 

+ AF'g + 2^ (p'g^ - F^g') j 

F( , o) + 2e F, (^ , o) = 0 

f'( ^ , O) = o; F'(i^ , oo) = A (^ ) 

G(^ , o) = 0 = G( ~). 

Here the quantity 

^ ^ ^ ^ Ug/U^ (5.120) 

is the ratio of displacement speed to velocity of slip of first 
order outer flow and 
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( 2 ^/\) (dX/d,e ) 

-i / r ds (5.121) 

will "be called the principal displacement speed function. 


= 2dVdx A 


low to study the displacement problem we need to know 
UgCs, o), which may be determined by various methods mentioned 
in the Chapter 1. Suppose that this has been done, so that the 
change in surface speed is 


uys, o') = x” E ;\ x” (5.121) 

n=0 2n 

this expression (5-121) includes the very important cases of: the 
flat plate (m = 0) for which it vanishes ( = o) , the plane 

or axisymmetric flow very near to the point of stagnation. In the 
stagnation region (m = 1) of parabolic cylinder Van I^ke (1964) 
found that A 2 ^ = - 0.61. 

low (5.121) and (5.19) gives 

) = Aq + + ••• (5.122) 

Slid cissuEiing "fchB followiriff ■foriiis fox* P cixid G 


® = Vo • 

the equations (5. '*19) gives 


(5. 123) 


2o(So) = 
ro(o) = 0 
0 ^( 0 ) = 0 


-2P, 


-F g 
0®0 


(5.124) 


J'o(O) , F^(oo) 


= 1 
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Y-j (P^ ) _2(1 + 

Z^(G^) = _3F^g^ 

P^(0) = 0 = P^(0); p;j(co) = -i 
G^(0) = 0 = G^(oo) 


(5.125) 


and 




T f 

o 


Z^CGiq) - - 3fiG' - 3Piog^ + 2 P^g^ 

CtiqCo) = 0 = 


f 


p 

0^1 

(5.126) 


It: is interesting to note that the equations (5. 124), (5.125) 
and (5.126) can he solved analytically in closed form in terms of 
j-irst order boundary layer solutions. The solution for equations 
(5.124) is 


^0 = (fo + T 4>/2 

(5.127a) 

''o = 1 «(/2 

(5.127h) 

for (5. 125) is 

Pl = {4-Po)f,/2 + (f„+3.'i4)A 

(5.128a) 

= (4-Po)gy2 + 

(5.128h) 

and for (5. 126) is 

^10 = 1 f^)/2 

(5.129a) 

^10 1 g-i/2 . 

(5.129h) 

I4e solutions (5. 127) , (5. 128) and (5. 129) are evaluated from 
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corresDonding first order solutions and are displayed in Fig. 5.36 
to Fig. 5.41. Ihe final form of fhe solution to F and G- is giveiiby 
by 

+ ^l(4-PQ)f.^+ A^p^(f^ + qf’)"| /2 + ... 

s r. , , (5.130) 

G = A^r| gy2 + [ A^T) gy2 + (4-p^)g^ 

+ A-o p^vj g' 'j/2 + 

and may be calculated in specific cases from the first order solu- 
tion. In this case various boundary layer characteristics are 




t2f2g/(ll^r''^) = f’'(^,0) (5.131) 

= ^"^0 ^-' 2,0 ^ 2,1 Pl ^ 0 ^ 2 , 10 ^'" ‘ 

G'(r,0) (5.132) 


n''(S^-h„)U^r3 ” 

■ ^0 “Is, 0 ^ *■'^1 ‘^2,1 * '^0 

°° X 

bg U.jrVf(2-?) = / 4 - F dv| (5-133) 


“ ^0 ^ 2,0 ^ ^ ^ 1 ^ 2,1 ^ 1^0 ^ 2 , 10 ^^^ 

f I . f 

U^rVfCSf^) = /(l-2f )F + ^‘f dr] (5.144) 


o 


“'^0^2.0'^ S Pi \^2,1 oA--'. 


Substituting the expansions (5.122) and (5.123) in above and collect 
ing various universal functions we get 
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and 


n,o 

II 

o 

'^2,0 

= <(0) 


C50 

= / 1 - Pq d"'') 

0 

^^^2,0 

oo ^ 

= / (i-2f;) p; + ar| , 

II 


^2,1 = 

G-^ (O) 

c* 

^2,1 

/ 1-P. dv| 

0 t 

A2,1 

{”(1-2^5 h < h 

o 

CM 

= F^o 

^2,10 

= So (0) 

*^2,10 

r°° ' 

= / - P^o 4’) 

0 

^2,10 

oo 

= f (1-2f;) P^o - 2f P„ 

0 


(5. 145) 


(5.146) 


f ^ dy|. 


These characteristics (5.145), (5.156) and (5.147) are also eva- 
luated from the first order boundary layer soltition, the results 
of calculation are displayed In Pigs. 5.42 to 5.44. 
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5 . 7 Discussion 

The higher order terms in the asymptotic expansion for 
large Reynolds number represents the- corre^ to the leading 

(first order boundary layer) approximation. The utility of 
these corrections depend upon the behavior of the first order 
boundary layer itself. In first order boundary layer, the func- 
tion f^ .satisfies the Palkner-Skan equa-tion whose solution has 
been studied extensively in the literature. Solutions to the 
Palkner-Skan eq_uation with exponential decay at infinity exists 
(see, Chapter 4) for -0.198858 < < “ with f^' (O) > 0, However, 

for most of the problems of interest <_ 2. The velocity profile 
! 

f^ vs *>'1 for various values of p^ is shown in Pig. 5.1. The Pig. 

5.5 displays the skin friction ^ (O) vs p^. The quantity 

f (0) decreases as p decreases and becomes zero at p = -0. 198838 
0 0 o 

(this later point will be referred later as the classical separation 
point and cited as CPP). The corresponding energy equation for 
g^ is also studied extensi":’’ely (see, Pai 1956 and Schlichting 1968). 
Some of the solutions for temperature profile g^ vs obtained in 
the present work are shown in Pig. 5.2, Heat transfer rates are 
shown in Pig. 5.6. The magnitude of heat transfer q., decreases 
as p or 1 decreases. 

The solutions for f.| and g.j (the second term of Gortler 
power series in the first order boundary layer theory) are shown 
in Pigs. 5.4 and 5.5 respectively. Hear CPP the functions f^ and 
becomes very large compared to f^ and g^. Purther, the f-unction 



f^(0) g^(0) are shown in Pigs. 5.5 and 5.6, indicate that 

these increases as decreases. The Pig. 5.5 also shows the 
displacement and momentum thicknesses due to similar and non- 
similar terms. Por few cases these equations are studied hy 
Gortler and Sparrow, the comparison of our results with their 
calculations is as follows 


Gortler (1957) f.,( 0 ) 

Present work 

Sparrow (1958) g^ (O) 
Present work 


Po = 0 Po = ^ 

1.032361 0.493840 

1.0323626 0.4938403 

= 1 and (T = 1 
- 0.062118 

- 0.0621187. 


Furthermore, the convergence of the series for f is good 
for positive values of as also shown hy Gortler) hut for nega- 
tive values of and particularly near GPP is poor. The similar 
comments also apply to series for g. It first sight one may expect 
that the expansions for f and g are not uniformly valid near GPP, 

It is shown later that near GPP (i) the convergence of second 
orde'r Gortler Power power series is also p#or (ii) even the con- 
vergence of leading terms (in the first and second order theory 
for = 0 ) is very poor in some of the second order effects. 

Hence there seems a more basic nonuniformity near GPP in inner 
asymptotic expansion itself. The reason is this that our inner 
expansion is not valid near GPP or I any separated flow - as the 
appropriate Euler limit of Navier-Stokes equations is not known 
(see Ghapter 2) . Thus we see that the houndary layer expansion 
Is not uniformly valid in a • involving separation. 
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Now we proceed do sdudy dTie various second order correc- 
tions to the first order boundary layer theory. The solution to 
transverse curvature problem governed by the sets of equations 

( 5 . 76 ), ( 5 . 77 ) and (5.78) are shown in Pig. 5.7 to 5.12. The 

% 

Fig. 5.7 for the solution of the first term F^ in the second 
order Gortler power series show that variation in velocity pro- 
file for various values of are small while in higher order 
terms F.| and F.|^ (see, Figs. 5.9 and 5.11) the variations are 
large, particularly near GPP. Similar remarks applies to the 
corresponding temperature profiles , and G.|q shown in Figs, 

5 . 8 , 5.10 and 5.12 respectively. The skin friction for the three 
gets of equations is shown in Fig. 5.13. As decreases the 
quantities 7 ' „ and 7'o^ increases where T o ^ first increases 
and then decreases to take large negative values near GPP. The 
Fig. 5.14 shows the heat transfer for the three sets. The func- 
tion - q^ , -lo ^ and q^ increases as p decreases. The effect 
d,o 2,1 2,10 0 

of transverse curvature on displacement and momentum thicknesses 
is shown in Fig. 5.15. Finally for a given pressure gradient p^, 
the transverse curvature increases the skin friction and heat 
transfer. 

The problem of longitudinal curvature has been contro- 
vertial (see, Ghapter 1). Let us first study the leading term 
described by equations (5.76). The velocity and temperature 
profile are shown in Figs. 3.15 and 3.16. The momentiun. equation 
( 5 , 76a) has also been studied by Narasimha and Ojha (1967). Our 
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results are in agreement wltli these authors except those of momen- 


turn thickness (see Pig. 5.23). The 

comparison of our 

various authors is as follows 

Source 

' e ■ 

1. 3.^=0 

Murphy (1953) 

1.753 

T an i ( 1 96 4 ) 

1.449 

Yen and T'oha (1961) 

- 0.975 

Hay as i (1963) 

1.449 

Murphy ( 1 96 5 ) 

0.997 

Schulta-Grunow and Breuer 

(1965) 1.442 

Narasimha and Ojha ( 1967 ) 

1.4469 

Present work 

1.4469667 

2- Po = ^ 

V an Dyke ( 1 96 2h ) 

1.913255 

larasimha and Ojha (1967) 

■ 1.9132 

Present work 

1.9132549 

The comparison of heat transfer results for -J* = 0.7 is 

Source 


Po = ^ 

i 

Schultz-Grunow and Breuer 

(1965) 0.245 

Present work 

0.2508955 

2. H = 1 

Yan Dyke ( 196 2h) 

0. 12811 

Present work 

0. 1281108 
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It IS to be noted that the velocity and temperature profiles 
sho» in Pigs. 5.15 and 5.16 become very large near OPP. fhe ve- 
locity and temperature profiles for nonslmllar sets of equations 
(5.77) and (5.78) are shown In Pigs. 5.17 to 5.20 and these also 
becomes large near OPP. ihe skin friction for the three sets 

(5.76) to (5.78) is shown in. Pig, 5 22 The fin- r 00 1 . 

° j.dc. me iig 5.22 shows that 

12 , 0 ’ ’^ 2,1 ®“d increases as decreases. The similar 

trend is observed for heat transfer results displayed In Pig. 

5.23, when Prandtl number is prescribed. If fa held fixed the 

heat transfer Increases as Prandtl number increases. The Pig. (5.24) 

shows the effects of longitudinal curvature on displacement and 

momentum thicknesses, lastly, regarding the effect of longitudinal 

curvature, it is clear thet convex curvature reduces both the skin 

fi?ic"bion “th.8 hocit; "tiTQ-iisfBir, ' 

The subject the effect of vortlclty In the basic flow upon 
the classical boundary layer has been the subject of prolonged 
controversy, (see. Chapter 1 ). let us first study our first, 
term in Sortler power series expansion for the vortlclty effects 
described by equations (5.97). The velocity and temperature pro- 
files for these equations are shown In Pigs. 5.25 and 5.26. The 
skin friction and heat transfer are shown In Pigs. 5.29 and 5.30. 

The oomparison with various authors is an follows: 

Skin friction: 

Source 


1 . Po = 0 


Van Dyke ( 1962 b) 


-<( 0 ) 


3.1260 
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* 

Glauert (1957) 


0. 795 

Murray ( 1 96 1 ) 


3 . 1259 

Present work 


3.125983 

2. p = 0.5* 

Rott and lenard (1959) 


0.649 

Kemp ( 1959 ) 


0.650 

Van Dyke (196.2b) 


1.76861 

Present work 


1. 76 8 5 9 7 

3. P. = 1.Q 

Stuart (1959) 


0,6078 

Van Dyke (1962b) 


1.40652 

Present work 


1.406543 

Heat transfer: 

Source 

c-yo) 



TD 

0 

II 

0 

0 

11 

0 

CO- 

1. •:)■= 0.7 

Tan Dyke ( 196 2b) 

0.91117 

0.38270 

Present work 

0.911168 

0.38269 


2. cr= 1.0 

*0vchinnikov (1960) 0.340 

Present work 1.04 26 45. 


The authors marked with above have neglected the induced pres- 
sure gradient by the interaction of external vorticity with classical 
boundary layer. So in some sense they have solved -the momentum 
equation (5,97a) with omitted. Thus if a particular integral 
IS added to the results of above authors the present results can be 
obtained. The Pig. 5.30 displays the heat transfer 


for 
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various values of As decreases ^ increases. Further, 

the velocity and temperature profiles for the second term in 
Crortler power series described by (5.98) are shown in Fig. 5. 27 
and 5.28. The skin friction and heat transfer results shown in 
Figs. 5.29 and 5.30. The effect of vorticity on displacement 
and momentum thicknesses is shown in Fig. (5.31). Finally the 
effect of vorticity is to decrease skin friction for all values 
of pQ, Further, the external vorticity also reduces the heat 
transfer except for p^ > 1 where it increases. Also near GPP the 
velocity and temperature becomes large. 

The effect of stagnation enthalpy on the boundary layer 

is described by equations (5.111) and (5.112). The Figs. 5.32 

and 5.33 shows the tem.perature profiles for the above equations. 

The heat transfer is shown in Fig. 5.34, which shows that as 

decreases the function qr, ^ decreases while q„ ^ increases. The 

v , 0 ^2,1 

effect of stagnation enthalpy gradients in basic flow is to reduce 
the heat transfer. However near GPP the temperature becomes very 
large. 

Lastly for the displacement problem (Sec. 5.6.5) the exact 
close form solutions ais given in terms of the first order bound- 
ary layer. The solutions (5.127), (5.128), (5.129) for the sets 
of equations, (5.124), ( 5 . 125 ) and (5.126) are shown in Figs. 

5.36 to 5.41. The veloeity and temperatures again increases near 
GPP. The skin friction is shown in Fig. 5. 42. The figure shows 
^2 o i^o <^®®neases while ^^2 1 "^ 2,10 


increases. 



142 


oimilar comments apply to the heat transfer results shown in Fig. 

5 . 43 . The Fig. 5.44 shows the effect of displacement speed on the 
displacement and iDomentum thicknesses. The effect of disnlacement 
speed is to increase the skin friction and heat transfer. 

5 . 8 Conclusions t 

The present solution for the second order effects show that 
the transverse curvature increases the skin friction and heat trans- 
fer. The convex longitudinal surface curvature decreases the skin .t , 
friction and heat transfer. The effect of external vorticity is to 
decrease the skin friction. The heat transfer is also decreased 
except for > 1 where it increases. The affect of stagnation ©n 
enthalpy is to decrease heat transfer*. Finally, the displacement 
speed increases the skin friction and heat transfer. 

Lastly, the results indicate that the convergence of Gortler 
power series for the first and the second order boundary layer 
functions becomes poor as favourable pressure gradient diminshes. 
Furthermore, in adverse pressure gradient (particularly near 
Pq = - 0.198838) some of the second order boundary layer quantities 
become^' very large compared with the corresponding first order boundary 
boundary layer quantities. This indicates that the boundary layer 
(inner) expansion is not uniformly valid near the separation. 
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APEENDIX L 


EYi'iUXTION OP INTEGRi'JjS 


In Chepers 3 and 4 we have used the following relations 


1. I. 


v/he re 


/ t^ exp(-at^‘/r! ) dt 


J. Or! ^ r rn+1 


r ' r 


^ (A.1) 


x) 


; e-’t t“-1 at 


is the incomplete gamma function. 

To show the equraity of this relation (1.1) we introduce 
new variable y(t) = at^/r! in I. to obtain 


n+1 / t n-fl-r 

T 1 /r! N r ^ r ”'-y , 

II = ^ (“) I y e dy 

n+1 ° ^ 

= i (f) , # -"). 

oo 2; . 

2. = f 2 ^ exp (-Pz) f t^exp(Pt-at)dt dz 


(A. 2) 


j'(m+n+2) Ymj.nj.o'l 

= (in+ 1 ^ m+n+2, n+2, 

where (a,b,c,d) is the well-known hype rgeome trie function 


defined by 


■ ■ OO . 

(a,b/c,x) = 

= ITa+n)/ IT?-) and [la) - (a, is the comp- 


where ~ I |(a) and }(«) = d («? '”) is ibe comp- 

lete garnmr, function. 

To prove the relation in (1.2) we first note that the integral 
in curly .brackets is rela.ted to incomplete gamma function and with 
the help of (l.'l) can be written as 
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/ (P't-S-’fc) dt - (g-P) ^ ^ (n+1 , az-Pz) 

o 

Peking use of the relation (see, Irdelyi 1954) 


A ) 

f x°“ e“f^^ ^(d, ax) dx 

0 


Rli+d) 

d(a+p)''+'^ 


2F^(1,c+d,d+1, 


(val.id when real part of a+j3, p, c+d are positive) the integral 
in (a. 2) ‘become 

h = 2h (1- “™+2’ ”+2’ l-P/®) 


further using the recurrence relation (see, Sneddon 1961) 

2^-1 (®^> Ps d, 1 — X ) — X ^ 2^1 1 — l/x) 

the above integral Ig become 

J . Ig - tg+2). m+n+2, n+2, 1-a/P) p-(^+^+2) 

X z 

3. = / z^ exp(^azV^! ) / exD (atV^! ) dz (A. 3) 

^0 0 

ip-m+ 2 

rT^TJ ^ r ’ r! 

r 

X 

where N^(a,x)= / t® I[(a,t) dt 

0 

-| 

ll(a, x)= a / e''^'^( 1-t)^“"* dt. 

0 

Here H(a, x) is a certain related gamma function studied by 
Tricomi (1950) and Haras imha (1964), 

The integral with in curly brackets is related to incomplete 
gamma function with a negative argument (See, Tricomi 1959) and 
using a procedure similar to in (i,.1) we get 
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z ^1 

/ exp (e.tVr! ) dr = ^ (r!/a) ^ '5^ (— 

n I I* 


az 

-r 


') (A. 4) 


wtiere. 


) ^ t a 1 

^^(a, x) = / e t dt. This f-unction is more conve- 

0 

nientl}/ studied in terms a related function N(a, ±) defined hy 




Oi(a» x) 


1 ,,a 
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x“ y(a, x) 


dt. 


or N(a, x) = ae~ x“ ^.j(a,x.) = a / 
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By means of (A. 4) and (1.5) the integral can he written 

. (Tl s 

r(n+1 ) ^ a. ^ 


( 1 . 5 ) 
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Some of the properties of the function N( a, x) are (see, . Tricomi 


1950) 


0 < N (a, x) < 1, U(a, 0) = 1, 


F(a, oo) = 0, N(0, x) 


-X 


further, it is easily shown (see, Narasimha 1964), that the asymp- 
totic expansions of N(a, x) for x - 0 with a fixed is 

2 


H - 1 

and for x — oo with a fixed 


X _r’ 

a+1 (a+iyra+^ 


g 

X 


_ a(a-l ) + «( ..«-. 1 I LoL-^ _ ... . • 
x 5 



APPEITDIX B 


EUIZ.RISATION; (E. a) METHOD 


To improve the convergence of a slowly convergent or even 
divergent series various methods are used in the literature which 
accelerate the convergence. Several interesting examples of such 
transformations are given hy Van Dyke (1964). One of the most 
widely used Being the process called Eulerisation (See, Meksyn 1961), 
It is found in the present work, however, that although the con- 
ventional Eulerisation procedure was not always successful, an 
extension described hy Hardy (1949) and called hy him the (E, q). 
process gave satisfactory results. 

Following Hardy (1949) to improve the convergence of a 

series 

= S a^ x”^’^ (B-1) 

^ n=0 


which is convergent for sufficiently small values of x, hy the 
(E, q) method, we proceed as follows. 

Let X = y/(1-qy) or y = x/(1+qx), q 0, (B-2) 


It is easily seen that the substitution of (B-2) in (B-1) gives 

m+1 

(B-3) 


S„ = P (y(l+q)) 

m=0 


where 


Q- 


(^1) = 


m 


n 


(q+1)®+'‘ «=G 


m-B _ 
IX' " 


2 L) 1" '■ a. 


If we choose X ^ 1 , then (1+q)y = 1 > and the series (B-4) reduces to 


= 2 a 

m=0 


(l) 


(B-4) 
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The series (B— 4) is said "to he ihe Eulerised form (hy q_) of 
n+7 

S = S a x . To Eulerise a given series, this q. has to 
n=0 

he determined hy a ftrial and error procedure (in the conventional 
Eulerization procedure q is taken as unity). Eurther, in practice 
it is found convenient to Eulerise a series leaving out some of the 
leading terms of the expansion, as Eulerization including in par- 
ticular (E,q) method, leads to unique results (See, Handy 1949), 
the final result does not depend in principle on the number of 
terms left out. It may he noted that it can he shows (see, Hardy 
1949), that if a series is summahle (E, q’) and q > q' , then it 
is also summahle (E, q) to the same sum. 

To plot the Eulerised results of lower series (4.45) for 
heat transfer in Eig, 4.4 and 4.5 we fotuad that Eulerisation with 
the following values of q gave rapid convergence,. 

1. D(-0, 198838, 0“ ) = 0.79788 .0501 7 'y 

+ 3.24851 - 6. 93062 + 0(cr^^^) 

cr~ q Leading terms left 

0.05 1 

0.1 1 

0.3 1 

0.6 2 

2. ,D(-0.16,0") = 0.79788 f?“V2 _ p. 086 49 CT 

+ 1.85581cr^/^ - 5,07590cr^+0(o-^5/^) 



0. 1 
0.5 

1.0 


1 

1 

2 


Eirst only 
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3, D(0,cT-) = 0.79788cr^/2 _ q. 77463 ■:?“+ 1.032l2cr^/^ 

- 1.31774<r^ + 0(<y-5/2) 

tr q. Leading terms left 

0.3 1 

0.5 1 

0.8 1 

1*0 1 ^ First only 

4.0 3 

9.0 4 

16.0 9 J. 

Further, to plot the Eulerised low ct result (4.4^) for the 
recovery factor in Fig. (4.6) the following values of q gave 
rapid convergence, 

Tp=0,<r) = 0.9255 + 0.1951 cr - 0.166lcr^/2 

+ 0.02%cr^ + o(cr^^^) 


0.9 

“I 1 



1.0 

1 f 

First 

two only 

4.0 

2 i 



9.0 

3 

First 

only , 


It is well known that the succes'^ of the usual Euleriza- 
tion procedure (q = 1) rests on the fact that the function re- 
presented by the series (say in x) possesses a simple pole singu- 
larity in the complex x - plane near x = -1 (See, Van Dyke 1964). 
Prom (B-2) it appears that (E,q) method is an extension to allow 
for the presence of the pole near x = -1/q. Ihe fact that in 
present work the appropriate q to use depends itself on x is 
perhaps due to the presence of more complicated singularities 
than simple poles in complex plane. Ihe main point is that 
one can determine a value of q which gives satisfactory conver- 
gence, so that the final results so obtained from a series 
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nominally constructed for low<rare excellent forCPso high. that' 
they agree with results obtained "by considering the opposite 
limit of large CT" » 



iiPEENDIX C 


FUIERIGAL SOLUTION OF HIGHER ORDER BOUITDARY LAY ER 

PROBLEMS 


A general programme to solve, any finite numler of total 
differential equations of fhe type studied in the present work, 
numerically hy Runge-Kutta method with Gill improvement is deve- 
loped for IBM 7044 Computer at Indian Institute of Technology, 
Kanpur. To illustrate the method of solution, let us consider 
the governing total differential equations for any one of the 
second order effects, say, transverse curvature. The problem 
of transverse curvature is defined by the sets of equations (5.55), 
(5.56) and (5.57) in addition to the corresponding first order 
boundary layer equations ( 5 - 25 ) and ( 5 . 29 ). Each of individual 
equation of above mentioned sets constitutes the two point bound- 
ary val-ue problem; since for example, for velocity profile two 
of the boundary conditions are known at = o, the third known 
at Y"! = 00 ^ To solve each of such equations it is necessary to 
guess a condition missing at the boundary = 0 and integrate 
towards = oo'to see if the boundary condition is meet here. The 
problem for f^ is nonlinear and therefore some what trouble some 
to handle. However, the solution to this equation is given by 

■ ■ . ft 

Tar ious aut Hors mertioBed iB OHapter 4 , tHe missing value f^(0) 
is taken from Smith (1954) in the present work. Rest of the 
equations are linear two point boundary value type and may be 
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integrated easily by some sort of iteration or interpolation process. 
Hov/ever, to solve each of these actuations separately (after obtain- 
ing the solutions to previous ones) is lonnecessar'^.r complicated and 
time consuming. We here have adopted a more efficient technique 
by solving at a time as many equations a.s possible. Prom the 
following chart one may find such number of equations which can 
bcC' solyed at a time 
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JroB the Chart It la cbcrlcus that once the solution for Is 
imown (for example by generating the profile from the InlMal 
values of tjo) obtained from Smith) we can solve simultaneously 
three equations f ^ at a time in the first Instance. In 
second we can solve g., , and In the third S., , G.,q, 

To accomplish the above, the programme is to be engineereci so 
as to study either in three different loops or all equations in 
a loop. The later one is more efficient and useful. To write 
up such a programme which can handle all equations at a time it 
was found convenient that the, equations be first expressed as a 
set of first order differential equations. The each set of equa- 
tions is therefore reduced to a set of five first order equations 
as follows, 

fo ; x( 2 ,j) = f^, X(3, J) = f’, x(4, J) = f^ 

X(2,0) = 0 = X(3,0), X(3, “) = 1 
: X(5,J) = g^, X(6,J) = g' 

X(5,0) = 1, X(5, ») = 0 
: X(7,J) = f^, X(8 J) = f^\ X(9.J) = f *' 

X(7,J) = X(8, 0) = X(7, oo) = 0 
V X(10,J) = g^, X(11, J) = g| 

X(10,0) = G = X(10, '»). 

xci2,j) = i- ; x(i3, J) J) = % 

X(12, 0) = 0 = X(13,0), X(t3, °°) = X(1 


: X(15,J) = G^, X(16,J) = g' 

X(15,0) = 0 = X(15, ~) 

: X(17,J) = P^, X(18, J) = p^, X(19, J) = p’ 

X(17,0) = 0 = X(18, 0), ,X(18, oo).= x(1, oo) 

G^ X(20,J) = G^, X(21, J) = g' 

X(20, 0) = 0 = X(20, °°) 

P^Q : X(22,J) = P^Q, X(23,J) = P^q, X(24, J) = P^'q 

X(22, 0) = X(23, 0) = X(23, °°) = 0 

G^q X(25,J) = G^q, X(26, j) = g’q 
X( 25, 0) = 0 = X(26, ~), 
where X( 1 , J ) = . 

A programme was developed so as to handle all 25 first order 
eq_uations in one loop. Assuming suitably the missing tboundary 
conditions, the equations are integrated numerically by Punge-Kutta- 
Gill method from = 0 to ^ ^ (where 7^ ^ was chosen to be 

between 5 to 9 depending upon and {J“ ) with a step size^’|= 0.05. 
After completion of first integration, the second integration is 
performed by giving a small perturbation to all the guessed quan- 
tities in the first integration. The two solutions were interpolated 
linearly so as to satisfy the boundary condition at The 

procedure of integration is repeated by setting up an iterative- 
loop. The initial values for each integration were corrected with 
the newly determined final values at the infinity. It was found 
that after few trials the procedure converged to the correct boundary 
conditions for the entire system. 
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Fig. 4.2_Function D( p,cr-) Vs ^ for low Prandtl number flows 





Eulerised result 



Fiq.4.4_Comparison 
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Fiq.5.1- Basic velocity profiles (Eqn. 5.23a) 





boundary layer ( Eqn. 5.24a ) 
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Fiq.5.7_ Change in velocity profiles due to transverse 
curvature (Eqn.5.55a) 
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x8_Chanqe in temperature profiles due to transverse 
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Fig.5.9_Change in ve^locity profiles du^ to tronsvc^rs^^ 
c.ury.ature_(EQn.^5^^^ 
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Fig.5.10_Chanqe in temperature profiles due to transverse 
curvature {Eqn.5. 56 b) 



11_Chan9e in velocity profiles due to transverse curvature 
(Eqn. 5.57a) 
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Fiq.5.12_Change in temperature profiles due to 








'.-ir 




Fig 5.13- Change in skin friction due to transverse curvature. 
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Fig. 5.17_Change in temperature profiles due to longitudinal 
curvature (Eqn. 5.76 b) 
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Fiq.540_ Change in velocity profile due to 

( Eqn. 5. 




20_Chanqe in velocity profiles due to longitudinal curvature 
(Eqn.5.70a) 











Chanqe in vdocity profile due to 
external vorticity (Eqn. 5.97a) 















Fig. 5.32_Chanqe in displacement & momentum thicknesses due to 
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Fiq.3.40_Change in velocity profiles due to displacement 
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Fig.5. 44- Change in displacement and momentum thickness du 




